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COMMUNICATION  SYSTEM  PLi EORMANCE  MODEL 
FOR  VHF  AND  HIGHER  FREQUENCIES 


R.D.  Jennings  and  S.J.  Paulson* 

Conmiunication  system  designers  and  engineers  often  must 
decide  whether  to  use  design-aid  models  not  well  suited  to  the 
problem  at  hand  or  to  perform  many  tedious  and  laborious  hand 
calculations  which,  coupled  with  certain  "rules  of  thumb,"  provide 
an  estimate  of  system  performance.  The  model  documented  in  this 
report  provides  relief  from  such  problems.  For  communication 
systems  which  operate  at  VHF  and  higher  frequencies,  the  character 
of  the  terrain  between  the  transmitting  and  receiving  antennas,  as 
well  as  the  heights  of  the  antennas  above  immediate  terrain,  will 
affect  significantly  the  attenuation  of  the  radio  signals.  The 
model  reported  herein  considers  these  terrain  influences,  along 
with  other  commonly  considered  influences,  by  automatically  accessing 
digitized  topographic  data  files  to  develop  the  path  data  needed 
to  compute  the  basic  transmission  loss.  Computation  of  the  basic 
transmission  loss  is  accomplished  using  the  Longley-Rice  formulation 
for  point-to-point  paths.  The  statistical  character  of  the  basic 
transmission  loss  results  from  long-term  fading  (time  availability) , 
path  to  path  differences  (location  variability),  and  prediction 
confidence  considered  in  the  model.  The  type  of  model  output  data 
is  selected  by  the  user.  Options  range  from  simple  path  terrain 
information  to  basic  transmission  losses  to  utilizations  of  the 
basic  transmission  loss  in  computing  power  density  or  received 
signal  level.  There  is,  in  addition,  an  output  option  which  portrays 
predicted  communication  reliability  as  the  probability  that  received 
signal  level  will  exceed  a specified  threshold  for  specified  condi- 
tions of  time  availcibility  and  location  variability.  All  output 
options  will  provide  data  in  either  tabular  or  plotted  form.  Plotted 
output  data  are  contoured  over  a geographical  area  not  to  exceed 
two  degrees  latitude  by  two  degrees  longitude. 

Key  Words:  Basic  transmission  loss;  communication;  communi- 
cation reliability;  communication  system  per- 
formance; con5)uter  model;  received  signal  level; 
path  profile;  power  density;  propagation; 
propagation  loss. 


*The  authors  are  with  the  Institute  for  Telecommunication  Sciences,  Office 
of  Telecommunications,  U.S.  Department  of  Commerce,  Boulder,  Colorado. 
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The  U.S.  Army  Communication  Command  (USACC) , with  general,  world-wide 
responsibilities  for  Army  communications,  includes  the  U.S.  Army  Communi-  ; 

cations  Electronics  Engineering  Installation  Agency  (CEEIA) . The  CEEIA 
has  the  technical  responsibilities  for  designing  and  installing  the  Army's 
strategic  and  tactical  communications,  so  that  required  system  performance 
criteria  are  satisfied.  Sophisticated  techniques  are  available  to  aid  j 

in  the  design  of  fixed-point- to-fixed-point  communication  links.  Such  j 

aids  usually  use  a propagation  loss  model  generally  referred  to  as  a "point- 
to-point"  loss  model.  Aids  also  are  available  for  predicting  the  performance  j 

of  land  mobile  and  broadcast  types  of  radio  communications.  By  contrast,  j 

these  aids  usually  employ  a propagation  loss  model  called  an  "area  prediction"  ] 

loss  model.  The  area  prediction  models  use  average  values  to  describe  | 

the  propagation  path  rather  than  specific  values  as  is  done  in  a point- 

to-point  loss  model.  Other  frequently  used  design  aids  involve  tedious  i 

and  laborious  hand  calculations  and/or  employing  "rules  of  thumb"  to  estimate 
the  performance  of  a communication  system. 

Application  of  the  more  specific  techniques  of  a point-to-point  propaga- 
tion loss  model  to  a large  number  of  paths  in  order  to  predict  the  performance 
of  mobile  and/or  broadcast  types  of  radio  systems  is,  in  essence,  the  approach 
of  the  automated  model  documented  by  this  report.  The  situation  can  be 
described  as  fixed-point-to-mobile-point  communications.  The  application 
of  a point-to-point  loss  model  to  such  situations  involves  the  definition 
of  a large  array  of  points  within  a defined  geographical  area.  Each  point 
linked  with  the  fixed-point  location  is  treated  as  a point-to-point  path 
for  purposes  of  estimating  the  propagation  loss  between  those  points. 

I 

The  CEEIA  has  identified  a need  for  this  type  of  communication  system  design  i 

aid  to  provide  evaluation  of  communication  system  performance  in  such  terms  ! 

as  received  signal  level  or  communication  reliability  (defined  in  Appendix 
F) . This  capability  is  defined  and  discussed  in  this  report.  i 

The  next  section  of  the  report  presents  a discussion  of  typical  applica- 
tion for  the  model,  section  3 is  a general  description  of  the  model,  section 
4 lists  a summary  of  important  model  features  and  output  data  options, 
and  section  5 discusses  recommendations  for  improvement  of  the  model. 

The  detailed  User’s  Guide  for  the  model  is  Appendix  A,  and  Appendix  B is  the 
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Software  Technical  Description.  Technical  details  of  the  point-to-point 
propagation  loss  model;  the  path  profile  and  path  parameter  calculation 
routines;  power  density,  received  signal  level,  and  communication  reliability 
calculations;  and  other  aspects  of  the  model  are  discussed  in  the  remaining 
appendices . 


i 2.  TYPICAL  APPLICATION 

In  a more  specific  context  than  the  preceding  introduction  has  revealed, 
the  Communication  System  Performance  Model  (CSPM)  has  been  developed  to 
provide  Army  communication  system  designers  and  engineers  an  ability  to 
predict  the  expected  performance  of  the  line-of-sight  (LOS)  and  tropo- 
spheric scatter  communication  systems  to  be  installed  in  support  of  the 
Tri-Service  Tactical  (TRITAC)  communication  system  and  various  other  local 
Ft.  Huachuca  communication  systems.  These  specific  applications  have  not 
dictated  any  model  limitations,  however. 

Suppose,  for  example,  that  the  communication  performance  of  mobile 
and  semi-mobile  units  linked  with  a well-sited,  fixed  installation  is  to 
be  evaluated.  An  engineer  performing  this  evaluation  will  know  the  transmitter 
output  power  and  the  antenna  heights  above  ground  and  gains  for  the  fixed 
installation  and  the  mobile  and  semi-mobile  units.  Furthermore,  he  will 
have  determined  that  satisfactory  system  performance  will  be  realized  when 
the  received  signal  level  (that  is,  the  signal  level  delivered  by  the  receiving 
antenna  to  the  input  port  of  the  receiver)  is  greater  than  some  threshold 
value  for  the  received  signal  level. 

Since  the  problem  involves  mobile  and  serai-mobile  systems,  the  evalu- 
ation must  consider  communication  between  any  point  within  some  geographical 
area  and  the  fixed-installation  location.  The  estimates  of  propagation 
losses,  then,  could  be  calculated  using  average  values  to  describe  certain 
essential  terrain  characteristics  required  by  the  propagation  loss  model. 

However,  suppose  the  engineer  decides  that  he  needs  better-defined 
estimates  than  can  be  calculated  using  average  characteristics.  He  decides 
that  estimates  of  propagation  loss  calculated  by  using  a point-to-point 
propagation  loss  model  will  be  satisfactory,  since  the  model  for  such  estimates 
' uses  specific  characteristics  rather  than  average  characteristics  for  each 

propagation  path.  With  sufficient  numbers  of  such  calculations  for  specific 
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point-to-point  propagation  paths  throughout  the  geographical  area,  the  com- 
munication system  performance  from  any  point  in  the  area  linked  with  the 
fixed  point  is  approximated  very  well.  The  CSPM  provides  this  type  of 
information. 

First,  however,  one  must  realize  that  propagation  loss  is  a variable 
phenomenon,  even  for  a defined,  specific  path,  due  to  meteorological 
variations.  Furthermore,  there  will  be  additional  variation  from  path  to 
path,  due  in  part  to  the  fact  that  different  paths  can  be  described  with 
the  same  model  input  path  descriptor  values.  Finally,  there  is  uncertainty 
associated  with  the  prediction  because  finite  sampling  has  supported  the 
development  of  the  prediction  model,  cuid  influences  which  do  exist  are  not 
modeled  explicitly. 

The  CSPM,  then,  provides  statistical  estimates  of  basic  transmission 
loss  for  a large  number  of  propagation  paths  throughout  a geographical  area. 
The  basic  transmission  loss  estimates  are  used  to  calculate  statistical 
estimates  of  power  density  and  received  signal  level.  Also,  the  statistical 
nature  of  basic  transmission  loss  is  utilized  to  provide  an  output  describ- 
ing communication  reliability  throughout  the  geographical  area  for  specified 
time  and  location  variabilities.  An  incidental  dividend  of  the  model  is 
the  option  for  output  data  which  describe  the  terrain  profile  of  a propagation 
path. 

In  order  to  use  the  CSPM  to  evaluate  the  problem  outlined  earlier, 
the  engineer  must  define  the  statistical  levels  at  which  the  satisfactory 
performance  criterion  must  be  met.  Suppose,  then,  that  received  signal 
level  (as  a function  of  specified  time  availability,  q^,  and  location  vari- 
ability, q ) greater  than  the  threshold  received  signal  level  will  provide 
L 

satisfactory  service.  The  communication  reliability  calculation  capability 
of  the  CSPM  can  be  used  to  provide  output  data  which  show  the  probability — 
comnunication  reliability — with  which  these  conditions  Ceui  be  expected 
throughout  the  geographical  area.  In  addition,  other  values  for  time  avail- 
ability and  location  variability  can  be  used  to  explore  the  impact  of  trade- 
offs that  may  be  necessary  to  consider  such  as  reducing  the  value  for  location 
variability  while  maintaining  the  value  for  time  availability.  The  statis- 
tical concepts  of  time  avail£d}ility  and  location  variability  are  discussed 
in  Appendix  C. 
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These  model  outputs  for  a hypothetical  problem  are  shown  in  Figures 
2-1  and  2-2.  This  type  of  data  is  typical  for  the  model  being  used  in  its 
"most  powerful"  form,  so  to  speak.  Several  "reduced-power"  applications 
also  are  available. 

For  example,  the  same  problem  may  be  considered  without  knowing,  or 
at  least  without  defining,  a threshold  for  received  signal  level.  That 
is,  no  criterion  for  system  performance  evaluation  is  known  or  defined. 

Now,  the  model  could  be  asked  to  calculate  and  plot  up  to  ten  equal-value 
contours  of  received  signal  level  throughout  the  geographical  area  of 
interest.  These  plots  would  be  provided  for  specified  values  for  the  vari- 
abilities— time  availability,  location  variability,  emd  prediction  confidence. 
Such  an  output  is  illustrated  by  Figure  2-3. 

A further  reduction  in  defining  the  input  data  for  the  model  can  be 
accommodated  by  asking  the  model  to  provide  equal-value  contours  of  power 
density.  These  data  may  be  obtained  without  knowing  (or  specifying)  the 
receiving  antenna  gain.  Or  perhaps  in  the  early  stages  of  defining  a problem, 
neither  the  transmitting  nor  the  receiving  antezina  gain  is  known.  Then 
the  model  user  may  be  interested  in  contours  of  basic  tr^msmission  loss. 

For  such  data,  the  transmitter  output  power  also  would  not  need  to  be  known 
or  specified. 

If  contours  of  power  density  or  basic  transmission  loss  were  produced 
for  the  same  problem  used  to  generate  the  received  signal  level  contours 
shown  in  Figure  2-3,  the  contours  would  have  very  similar  shapes.  This 
would  occur  because  the  calculated  values  for  received  signal  level,  power 
density,  cmd  basic  transmission  loss  at  each  point  have  constant  differences. 
Power  density  is  different  from  received  signal  level  by  a constant  value 
which  accounts  for  the  receiving  antenna  gain,  appropriate  aperture-to- 
medium  coupling  loss  adjustments,  and  the  free-space  propagation  loss. 

Basic  transmission  loss  is  different  from  received  signal  level  by  a constant 
value  equaling  the  sum  of  the  transmitter  output  power,  the  transmitting 
and  receiving  antenna  gains,  and  the  appropriate  aperture- to-raedium  coupling 
loss  adjustments.  The  concepts  cf  communication  reliability,  receivT-il  signal 
level,  power  density,  and  basic  transmission  loss  are  discussed  In  Appendix  P. 

Any  of  these  types  of  ou^ut  data  may  be  requested  by  a user  If  he 
specifies  only  the  trzuismltter  coordinates  (euid  receiver  coordinates,  if 
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Figure  2-3.  Received  signal  level  throughout  a geographical  area  for 
conditions  of  time  availedsility  = 0.90  and  location 
variability  = 0.90,  and  confidence  = 0.50. 
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a path  profile  is  requested) , since  the  model  has  built-in  default  values 
for  all  other  input  data.  Of  course,  the  use  of  the  model  must  evaluate 
the  output  data,  making  some  subjective  judgment  about  how  well  the  default 
conditions  represent  his  problem.  In  addition  to  the  plotted  data,  each 
type  of  output  data  may  be  requested  in  a tabulated  form  as  well. 


3 . MODEL  DESCRIPTION 

Though  designed  for  machine  independence,  a CDC-6500  system  will  be 
employed  by  the  Army  for  applying  the  model  to  communication  system  perform- 
ance evaluations.  The  model  has  been  designed  to  provide  flexibility  for 
the  user  in  his  study  of  existing  or  proposed  communication  systems.  The 
model  user  may  define  a single  communication  path  or  a geographical  area 
throughout  which  communication  is  desired.  The  model  will  extract  terrain 
data  from  a digitized  topographic  data  file,  use  the  extracted  data  to  con- 
struct path  profiles,  and  calculate  appropriate  path  parameter  data.  The 
path  parameter  data  are  required  by  the  point-to-point  propagation  loss 
model  which  calculates  estimates  of  basic  transmission  loss  (I^)  which  may 

be  used  in  subsequent  calculations  of  power  density  (P ,) , received  signal 

d 

level  (RSL) , and/or  communication  reliability  (REL) . The  user  selected 
outputs  of  Lj^,  P^,  RSL,  and/or  REL,  as  well  as  path  profile  data,  are 
provided  in  tabular  and/or  plotted  form.  When  communication  reliability 
is  specified  as  a desired  output,  a received  signal  level  threshold 
for  system  performance  also  must  be  specified.  Detailed  technical  description 
of  the  model  is  provided  in  Appendix  B,  Software  Technical  Description, 
and  Appendices  C,  D,  E,  and  F. 

3.1  Model  (Program)  Organization 

The  model  is  divided  into  four  program  overlays,  so  that  reasonable 
requirements  for  computer  core  are  not  exceeded.  The  organization  is  illus- 
trated in  Figure  3-1.  The  primary  (0,0)  overlay,  with  program  CSPM  as  the 
overlay  controller,  provides  overall  program  control.  The  source  code  for 
program  CSPM  contains  a commentary  on  the  model  and  the  coomon  block  variable 
definitions  for  /OPTIONS/  and  /LOSSMOD/.  Program  CSPM  calls  the  appropriate 
secondary  overlays  to  accomplish  program  execution. 


Figure  3-1.  CSPM  overlay  structure 
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The  (1,0)  secondary  overlay,  with  program  INPUT  as  the  overlay  control- 
ler, reads  the  input  in  free  format,  then  lists,  decodes,  verifies,  and 
limit  checks  the  input  data.  An  input  data  summary,  preceded  by  error  mes- 
sages as  appropriate,  also  is  printed. 

Path  profile  generation  and  path  parameter  calculation  routines,  along 
with  the  basic  transmission  loss  model,  are  contained  in  the  (2,0)  overlay 
with  program  LOSSMOD  as  the  overlay  controller.  A topographic  data  file, 
containing  terrain  elevation  at  30"  intervals  in  latitude  and  longitude 
for  the  continental  United  States  (CONUS) , is  accessed  to  obtain  data  to 
generate  the  path  profile.  The  path  profiles  provide  the  necessary  informa- 
tion for  calculating  the  path  parameter  data  (elevation  of  and  distance 
to  the  radio  horizon  for  each  antenna  and  the  average  elevation  of  the  terrain 
foreground  for  each  antenna)  required  as  input  to  the  point-to-point  propa- 
gation loss  model.  The  loss  model  is  the  Longley-Rice  model  used  in  the 
point-to-point  mode.  The  loss  model  computes  statistical  estimates  of  basic 
transmission  loss  as  a function  of  distance.  The  final  computation  providing 
the  basic  transmission  loss  is  done  in  the  (3,0)  secondary  overlay,  where 
all  output  calculations  are  performed.  But  all  parameters  required  for 
that  calculation  are  computed  in  the  (2,0)  overlay  and  written  to  a file 
for  use  in  the  calculations  of  not  only  basic  transmission  loss  but  also 
power  density,  received  signal  level,  and  communication  reliability  in  the 
(3,0)  overlay.  Single  path  or  geographical  area  data  are  plotted  and/or 
tabulated  in  the  (3,0)  overlay,  according  to  specifications  by  the  user. 

3.2  Description  of  Input  Data 

The  input  data  serve  to  define  the  type  of  output  desired  by  the  model 
user  as  well  as  the  conditions  to  be  applied  in  the  calculations  of  basic 
transmission  loss,  power  density,  received  signal  level,  and  communication 
reliability.  No  attempt  is  made  in  this  section  to  discuss  the  physical 
basis  or  rationale  of  the  propagation  loss  model.  Rather,  some  familiarity 
with  the  model  is  assumed.  If  that  familiarity  is  lacking,  the  reader  is 
referred  to  Appendix  C and  the  references  by  Longley  and  Rice  (1968) , Rice 
et  al.  (1967),  Longley  (1976),  private  communications  with  Longley  et  al. 
(1971),  and  Hufford  (1972)  (available  from  the  authors)  which  discuss  the 
model  in  considerable  detail. 
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The  first  input  data  card  must  be  a title  for  the  run.  Following 
the  title  and  starting  with  the  next  card  of  the  input  data  deck,  the 
input  data  may  be  specified  in  any  order  with  multiple  items  of  data  on 
each  card.  However,  an  item  of  data  may  not  begin  on  one  card  and  end 
on  another  card.  The  data  are  specified  by  using  keywords  (shown  in  capital 
letters  in  this  discussion)  for  each  item  of  data.  Default  values  for 
most  of  the  data  items  are  built  into  the  model. 

Title.  The  title,  provided  by  the  first  input  data  card,  will  be 
printed  on  the  output  data  tabulations  and  plots.  Up  to  60  characters, 
including  blanks,  may  be  used. 

Type  of  output  data.  Various  selections  for  the  output  data  are 
available.  The  following  list  indicates  most  of  the  possible  selections. 

-Tcdsulated  path  profile  data. 

-Plotted  path  profile  data. 

-Tabulated  values  of  LB,  PD,  RSL,  and/or  REL  for  a geographical  area, 
with  or  without  path  profile  data  for  selected  radials. 

-Plotted  values  of  LB,  PD,  RSL,  and/or  REL  for  a geographical  are^, 
with  or  without  path  profile  data  for  selected  radials. 

A specification  of  PATH  will  provide  path  profile  data  for  a single 

propagation  path.  Basic  transmission  loss  also  may  be  obtained  for  the 

path  by  specifying  LB.  Specifying  MAP  will  result  in  the  production  of 

o o 

I£,  data  for  a geographical  area  1 in  latitude  by  1 in  longitude  with 

the  transmitter  located  in  the  center  of  the  area.  Power  density,  received 

signal  level,  or  communication  reliability  also  may  be  obtained  for  the 

geographical  area  by  appropriately  specifying  PD,  RSL,  or  REL.  Specifying 

MAP  = (1st  latitude/2nd  latitude/lst  longitude/2nd  longitude)  will  produce 

LB,  PD,  RSL,  and/or  REL  data  for  the  geographical  area  thus  defined. 

o o 

The  maximum  allowed  area  size  is  2 latitude  by  2 longitude.  The  transmitter 
may  be  specified  at  any  location  within  the  area. 

Two  additional  keywords  which  pertain  to  the  output  of  path  profile 
data  are  PRO,  denoting  plotted  profile  data,  and  PTAB,  denoting  tcdsular 
profile  data.  When  the  input  data  contain  PATH,  PRO,  only  plotted  profile 
data  are  provided.  Specifying  PATH  means  that  PTAB  will  result  in  only 
tcd>ulated  profile  data.  However,  specifying  only  PATH,  the  output  will 
be  both  plotted  and  tabulated  path  profile  data. 


r 


PLOT  and  TABL  are  keywords  applying  to  the  geographical  area  data  in  a 
way  functionally  similar  to  the  use  of  PRO  and  PTAB  for  path  profile  data. 

The  desired  type(s)  of  output  data  for  the  geographical  area  are  obtained 
by  appropriately  specifying  LB  (for  basic  transmission  loss) , PD  (for  power 
density) , RSL  (for  received  signal  level) , and/or  REL  (for  communication 
reliability) . 

One  also  may  obtain  path  profile  data  for  the  geographical*area  (that  *■ 
is,  MAP  has  been  specified)  by  specifying  PRO  = a number  and  PTAB  = a number. 
Only  certain  numbers  are  acceptable  as  shown  in  Table  3-1.  For  exan^le, 
if  profile  data  are  desired  for  radial  paths  extending  from  the  transmitter 
location  to  the  area  ]x>undary  at  uniformly  spaced  azimuth  angle  increments 
of  22.5°,  the  numijer  16  must  be  entered  with  the  keyword,  i.e.,  PRO=16  or 
PTAB=16. 


Table  3-1.  Valid  Numbers  for  TABL,  PRO,  and  PTAB  (Area  Option) 


Number  Following 

Bearing  Angle  Increment 

TABL,  PRO,  and  PTAB 

Between  Radials  Degrees 

0 

... 

2 

180 

4 

90 

8 

45 

16 

22.5 

32 

11.25 

64 

5.63 

128 

2.81 

256 

1.41 

512 

0.70 

1024 

0.35 

2048 

0.18 

Trcmsmitter  location.  Two  keywords  provide  information  on  transmitter 
location.  These  keywords  are  XMTRLOC  and  XMTRCOOR.  XMTRLOC  is  used  to 
define  word  identification  of  the  tremsmitter  location  on  the  output  data 
simply  as  a convenience  to  the  user.  The  information  otherwise  is  not 
used  by  the  model.  XMTRCOOR  is  used  to  specify  data  which  define  latitude 
euid  longitude  coordinates  for  the  treuismitter . Transmitter  location  coordinates 
are  required  data,  and  there  is  no  default  for  this  item.  When  a geographical 


area  is  defined,  the  XMTRCOOR  data  must  be  on  or  within  the  area  Isoundary. 
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Receiver  location.  The  keywords  defining  receiver  location  are  RCVRLCX: 
and  RCVRCCX)R.  RCVRLOC  provides  word  identification  of  the  receiver  location 
on  output  data  as  a user  convenience.  RCVRCOOR  is  required  data  which 
define  latitude  and  longitude  coordinates  for  the  receiver.  There  is  no 
default  for  this  item.  However,  receiver  location  data  are  required  and 
used  only  when  PATH  has  been  specified  to  obtain  path  profile  data  for 
a specific  propagation  path. 

Frequency . The  carrier  frequency  of  a radio  transmission  must  be 
specified  through  use  of  the  keyword  FREQ.  The  input  specification  for 
carrier  frequency  must  be  within  the  range  of  20  MHz  to  20,000  MHz.  The 
model  contains  a default  value  of  300  MHz. 

Transmitter  power.  A specification  of  transmitter  power  - actually 
power  delivered  to  the  input  port  of  the  transmitting  antenna  - is  provided 
by  the  keyword  PWR.  Transmitter  power  may  be  specified  using  either  watts 
or  decibels  referenced  to  a milliwatt  (dBm)  as  the  power  units.  The  model 
default  value  is  1 W (30  dBm) . 

Antenna  heights.  The  structural  antenna  heights  are  specified  using 
the  keywords  HGl  for  the  transmitting  antenna  and  HG2  for  the  receiving 
antenna.  These  heights  may  be  specified  using  units  of  either  meters  or 
feet.  The  default  value  for  each  antenna  is  10  m. 

Antenna  gains.  The  gain  of  each  antenna  in  the  direction  of  the  other 
antenna  (no  general  antenna  gain  model  is  used)  is  an  input  to  the  model. 
Keywords  for  the  gain  of  each  antenna  are  G1  for  the  transmitting  antenna 
and  G2  for  the  receiving  antenna.  These  gain  values  must  be  specified 
in  decibels  referenced  to  the  gain  of  an  isotropic  radiator  or  receptor 
(dBi)  of  electromagnetic  energy.  Default  values  are  10  dBi  for  the  transmit- 
ting antenna  and  0 dBi  for  the  receiving  antenna. 

Polarization.  Polarization  of  the  energy  radiated  by  the  transmitting 
antenna  is  specified  through  use  of  the  keyword  POL.  Either  horizontal 
or  vertical  polarization  may  be  selected.  The  default  condition  assumes 
horizontal  polarization. 

Ground  description.  The  model  provides  two  schemes  for  defining  the 
required  ground  constants.  One  scheme  allows  the  user  to  input  a ground 
descriptor  associated  with  the  keyword  GND.  The  choices  for  ground  descriptor 
are  shown  in  Tctble  3-2. 
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Table  3-2.  Ground  Descriptors  and  Electrical  Constants 


Descriptor 

GND 

Conductivity 

(S/m) 

Dielectric 

Constant 

Good  Ground 

GOOD 

0.02 

25 

Average  Ground 

AVG 

0.005 

15 

Poor  Ground 

POOR 

0.001 

4 

Sea  Water 

SEA 

5 

81 

Fresh  Water 

WATER 

0.01 

81 

Concrete 

CONCRETE 

0.01 

5 

The  second  scheme  requires  the  user  to  define  the  conductivity  using  the 
keyword  SGM  and  the  dielectric  constant  using  the  keyword  EPS.  The  default 
condition  assiMes  average  ground  with  SGM  = 0.005  S/m  and  EPS  = 15.  The 
allowable  ranges  are  0 to  100  S/m  for  conductivity  and  1 to  200  for  the 
dielectric  constant. 

Atmospheric  refractivity . Two  forms  of  refractivity  data  are  accepted 
by  the  model.  The  user  may  specify  surface  refractivity  (N^)  using  the 
keyword  NS,  or  he  may  specify  surface  refractivity  referenced  to  mean  sea 
level  (Nq,  sometimes  termed  reduced  surface  refractivity)  using  the  keyword 
NO.  When  refractivity  referenced  to  mean  sea  level  is  specified,  it  is 
converted  to  a surface  refractivity  value  as  described  in  Appendix  C and 
the  propagation  loss  model  references  cited  earlier  in  this  section. 

Climate.  An  appropriate  descriptor  for  climate  in  the  geographical 
region  containing  the  communication  system (s)  being  designed  or  evaluated 
is  specified  using  the  keyword  CLIM.  The  available  choices  are  shown  in 
Table  3-3. 


Table  3-3.  Climate  Descriptors 


Equatorial 

EQUATOR 

Continental  Sub-Tropical 

CNTROP 

Maritime  Sub-Tropical 

MRTROP 

Desert 

DESERT 

Continental  Temperate 

CNTEMP 

Maritime  Temperate,  Overland 

MRTMPLND 

Maritime  Temperate,  Oversea 

MRTMPSEA 

The  default  condition  for  climate  is  continental  temperate  (CNTEMP) . 
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Propagation  loss  variabilities.  There  are  three  contributors  to  the 
uncertainty  of  propagation  loss  estimates  calculated  by  the  point-to-point 
propagation  loss  model.  (Since  basic  treinsmission  loss  is  used  in  the  calcu- 
lation of  power  density  and  received  signal  level,  the  variabilities  also 
influence  these  estimated  quantities;  see  Appendix  F.)  The  first  influence 
considered  by  the  model  is  that  of  long-term  (greater  than  one  hour)  fading, 
and  the  terminology  usually  applied  to  this  variability  is  time  availability, 
denoted  by  q^.  Time  availability  expresses  the  fraction  of  time  during 
which  hourly  median  propagation  loss  will  not  exceed  a given  value  due  to 
long-term  fading  effects.  The  second  variability  considered  is  the  path- 
to-path  variation  in  loss.  The  term  applied  to  this  influence  is  location 
variability,  denoted  by  q . This  variability  expresses  the  fraction  of 

Li 

similar  paths  for  which  propagation  loss  will  not  exceed  a given  value  for 
a specified  fraction  of  the  time.  The  third  contribution  to  variability 
often  is  termed  the  confidence  (sometimes  prediction  error) , denoted  by  the 
symbol  Q.  Including  the  confidence  in  the  model  is  an  aclcnowledgment  of 
influences  which  are  not  included  explicitly  in  time  availaJaility  and 
location  variability. 

This  brief  bac)cground  to  variabilities  affecting  propagation  loss  allows 

one  to  realize  that  each  prediction  of  propagation  loss  is  an  estimate  of 

loss  which  will  not  be  exceeded  (for  at  least)  a fraction  q^  of  the  time 

at  a fraction  q of  the  locations,  estimated  with  confidence  Q.  These 
L 

Vctriabilities  are  not  independent,  hence  they  are  specified  in  sets.  The 
CSPM  input  data  can  include  up  to  three  such  sets  of  specified  variabilities. 
The  )ceywords  assigned  to  these  sets  of  variabilities  are  VI,  V2,  and  V3. 

The  model  default  condition  defines  one  set  of  variabilities  by  assuming 
median  values  for  each  parameter  of  the  set,  i.e.,  q^  = 0.50,  q^  = 0.50, 
and  Q = 0.50. 

Comnunication  reliability  output.  An  output  option  of  the  CSPM  is 
plotted  and/or  tabular  data  reflecting  estimated  communication  reliability. 
(See  Appendix  F for  definition  and  discussion  of  communication  reliability.) 
Four  levels  of  communication  reliability,  to  be  plotted  as  equal-value 
contours,  may  be  specified  as  part  of  the  input  data  by  using  the  )ceyword 
RQ.  Any  values  greater  than  zero  cind  less  than  one  may  be  selected.  The 
default  condition  provides  output  data  for  communication  reliability  values 
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or  0.10,  0.50,  0.90,  and  0.99. 

Threshold  for  received  signal  level.  The  concept  of  conununication 
reliability,  identified  in  the  preceding  paragraph,  makes  sense  only  in 
the  context  of  communication  system  performance  exceeding  some  threshold 
of  performance.  In  the  CSPM  the  index  or  measure  of  expected  system  perform- 
ance is  received  signal  level.  Therefore,  when  received  signal  level  exceeds 
some  threshold  for  received  signal  level  (for  specified  conditions  of  time 
availability  and  location  variability) , the  condition  is  satisfied  with 
some  probability  termed  communication  reliability  in  this  report.  A CSPM 
user  either  must  specify  the  received  signal  level  threshold  for  his  problem 
or  accept  the  default  value  of  -88  dBm.  The  keyword  for  specifying  received 
signal  level  threshold  is  RSLTH. 

3.3  Description  of  Output  Data 

Output  data  will  describe  a single  propagation  path  or  a geographical 
area  for  which  propagation  and  system  performance  data  are  desired.  The 
various  possible  types  of  output  data  were  listed  in  the  preceding  section: 
they,  also,  are  shown  in  Figure  3-1.  In  this  section  the  data  will  be 
described  first  for  a single  path,  since  such  data  are  required  for  develop- 
ing basic  transmission  loss,  power  density,  received  signal  level,  and  commu- 
nication reliability  data  for  a geographical  area. 

3.3.1  Single  Path  Data 

As  has  been  stated,  the  propagation  path  profile,  generated  by  retrieving 
terrain  elevation  data  from  the  digitized  topographic  data  file,  is  available 
in  plotted  and  tabular  forms.  These  data,  shewing  terrain  elevation  at 
points  along  the  path,  are  calculated  at  approximately  250  m intervals. 

At  each  point  along  the  path,  the  topographic  data  which  provide  elevation 
at  30  sec  intervals  in  latitude  and  longitude  are  bilinearly  interpolated  to 
determine  path  elevation.  Figure  3-2  is  a typical  plotted  path  profile, 
and  Table  3-4  shows  the  tabular  form  of  the  same  path  profile  data. 

One  will  observe  that  the  abscissa  of  Figure  3-2  is  not  linear.  Rather, 
the  plot  is  produced  with  a curvature  which  corresponds  with  the  radio  ray 
bending  that  would  occur  for  the  specified  atmospheric  refractivity . The 


17 


convenience  of  such  a plot  is  that  straight  lines  can  be  dravm  on  the  plot 
to  represent  the  radio  rays.  For  example,  the  radio  horizon  and  elevation 
angle  can  be  determined  graphically  by  drawing  a straight  line  from  an  antenna 
to  the  highest  visible  terrain  point  in  the  foreground  of  that  antenna. 

Another  important  feature  of  Figure  3-2  is  that  the  left-hand  scale  for 
elevation  is  in  meters,  while  the  right-hand  scale  is  in  feet. 

As  mentioned,  the  same  data  in  tabular  form  are  shown  in  Table  3-4. 

For  each  path  profile  point,  distance  along  the  geodesic  and  elevation  are 
presented  in  both  International  and  English  units.  Finally,  note  that  infor- 
mation about  the  path  is  printed  as  a heading  to  the  tabulated  data.  Basic 
transmission  loss  for  the  path  follows  the  tabulated  path  profile  data.  This 
is  an  optional  output  features. 

3.3.2  Geographical  Area  Data 

The  geographical  area  data,  like  the  path  profile  data,  are  produced 
in  plotted  emd  tabular  forms.  These  data  reflect  basic  transmission  loss, 
power  density,  received  signal  level,  and  communication  reliability  for 
the  geographical  area  containing  a defined  transmitter  (or  receiver)  location. 

These  data  are  developed  by  defining  a large  number  of  points  equally 
spaced  along  radial  paths  extending  from  the  defined  location.  The  ^mgular 
increment  for  radial  paths  has  been  defined  in  a way  to  achieve  approximately 
uniform  density  of  points  throughout  the  geographical  area.  The  reader 
is  referred  to  Appendix  F for  definition  and  discussion  of  the  algorithm 
used  in  selecting  the  points  at  which  calculations  are  made. 

A path  profile  (such  as  discussed  in  the  preceding  section)  is  developed 
for  each  radial  path.  These  profile  data  are  used  to  compute  the  necessary 
path  parameter  data  used  by  the  Longley-Rice  propagation  loss  model  in  the 
point-to-point  mode  to  conpute  basic  transmission  loss  to  each  point  thus 
defined.  The  basic  transmission  loss  data  then  are  used  to  calculate  power 
density,  received  signal  level,  and  comnunication  reliability  data,  discussed 
in  Appendix  F. 

Plotted  output  data  provide  equal-value  contours  of  basic  transmission 
loss,  power  density,  received  signal  level,  or  communication  reliidsility . 

The  contour  levels  for  basic  transmission  loss,  power  density,  and  received 
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Table  3-4.  Tabulated  Path  Profile  Data  (for  Path 
Profile  Shown  in  Figure  3-2) 
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Table  3-4.  Continued 


Table  3-4 . Continued 
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signal  level  are  selected  to  span  the  minimum-value-to-maximum-value  interval 
using  up  to  ten  subintervals.  Figure  3-3  is  a typical  plot  of  contour  data. 
Four  user-selected  values  of  communication  reliability  are  contoured  when 
that  option  is  selected. 

Tabulated  output  data  provide  the  calculated  values  of  basic  transmis- 
sion loss,  power  density,  received  signal  level,  and  communication  reliability 
at  points  spaced  1 km  apart  along  selected,  uniformly  spaced  radials.  Table 
3-5  is  a typical  example  of  tabulated  data.  This  form  of  output  data  can 
become  voluminous,  consequently  the  user  is  cautioned  to  not  specify  need- 
lessly large  numbers  of  radials  for  tabular  data. 

3.4  Limitations  of  the  Model 

The  CSPM  has  been  designed  with  the  fewest  possible  constraints,  so 
that  a broad,  general  utility  can  be  realized.  Several  of  the  model  limi- 
tations are  imposed  by  the  propagation  loss  model.  These  limitations 
include: 

— Carrier  frequencies  between  20  and  20,000  MHz. 

— Antenna  heights  between  0.5  and  3000  m. 

— Vertical  or  horizontal  polarization  for  radiated  energy. 

In  addition,  the  following  constraints  have  been  programmed  into  the  model; 

— A maximum  of  three  sets  of  variabilities  (for  time  availability, 
location  variability,  and  confidence)  for  each  analysis  run. 

— A maximum  of  four  levels  of  communication  reliability  plotted  per 
run . 

— A maximum  geographical  area  of  2°  latitude  by  2°  longitude  for 
contoured  output  data. 

Finally,  the  topographic  data  impose  a restriction  on  the  model's 
capability  to  consider  the  fine-grain  character  of  the  terrain,  since  granu- 
larity of  the  digitized  topographic  data  is  30"  in  latitude  and  longitude 
and  precision  of  the  associated  elevation  values  is  20  ft  (6.1  m). 
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4.  SUMMARY  OF  MODEL  FEATURES  AND  OUTPUT  OPTIONS 

1.  A con^uter-based  method  has  been  developed  which  uses  a point-to-point 
propagation  loss  model  to  evaluate  communication  system  performance 
within  a defined  geographical  area.  In  addition,  propagation  path 
characteristics  for  specific  paths  can  be  obtained  from  the  model. 

2.  The  model  is  user-oriented,  which  means  that  input  information  may 

be  supplied  in  free  format,  and  various  options  are  allowed  for  describing 
and  specifying  the  input  characteristics  and  data. 

3.  The  model  provides  a capability  for  plotting  equal-value  contours 
over  a geographical  area.  The  data  which  may  be  plotted  include  basic 
transmission  loss,  power  density,  received  signal  level,  and  communication 
reliability.  Tabulated  output  data  also  may  be  specified  for  each 
output . 

4.  The  model  includes  a file  of  digitized  topographic  characteristics 

for  the  contiguous  United  States.  Terrain  elevation  data  are  retrieved 
from  this  file  and  used  in  developing  required  input  data  for  the 
point-to-point  propagation  loss  model. 

5.  The  model  is  structured  so  as  to  facilitate  the  use  of  additional 
and/or  improved  data  as  they  become  available.  This  feature  is  par- 
ticularly important  when  considering  the  use  of  an  improved  propagation 
loss  model,  new  or  improved  terrain  data  which  would  provide  greater 
detail  in  describing  the  topography  of  an  area,  new  statistical  charac- 
terizations within  the  model,  and  new  indices  (or  measures)  of  satis- 
factory system  performance. 

6.  The  model  is  capcdsle  of  very  general  application  and  utility.  Limi- 
tations of  the  model  are  those  imposed  by  the  propagation  loss  model, 
the  granularity  of  the  topographic  data,  and  practical  restrictions 
arbitrarily  progr£unmed  into  the  model  to  maintain  reasonable  limits 
on  the  volume  of  output  data. 

7.  The  propagation  loss  model  calculates  basic  trctnsmission  loss  as  a 
random  function  of  time  and  space  by  modeling  atmospheric,  climate, 
and  terrain  irregularity  influences  upon  propagation.  The  model  is 
good  for  frequencies  between  20  and  20,000  MHz. 
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5 . RECOMMENDATIONS 


Although  the  capability  provided  by  the  Comnunication  System  Performance 
Model  is  a valuable  addition  to  the  automated  design  evaluation  aids  available 
to  a communication  system  engineer,  to  further  improve  and  extend  the  useful- 
ness of  the  CSPM,  it  is  recommended  that; 

1.  The  feasibility  be  studied  for  adapting  the  CSPM  for  interactive 
operation  from  a remote  terminal. 

2.  The  influence  of  rainfall  upon  radio  signal  attenuation  be  incor- 
porated into  the  statistics  of  the  point-to-point  propagation  loss 
model  used  by  the  CSPM. 

3.  An  improved  terrain  characteristics  data  base  be  added  to  the 
model,  so  that  considerably  greater  detail  of  the  terrain  character- 
istics and  man-made  obstructions  can  be  considered  in  generating 
the  propagation  path  profiles  used  by  the  point-to-point  propaga- 
tion loss  model. 

4.  A revised  point-to-point  propagation  loss  model  be  used  to  replace 
the  Longley-Rice  Model,  used  in  the  point-to-point  mode  in  the 
present  CSPM. 

5.  A map  of  surface  refractivity  referenced  to  mean  sea  level  (NO) 
be  developed  for  the  CSPM. 

Recommendation  three  may  involve  handling  much  larger  files  of  terrain 
information  than  the  present  data  base  requires.  It  may  be  necessary,  there- 
fore, to  consider  use  of  some  mass  storage  medium  other  than  magnetic  tape — 
for  example,  permanent  disc  pack.  Similarly,  recommendation  one  could  involve 
requirements  for  a callable  storage  medium  other  than  magnetic  tape,  if 
magnetic  tape  calls  from  a remote  terminal  are  not  allowed. 
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APPENDIX  A.  USER’S  GUIDE 
COMMUNICATION  SYSTEM  PERFORMANCE  MODEL 


1 

I 


This  appendix  describes  how  to  use  the  Communication  System  Perform- 
ance Model  (CSPN)  on  a CDC  6000  series  computer  system.  The  model  has  been 
designed  with  free-format  input  and  mnemonic  keywords  to  aid  the  user  in 
describing  the  system  to  be  analyzed  in  familiar  engineering  terms.  Options 
are  available  to  select  the  type  of  output  desired,  as  well  as  to  describe 
the  characteristics  of  the  geographical  area  and  communication  system  being 
studied.  A run  title  card  is  required  the  first  of  the  data  cards. 

Following  the  title  card,  the  CSPM  input  data  may  be  given  in  any  order 
on  the  subsequent  data  cards.  If  more  than  one  item  of  input  data  is  specified 
per  card,  commas  must  be  used  as  separators.  Blanks  always  are  ignored 
except  on  the  title  card.  Error  messages  or  warning  diagnostics  are  issued 
if  conflicting  options  are  selected,  if  em  option  is  specified  more  than 
once,  or  if  values  selected  for  a parameter  are  out  of  range  for  the  model. 
Default  values  have  been  included  in  the  model  for  all  parameters  except 
the  transmitter  and  receiver  coordinates. 

The  CSPM  may  be  used  in  two  different  modes  - for  a single  path  or 
for  an  area  defined  around  the  transmitter.  For  either  mode,  the  user  may 
select  output  options  that  produce  plotted  path  profiles,  tabulated  path 
profile  data,  and  a t^d>ulated  summary  of  basic  transmission  loss  for  the 


path  or  area.  When  the  analysis  is  done  for  fui  area,  contour  plots  of  basic 
transmission  loss,  power  density,  received  signal  level,  and  communication 
reliability  over  the  desired  area  are  available.  It  also  is  possible  to 
have  tabular  printed  summaries  of  each  type  o'f  output  data  throughout  the 
area.  The  single  path  mode  is  selected  by  the  keyword  PATH;  the  geographical 
area  mode  is  selected  by  MAP.  No  specification  of  mode  results  in  the  default 
condition  which  is  MAP.  A list  of  the  options  that  a user  may  select,  appro- 
priate to  each  mode,  is  given  below.  Examples  of  CSPM  input  and  output 
follow  later  in  this  appendix. 

A.l  Single  Path  Mode 

As  discussed,  specifying  the  keyword  PATH  selects  the  single  path  mode 
of  the  CSPM.  Input  data  pertinent  to  this  mode  are: 


mCSDINO  PAOl  BUMC-NDT  fl 


XMTRLOC  and  RCVRLOC 

These  optional  keywords  allow  the  user  to  specify  an  alphanumeric  name 
of  up  to  10  characters  to  be  associated  with  the  transmitter  and  receiver 
locations  on  the  plotted  and  tabulated  output. 

PRO  and  PTAB 

The  keyword  PRO  indicates  plotted  path  profile  data  are  desired.  The 
keyword  PTAB  indicates  tabulated  path  profile  data  are  desired.  In 
the  single  path  mode,  the  default  condition  provides  both  forms  of 
output.  To  select  one,  but  not  the  other,  the  user  should  include 
the  desired  keyword  in  his  input  data.  This  selection  will  cause  the 
other  option  to  be  defaulted  off.  An  alternative  specification  for 
either  option  is  PRO=0  or  PTAB=0  which  means  that  the  indicated  form 
of  output  will  not  be  provided. 

LB 

The  user  may  request  that  basic  transmission  loss  for  the  path  be  printed 
by  including  this  keyword.  In  the  single  path  mode,  the  default  condition 
is  no  basic  transmission  loss  calculation.  If  LB  is  selected,  then 
the  user  also  may  enter  any  values  tliat  he  desires  for  FREQ,  POL, 

HGl,  HG2,  CLIM,  GND,  or  EPS/SGM,  NS  or  NO,  and  the  variabilities  VI, 

V2,  or  V3.  Otherwise,  the  CSPM  default  values  will  be  used.  Please 
refer  to  the  next  section  of  this  appendix  for  discussion  of  these 
parameters . 


A. 2 Geographical  Area  Mode 

Specifying  the  keyword  MAP  selects  the  mode  of  the  CSPM  that  analyzes 
communication  system  performance  over  an  area  around  the  transmitter  coordi- 
nates. If  only  the  keyword  MAP  is  specified  or  if  the  default  option  is 
desired,  the  area  to  be  analyzed  is  defined  as  a 1°  by  1°  "box"  with  boundaries 
formed  by  the  latitude  and  longitude  values  0.5°  In  each  direction  from 
the  transmitter  coordinates.  When  a particular  area  Is  desired  (the  transmit- 
ter may  be  located  anywhere  within  the  area) , the  user  must  specify  MAP 
followed  by  two  latitude  values  and  two  longitude  values  which  define  the 
area  boundaries,  i.e.,  MAP  « (31, 27, 00,N/31, 34, 00,N/110, 26, 30,W/  110, 17, 30, W) . 
The  two  latitudes  are  given  first,  and  the  values  are  in  positive  degrees, 
minutes,  and  seconds  with  the  hemisphere  given  as  N or  S and  the  longitude 


with  respect  to  the  Prime  Meridian  given  as  E or  W.  In  all  cases,  the  area 
must  lie  within  the  bounds  of  the  data  available  in  the  topographic  data 
base,  which  is  an  area  from  23°  to  51°  N latitude  and  from  60°  to  130°  W 
longitude.  Other  input  data  generally  applicable  to  the  area  mode  of  model 
operation  are: 

PLOT  and  TABL 

The  keyword  PLOT  indicates  that  plotted  (contoured)  output  is  desired. 

The  keyword  TABL  indicates  that  tabulated  output  describing  basic  trans- 
mission loss,  power  density,  received  signal  level,  or  communication 
reliability  is  described.  The  default  condition  is  PLOT.  A contour 
plot  will  be  produced  for  each  type  of  data  selected  (LB,  PD,  RSL, 
and  REL)  and  for  each  set  of  variabilities  (VI,  V2,  and  V3)  that  was 
specified  in  the  input  data  for  this  run.  If  no  plots  are  desired, 
specify  PLOT=0.  When  TABL  is  specified,  the  tabulated  data  will  be 
provided  only  for  selected  radials  extending  at  equal-value  bearing 
angle  increments  from  the  transmitting  antenna  location  throughout 
the  area.  A separate  column  is  produced  for  each  type  of  data  for 
each  selected  variability  set  (Vl,  V2,  and  V3) . The  radial  paths  are 
selected  by  a number  following  the  keyword,  i.e.,  TABL=32.  The  only 
acceptable  numbers  are  powers  of  two.  Table  A-1  shows  the  acceptable 
numbers  and  the  bearing  angle  increment  between  radials  that  are  realized 
for  each  number.  When  a number  is  used  which  is  not  a power  of  two, 
the  next  lower  value  corresponding  to  a power  of  two  will  be  used. 

The  default  condition  for  TABL  is  16.  The  model  user  is  cautioned 
to  consider  carefully  the  value  used  with^ TABL,  for  the  output  could 
become  voluminous  for  large  numbers. 

PRO  and  PTAB 

These  keywords  are  valid  for  the  area  mode  of  operation  as  well  as 
for  the  single  path  mode,  and  the  same  meaning  accompanies  each  keyword. 
However,  the  plotted  or  tabulated  path  profile  data  will  be  provided 
only  for  selected  radial  paths  designated  by  a number  following  the 
keyword  as  discussed  for  TABL.  The  same  numbers  are  valid  as  for  TABL, 
but  the  default  value  for  PRO  and  PTAB  is  eight.  In  the  area  mode 
of  operation,  no  path  profile  data  will  be  output  unless  the 
specification  PRO”n  or  PTAB=n  is  given,  where  n denotes  a number. 
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The  same  rule  discussed  for  TABL  for  values  other  than  powers  of  two 
is  followed  by  PRO  emd  PTAB. 

Table  A-1.  Valid  Numbers  for  TABL,  PRO,  and  PTAB  (Area  Option) 


NumJaer  Following 

TABL,  PRO,  and  PTAB 

Bearing  Angle  Increment 
Between  Radials  Degrees 

0 

2 

180 

4 

90 

8 

45 

16 

22.5 

32 

11.25 

64 

5.63 

128 

2.81 

256 

1.41 

512 

0.70 

1024 

0.35 

2048 

0.18 

LB,  PO , RSL,  and  R£jL 

These  keywords  define  the  CSPM  outputs  of  basic  transmission  loss 

(LB) , power  density  (PD) , received  signal  level  (RSL) , and  communication 

reliability  (REL) , throughout  the  geographical  area  which  may  be  speci- 

i 

fied  by  the  model  user.  Since  basic  transmission  loss  is  required 
for  calculation  of  each  of  the  other  type  of  output,  LB  is  the  default 
type  of  output  data.  Any  combination  including  all  four  types  of 

( 

output  may  be  specified.  Basic  transmission  loss  (LB)  always  is  tabulated  ! 
when  TABL  is  specified.  i 

Discussion  of  the  remaining  input  data  pertinent  to  the  area  mode 
is  organized  to  present  first  the  input  required  to  calculate  Isasic  transmis-  1 

sion  loss.  The  additional  input  used  to  calculate  power  density  is  presented  ; 

next,  followed  by  the  additional  data  required  to  calculate  received  signal  I 

level,  and  finally  the  additional  input  required  for  calculating  communica-  ^ 

tion  reliability  is  discussed.  [ 
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A. 2.1  Input  Required  for  Basic  Transmission  Loss  (LB) 

XMTRCOOR 

Coordinates  of  the  transmitter  (actually  the  transmitting  antenna) 
location  are  required  input;  there  is  no  default.  XMTRCOOR  is  used 
to  specify  the  latitude  and  longitude  of  the  transmitter  in  positive 
degrees,  minutes,  and  seconds.  The  latitude  must  be  given  first  followed 
by  N or  S,  and  the  latitude  must  be  followed  by  E or  W,  i.e.,  XMTRCOOR= 
(31, 31, 51,N/110, 24, 10,W) . The  transmitter  coordinates  must  be  within  or 
on  the  boundaries  defined  by  deyword  MAP. 

XMTRLOC 

This  optional  keyword  may  be  used  to  assign  an  alphanumeric  name  of 
up  to  10  characters  to  the  transmitter  location,  i.e.,  XMTRLOC=HUACHUCA . 
This  provision  is  simply  a user  convenience  for  quick  identification 
of  plotted  and  printed  output. 

FREQ 

The  transmitter  radio  frequency  in  megahertz  is  specified  using  this 
keyword,  i.e.,  FREQ=237.865  MHz.  Limits  for  the  specification  of 
frequency  are  20  to  20,000  MHz.  The  default  for  frequency  is  300 
MHz.  Any  value  within  the  above  limits  and  up  to  10  characters  (includ- 
ing a decimal  point)  in  length  is  accepted  by  the  model.  Printed 
summaries  of  the  input  will  show  values  to  four  decimal  places. 

POL 

Polarization  of  the  radiated  energy  may  be  specified  using  the  keyword 
POL.  The  options  are  horizontal,  P0L=H,  and  vertical,  POL=V,  polari- 
zation. The  default  condition  is  horizontal  polarization. 

HGl  and  HG2 

These  keywords  are  used  to  specify  antenna  heights  2d>ove  ground  for 
the  transmitter  (HGl)  and  receiver  (HG2) . Heights  may  be  specified 
in  either  meters  or  feet,  i.e.,  HG1*‘15.3  M or  HG2=12.6  FT.  When  units 
(M  or  FT)  are  not  included  in  the  specification,  meters  are  assumed. 

Valid  limits  for  either  antenna  height  are  0.5  to  3000  m (1.64  to 
9843  ft) . The  default  for  each  height  is  10  m.  Values  within  the 
above  limits  up  to  10  characters  (including  the  decimal  point)  in 
length  are  accepted  by  the  model.  Printed  summaries  of  the  input 
will  show  values  to  one  decimal  place  for  meters  and  to  units  for  feet. 


Selection  of  a climate  code  describing  the  geographical  area  of  interest 
is  allowed  through  use  of  this  keyword,  i.e.,  CLIM=CNTEMP.  The  choices 
are  shown  in  Table  A-2  where  a descriptive  word  or  phrase  is  keyed 
to  a single  keyword  descriptor  to  be  used  with  CLIM.  The  default 
condition  is  continental  ten^jerate  climate. 

Table  A-2.  Climate  Descriptors 


Climate  Keyword 

Climate  Descriptor 

EQUATOR 

Equatorial 

CNTROP 

Continental  Sub-Tropical 

MRTROP 

Maritime  Sub-Tropical 

DESERT 

Desert 

CNTEMP 

Continental  Teii:5>erate 

MRTMLND 

Maritime  Temperate,  Overland 

MRTMSEA 

Maritime  Temperate,  Oversea 

EPS  and  SGM  or  GND 

The  ground  dielectric  constant  and  conductivity  may  be  specified  by 
using  keywords  EPS  and  SGM,  respectively,  i.e.,  EPS=11.9  and  SGM=0.0385. 
Valid  limits  for  the  dielectric  constant  (EPS)  are  1 to  200.  Valid 
limits  for  conductivity  (SGM)  are  0 to  100  S/m.  Values  within 
the  above  limits  up  to  10  characters  (including  the  decimal  point) 
in  length  for  each  variable  are  accepted  by  the  model.  Printed  summaries 
of  the  input  will  show  the  dielectric  constant  to  one  decimal  place 
and  the  conductivity  to  four  decimal  places.  An  alternative 
to  specifying  EPS  and  SGM  is  the  specification  of  a ground  descriptor 
using  the  keyword  GND,  i.e.,  GND=POOR.  Teible  A-3  shows  the  acceptable 
ground  descriptors  and  keywords  with  associated  values  for  dielectric 
constant  and  conductivity.  The  default  condition  is  average  ground. 


42 


Table  A-3.  Ground  Descriptors  and  Electrical  Constants 


Ground  Keyword 

Ground  Descriptor 

SGM  (S/m) 

GCX3D 

Good  Ground 

0.02 

AVG 

Average  Ground 

0.005 

POOR 

Poor  Ground 

0.001 

SEA 

Sea  Water 

5 

WATER 

Fresh  Water 

0.01 

CONCRETE 

Concrete 

0.01 

NS  or  NO 

Surface  refractivity  denoted  by  the  keyword  NS  or  surface  refractivity 
referenced  to  mean  sea  level  (reduced  surface  refractivity)  denoted  by 
NO  may  be  specified,  i.e.,  NS=327  or  N0=289.  Valid  limits  for  each 
variable  are  250  to  400  N-units.  The  default  value  is  301  N-units 
which  corresponds  to  the  common  assun^tion  of  effective  earth  radius 
equal  to  4/3  actual  earth  radius.  Printed  summaries  will  show  the 
value  to  one  decimal  place.  When  surface  refractivity  (NS)  values  are 
not  known  for  the  area  of  interest,  the  model  user  may  wish  to  use 
reduced  surface  refractivity  (NO)  values  from  world  maps  of  that  parame- 
ter. Such  maps  are  available  from  the  International  Radio  Consultative 
Committee  (CCIR,  1974)  or  in  the  NBS  Monograph  by  Bean  et  al  (1960) . 

VI,  V2,  and  V3 

The  time  variability,  location  variability,  and  prediction  confidence 
to  be  used  in  calculating  basic  transmission  loss  are  specified  at 
a variability  set,  i.e.,  Vl=( .99, .95, .90) . Time  variability  alvays 
relates  to  the  first  value  in  the  set,  location  variability  re  .a ted 
to  the  second  value,  and  prediction  confidence  relates  to  the  third 
value.  Up  to  three  such  sets  may  be  specified  per  run.  Values  must 
be  between  .001  and  .999.  If  a 0 is  specified,  it  will  be  changed 
to  .001,  and  1 will  be  changed  to  .999.  The  default  condition  is 
one  variability  set  using  median  values.  If  a value  is  omitted  from 
a set,  the  corresponding  value  from  the  last  set  (or  the  default  in 
the  case  of  VI)  will  be  used.  For  example:  VI  « (.2,. 3),  V2  = (,,.9) 
would  result  in  VI  = (.2,  .3,. 5),  and  V2  > (.2,. 3, .9)  being  read. 
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A. 2. 2.  Additional  Input  Required  for  Power  Density  (PD) 

All  input  discussed  in  the  preceding  parts  of  Section  A. 2 is  used  in 
the  calculation  of  power  density.*  In  addition,  the  following  input  items 
also  are  used. 

PWR 

The  specification  of  transmitter  output  power  in  watts  (W)  or  decibels 
referenced  to  a milliwatt  (dBm)  is  provided  using  the  keyword  PWR, 
i.e.,  PWR=100  W or  PWR=30.  When  no  units  (W  or  dBm)  are  specified, 
decibels  referenced  to  a milliwatt  are  assumed.  Limits  for  specif ica- 

9 

tion  of  transmitter  power  are  0 to  120  dBm  or  0.001  to  1x10  W.  The 
default  value  is  30  dBm.  Any  value  within  the  above  limits  up  to  10 
characters  (including  the  decimal  point)  in  length  is  accepted  by  the 
model;  however,  printed  summaries  of  the  input  will  show  values  to 
one  decimal  place. 

G1 

Transmitting  antenna  gain  in  the  direction  of  the  receiving  antenna 
may  be  specified  in  decibels  referenced  to  an  isotropic  antenna  (dBi) 
using  this  keyword,  i.e.,  Gl=23.  Limits  for  specifying  the  transmitting 
antenna  gain  are  0 to  100  dBi . Any  value  within  these  limits  up  to 
10  characters  (including  the  decimal  point)  in  length  is  accepted  by 
the  model.  Printed  summaries  of  the  input  will  show  the  value  to  one 
decimal  place.  The  default  value  for  the  transmitting  antenna  gain 
is  10  dBi. 

A. 2. 3.  Additional  Input  Required  for  Received 
Signal  Level  (RSL) 

All  input  discussed  thus  far  in  Section  A. 2 is  used  in  the  calculation 
of  received  signal  level.  In  addition,  the  following  input  item  also  is 
used. 

•As  shown  in  Appendix  F (F-3) , power  density  is  independent  of  frequency. 
However,  the  model  calculations  of  power  density  are  made  using  basic  trans- 
mission loss  data  computed  earlier  in  the  program  execution.  Frequency 
dependence  is  removed  as  shown  in  (F-4) . 
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G2 

Receiving  antenna  gain  in  the  direction  of  the  transmitting  antenna 
may  be  specified  in  decibels  referenced  to  an  isotropic  antenna  (dBi) 
using  this  keyword,  i.e.,  G2=6.  The  same  limits  and  discussion  of 
values  presented  in  the  preceding  section  for  the  transmitting  antenna 
gain  apply.  The  default  value  for  the  receiving  antenna  gain  is  0 dBi. 

A. 2. 4.  Additional  Input  Required  for  Communication 
Reliability  (REL) 

All  input  discussed  in  section  A. 2 is  used  in  the  calculation  of  com- 
munication reliability.  The  following  input  items  are  added  to  those  items 
already  discussed. 

RSLTH 

The  threshold  for  received  signal  level — a satisfactory  system  perform- 
ance crit  rion — may  be  specified  in  decibels  referenced  to  a milliwatt 
using  this  keyword,  i.e.,  RSLTH=-83.  Limits  for  specifying  the 
received  signal  level  threshold  are  -200  and  +50  dBm.  Any  value  within 
these  limits  up  to  10  characters  (including  the  decimal  point)  in  length 
is  acceptable  to  the  model,  but  printed  summaries  of  the  input  will 
show  the  value  to  one  decimal  place.  The  default  for  received  signal 
level  threshold  is  -88  dBm. 

RQ 

Up  to  four  contour  levels  for  plotting  communication  reliability  may 
be  specified  using  this  keyword,  i.e.,  RQ= ( .99, .95, .90, .50) . Values 
must  be  between  0.001  and  0.999;  however,  zero  will  be  changed  to  0.001 
and  one  will  be  changed  to  0.999.  The  default  values  for  communication 
reliability  are  .99, .90, .50,  and  .10.  If  a value  is  omitted  by  leaving 
its  space  vacant,  the  default  value  is  used.  For  example:  RQ  = (.2,,. 6) 
would  be  read  as  RQ  = (.2, .95,. 6). 

A. 3 Sample  Input  and  Output 

The  input  and  output  for  three  sample  applications  are  presented  in 
this  section  as  further  aid  to  users  of  the  CSPM.  The  first  sample  appli- 
cation is  use  of  the  model  to  generate  a single  path  profile  and  to  calculate 
basic  transmission  loss  for  the  path.  The  second  sample  illustrates  use 
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of  the  CSPM  to  produce  plotted  and  tabulated  output  of  received  signal  leve’ 
throughout  a geographical  area.  The  third  sample  application  illustrates 
plotted  communication  reliability  for  the  same  area,  but  some  different 
input  values  are  used. 

Input  data  (provided  by  the  user)  for  the  single  path  are  as  follows, 
where  each  line  contains  one  item  of  input:  PATH  PROFILE  TEST  (Run  title 
which  may  contain  up  to  60  characters) . 

XMTRCOOR  = (31,31,51,N/110,24,10,W) 

XMTRLOC  = HUACHUCA 

RCVRCOOR  = (32,12,38,N/110,58,31,W) 

RCVRLOC  = TUCSON 

PRO 

LB 

FREQ  = 376.125 
HGl  = 30  FT 
HG2  = 75  FT 
CLIM  = CNTEMP 
NO  = 290 
VI  = (.9, .9, .9) 

Note  that  POL  and  GND  (or  EPS  and  SGM)  have  not  been  specified;  hence,  the 
model  will  use  default  conditions  for  these  items.  The  default  conditions 
are  POL=H  and  GND=AVG  (EPS=0.005  and  SGM=15) . The  printed  input  summary 
is  shown  in  Table  A-4,  and  Figure  A-1  is  the  output  path  profile. 

Input  data  (provided  by  the  user)  to  obtain  plotted  and  tabulated  output 
indicating  received  signal  level  throughout  a geographical  area  are  as 
follows,  where  each  line  contains  one  item  of  input: 

RECEIVED  SIGNAL  LEVEL  TEST  (Run  title — up  to  60  characters)  j 

MAP  = (31,27.00,N/31,34,00,N/110,26,30,W/110,17,30,W)  j 

PLOT  I 

TABL  = 8 I 

RSL 

XMTRCOOR  = (31,31,51,N/110,24,10,W) 

XMTRLOC  = HUACHUCA 
FREQ  = 394.675 
POL  = V 
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Table  A-4.  Printed  Sunanary  of  Output  Data  for  Single  Path  Profile 


Table  A-4.  Continued 


Table  A-4 . Continued 


Table  A- 4 . Continued 


PWR  = 100  W 
HGl  = 20 
HG2  = 6 FT 
G1  = 15 
GND  = PCXDR 

VI  = (.90, .90, .90) , V2  = (.90, .90, .50) , V3  = ( . 90, . 50, . 50) 

Items  CLIM,  NS  (or  NO) , and  G2  also  are  used  in  calculating  received 
signal  level;  therefore,  the  default  values  are  used.  The  default  values 
are  CLIM  = CNTEMP,  NS  = 301,  and  G2  = 0.  Though  HGl  is  specified,  no  units 
are  included;  the  default  units  for  HGl  are  M.  The  printed  input  summary 
is  shown  in  Table  A-5.  Figures  A-2,  A-2,  and  A-4  are  the  plotted  received 
signal  level  outputs.  The  tabulated  received  signal  level  outputs  are 
included  as  Tables  A-6,  A-7,  A-8,  A-9,  A-10,  A-11,  A-12,  and  A-13. 

Input  data  (provided  by  the  user) , to  obtain  plotted  output  for  communi- 
cation reliability  throughout  the  same  geographical  area  as  used  for  the 
second  sample  application,  now  are  shown  with  each  line  containing  one 
item  of  input. 

COMMUNICATION  RELIABILITY  TEST  (Run  title—up  to  60  characters) 

MAP  = (31,27,00,N/31,34,00,N/110,26,30,W/110,17,30,W) 

REL 

XMTRCOOR  = (31,31,51,N/110,24,10,W) 

XMTRLOC  = HUACHUCA 
FREQ  = 245.798 
PWR  = 100  W 
HGl  = 10  M 
HG2  = 6 FT 
CLIM  = CNTEMP 
GND  = POOR 

VI  = (.95,. 95),  V2  = (.95,. 50),  V3  = (.50, .50) 

RSLTH  = -83.5 

RQ  = (.99,. 90,. 50,. 10) 

Note  that  output  data  form  (PLOT  or  TABL)  is  not  specified;  this  means 
that  the  default  provision  will  be  invoiced  which  provides  only  plotted 
output  for  the  geographical  area  mode  of  operation.  Also,  it  will  be 
observed  that  POL,  Gl,  and  G2  are  not  specified.  The  default  values  that 
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Figure  A-3,  Received  signal  level  contours  for  V2  • (.90,  .90,  .50). 

Typical  application  of  the  CSPM  in  the  geographical  area  mode. 
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Figure  A-4.  Received  signal  level  contours  for  V3  = (.90,  .50,  .50). 

Typical  application  of  the  CSPM  in  the  geographical  area  mode 
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Table  A-12.  Tabulated  Received  Signal  Levels  for  Radial  at  270 


Table  A-14.  Printed  Summary  of  Input  Data  for  Communication  Relieibility 
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are  used  are  POL  = H,  G1  = 10,  and  G2  = 0.  A printed  summary  of  the  input 
used  to  produce  the  communication  reliability  output  is  shown  in  Table 
A-4.  Figures  A-5,  A-6,  and  A-7  are  the  plotted  output  of  communication 
reliability. 
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APPENDIX  B.  COMMUNICATION  SYSTEM  PERFORMANCE  MODEL 
SOFTWARE  TECHNICAL  DESCRIPTION 


Technical  description  of  the  Communication  System  Performance  Model 
(CSPM)  software  given  in  this  section  is  intended  to  serve  as  a guide 
for  program  maintenance  and  modification.  The  first  section  includes 
a general  overview  of  progreim  structure,  processing  flow,  file  and  record 
descriptions,  and  common  block  definitions.  Documentation  for  each  individual 
routine  is  given  in  the  second  section.  The  third  section  describes  how 
to  install  the  CSPM  on  a CDC  6500  computer  system. 


B.l  General  Overview 

The  CSPM  has  been  developed  for  the  Army's  CEEIA  as  em  aid  in  designing 
and  evaluating  Army  communication  systems.  The  Longley-Rice  propagation 
prediction  loss  model  is  used  in  the  point-to-point  mode  to  determiiie 
basic  transmission  loss  statistics  which  are  used  to  calculate  power  density, 
received  signal  level,  and  communication  reliability.  These  values  may 
be  tckbulated  for  selected  paths  within  a geographic  area  or  used  to  produce 
contour  plots  over  the  area.  Tabulated  or  plotted  path  profiles  also  may 
be  requested  for  selected  paths  within  the  area.  If  a single  (discrete) 
path,  rather  than  a geographic  area,  is  to  be  analyzed,  the  CSPM  will  produce 
a tabulated  or  plotted  profile  and  calculate  the  basic  tmasmission  loss 
for  the  path.  A restart  option  is  available  to  allow  the  user  to  input 
new  statistical  values  and  output  options  to  be  used  with  the  data  from 
a LOSS  file  created  in  a previous  run.  The  CSPM  then  calculates,  tabulates, 
and  plots  another  set  of  values  for  basic  trcuismission  loss,  power  density, 
received  signal  level,  and  communication  reliability  for  the  same  geographic 
area. 

The  process  diagram  for  the  CSPM  (Figure  B-1)  shows  the  organization 
and  flow  of  the  four  major  parts  of  the  program.  The  model  is  divided 
into  a primary  overlay  and  three  secondary  overlays  to  minimize  the  core 
required  for  execution.  The  primary  overlay  (CSPM)  defines  the  files  and 
common  blocks  used  throughout  the  program.  It  controls  execution  by  calling 
the  secondary  overlays  as  they  are  needed.  The  (1,0)  overlay  (INPUT)  reads 
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Figure  B-1.  CSPM  program  orgctnization. 


the  free-format  input  data  cards,  checks  for  a restart  file,  determines 
if  any  errors  or  conflicts  exist  in  the  input  options,  stores  the  input 
values  in  the  /OPTIONS/  common  block,  and  prints  the  input  data  summary. 

The  (2,0)  overlay  (LOSSMOD)  extracts  path  profiles  from  the  topographic 
data  base  and  uses  the  Longley-Rice  propagation  model  to  calculate  ACR 
values  (attenuations  relative  to  free-space  losses)  and  statistical  values 
for  data  points  on  that  path.  The  (3,0)  overlay  (OUTPUT)  reads  the  profile 
and  ACR  data  from  file  LOSS  or  RESTRT,  if  it  is  a restart  run,  and  prints 
tabulated  values  or  produces  CALCOMP  plots  of  the  path  profile  or  contour 
plots  over  the  geographic  area. 

Figure  B-2  gives  a more  detailed  description  of  the  primary  overlay 
(CSPM)  execution  flow.  Program  CSPM  initially  issues  the  starting  dayfile 
messages  and  calls  a system  routine  to  return  the  starting  CP  time  which 
is  used  later  with  the  final  CP  time  to  calculate  the  execution  time  for 
the  run.  Next,  CSPM  calls  the  INPUT  overlay  to  read  the  user  input  data 
cards  and  check  for  a restart  file.  If  any  fatal  errors  are  detected  during 
input  processing,  CSPM  displays  the  number  of  errors  in  a dayfile  message 
and  aborts  the  run.  If  no  errors  are  found,  CSPM  checks  the  restart  flag 
to  determine  if  the  LOSSMOD  is  called  to  extract  path  profiles  and  generate 
ACR  and  statistical  data  for  a single  path  or  over  a geographic  area. 

If  fatal  errors  are  encountered  during  the  profile  extraction  due  to  missing 
data  on  the  TOPO  data  base,  CSPM  will  abort  the  run.  When  all  the  data 
have  been  calculated  and  written  to  file  LOSS,  CSPM  calls  the  OUTPUT  overlay 
to  read  the  LOSS  file  and  produce  the  requested  tabulated  and  plotted  output. 

If  this  is  a restart  run,  OUTPUT  reads  the  data  from  file  RESTRT,  since 
no  new  LOSS  file  has  been  generated.  When  all  the  requested  output  has 
been  produced,  CSPM  issues  the  final  dayfile  messages  to  signify  that  the 
run  completed  normally.  When  CSPM  has  finished  executing,  it  is  the  user's 
responsibility  to  dispose  the  PLOT  file  to  the  CALCOMP  plotter  and  (save 
the  LOSS  file  for  use  in  future  restart  runs) . Table  B-1  lists  the  files 
that  are  defined  by  CSPM  and  gives  a brief  description  of  the  type  of  data 
contained  in  each.  The  record  stirjctures  for  files  LOSS  and  CNTR  are  described 
in  Tables  B-2  and  B-3.  CSPM  also  defines  two  comnon  blocks  that  are  used 
throughout  the  program.  The  variedjles  in  common  block  /OPTIONS/  are  described 
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Figure  B-2 • Process  diagram  for  overlay  (0,0)  - CSPM. 
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Table  B-1.  Files  Used  by  CSPM 


File  Name 

Logical  Unit  Number 

Description 

INPUT 

5 

INPUT  file  for  the  CSPM  contains  the  1 

user  data  cards.  j 

OUTPUT 

6 

OUTPUT  file  for  the  CSPM  is  used  for 

all  printed  output.  ! 

PLOT 

7 

PLOT  file  contains  CALCCMP  plotter  ; 

commands  to  produce  profile  and  i 

contour  plots.  Each  plot  is  a i 

separate  file.  This  file  should 

be  disposed  to  a 12  inch  CALCOMP 

plotter  after  the  CSPM  run  is  i 

finished. 

RESTFT 

8 

LOSS  file  created  in  a previous  j 
run.  The  data  on  this  file  is  j 
read  when  the  restart  option  is  j 
selected.  j 

LOSS 

9 

LOSS  file  contains  topographic  path 
profile  data  and  ACR  and  statistical 
data  for  a single  path  or  a geo- 
graphic area.  These  values  are  used 
to  produce  tabulated  and  plotted  path 
profiles  and  to  calculate  LB,  PD, 

RSL,  and  REL  for  tabulated  printouts 
and  contour  plots. 

CNTR 

10 

CNTR  is  an  intermediate  scratch  file 
used  to  store  values  of  basic  trans- 
mission loss  and  communication 
reliability  for  a geographic  area 
for  the  contour  routines  to  use. 
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Table  B-2.  Record  Structure  for  File  LOSS 


File  LOSS  contains  three  types  of  records: 

OPTIONS  - the  values  in  common  block  /OPTIONS/  are  written  to  file 
LOSS  as  the  first  record.  If  this  file  is  used  later 
with  the  restart  option,  all  the  options  and  data  values 
used  in  the  initial  run  are  available. 

PROFILE  - contains  topographic  elevation  data  for  the  path  profile 
along  with  the  path  parameter  information. 

ACR  DATA  - contains  the  eight  ACR  and  statistical  values  needed  to 
calculate  basic  treuismission  loss  for  each  data  point 
along  the  path.  The  bearing,  path  length,  and  number  of 
data  points  on  the  path  also  are  given. 

For  a single  path,  the  LOSS  file  will  contain  only  the  OPTIONS  record 
followed  by  a PROFILE  record  for  the  path  and,  optionally,  an  ACR  DATA 
record  if  LB  has  been  requested  for  the  path.  VThen  a geographic  area  is 
being  analyzed,  the  paths  (map  radials)  are  processed  in  clockwise  order 
from  north.  If  the  radial  being  used  has  been  selected  for  a tabulated 
or  plotted  path  profile,  the  PROFILE  record  will  be  written  to  file  LOSS. 
If  tabulated  data  values  (LB,  PD,  RSL,  PEL)  have  been  requested  for  this 
radial  or  contour  plots  have  been  requested,  an  ACR  DATA  record  will  be 
written  to  file  LOSS  following  the  PROFILE  record.  Thus,  for  LOSS  files 
created  under  the  MAP  option,  the  PROFILE  and  ACR  DATA  records  are  inter- 
spersed depending  on  the  other  options  selected. 

The  structures  of  these  records  are  shown  below: 


Record  Type  Word  Number 


Description 


PROFILE 


Path  profile  "records  contain  the 
topographic  elevation  data  and  path  param- 
eters necessary  to  produce  tabulated  and 
plotted  path  profiles. 

1(N)  Number  of  words  to  follow  in  the  record. 

2 (FLAG)  Set  to  -1  to  flag  this  record  as  a 

profile. 

3 (TYPE)  Set  to  -1  if  this  profile  is  to  be 

plotted  only,  to  0 if  it  is  to  be 
tabulated  only,  and  to  1 if  both 
tabulated  and  plotted  output  are 
requested. 
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I’able  D-2.  (Continued) 


Record  Type  Word  tiumber  Description 

4 (AZIM)  Bearing  of  the  path  in  decimal  degrees, 

clockwise  from  north. 

5 (DLTH)  Terrain  irregularity  in  meters  for  the 

path. 

6 (ENS)  Surface  refractivity  for  the  path. 

7 (HSE(l))  Effective  surface  elevation  in  meters  at 

the  transmitter. 

8 (HSE(2))  Effective  surface  elevation  in  meters  at 

the  receiver. 

9 (DL(1))  Distance  in  meters  from  the  transmitter 

to  its  horizon. 

10  (DL(2))  Distance  in  meters  from  the  receiver  to 

its  horizon. 

11  (HL(1))  Transmitter  horizon  elevation  in  meters. 

12  (HL(2))  Receiver  horizon  elevation  in  meters. 

13  (DIST)  Path  length  in  meters. 

14  (XI)  Distance  increment  in  meters  between 

profile  points. 

15  (EP)  Number  of  points  on  the  profile  including 

the  point  at  the  transmitter  but  not  the 
point  at  the  receiver. 

16  to  K-15  NP+1  words  containing  elevations  in  meters 

for  each  profile  point. 

ACP,  DATA  ACR  records  contain  the  attenuation  and 

statistical  values  necessary  to  calculate 
LB,  PD,  RSL,  and  PEL  for  the  data  points 
on  the  path. 

1 (N)  Number  of  words  to  follow  in  the  record. 

2 (FLAG)  Set  to  1 if  tabulated  data  values  are  to  | 

be  printed  for  this  path  and  to  0 if  data  i 

arc  only  to  be  used  for  contour  plots.  j 

I 
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able  B-2.  (Continued) 


Record  Type  Word  Number  Description 

3 (AZIM)  Bearing  of  the  path  in  decimal  degrees, 

clockwise  from  north. 

4 (STDIST)  Distance  in  kilometers  from  the  trans- 

mitter to  the  first  data  point  on  the 
path. 

5 (PATHDST)  Path  length  in  meters. 

fi  (DLTH)  Terrain  irregularity  in  meters  for  the 

path. 

7 (EKE)  Surface  refractivity  for  the  path. 

8 (I'lUM)  Number  of  data  points  on  this  path  for 

which  a set  of  ACR  values  follow. 

A set  of  eight  values  follows  for  each  data  point  on  the 
path  (repeat  for  K=1,NUM) 

(K-l)*8+9  Free  space  loss  in  dB  for  the  path. 

(XLBF) 

(K-l)*6+10  A predicted  reference  value  of  attenuation 

(ACRV)  below  free  space  in  dB. 

(K-1) *8+11  The  deviation  of  calculated  reference 

(VMD)  attenuation  from  yearly  median  attenuation  - 

a function  of  climatic  region  and  effective 
distance. 

(K-1) *8+12  Standard  deviation  of  prediction  error. 

(SGC) 

(K-1) *8+13  Standard  deviation  of  predicted  attenuation 

due  to  location  variability. 

(K-l*8+14  Standard  deviation  of  predicted  attenuation 

(SGTM)  for  time  variabilities  less  than  the  median 

value. 

(ic-l)*e  + 15  Standard  deviation  of  predicted  attenuation 

(SC'T)  for  time  variabilities  greater  than  the 

median  value. 


Table  B-2.  (Continued) 


Record  Type  Word  Number  Description 

(K-l)*e+16 

(SGTD) 

Standard  deviation  of  predicted  attenuation 
for  small  time  variabilities,  such  as  are 
associated  with  the  ducting  phenomenon. 

Table  B-3. 

Record  Structure  for  File  CNTR 

Progreim  OUTPUT  writes  a record  in  the  following  format  to  file  CNTR  for 
each  path  (map  radial)  in  the  geographic  area  being  analyzed.  The  routine 
CONTOUR  reads  these  data  and  uses  them  to  generate  the  requested  contour 
plots . 

Word  Number 

Description 

1 (STDIST) 

Distance  in  kilometers  from  the  transmitter  to 
the  first  data  point  on  the  path. 

2 (AZIM) 

Bearing  of  the  radial  in  decimal  degrees 
clockwise  from  north. 

3 (NP) 

Number  of  data  points  on  the  path. 

4 to  NP+3 

(XL0SS(V1,J) ,J=1,NP) 

Basic  transmission  loss  values  for  all  data 
points  for  variability  set  Vl. 

NP+4  to  2*NP+3  . 
(XLOSS(V2,J) ,J=1,NP) 

Basic  transmission  loss  values  for  all  data 
points  for  variability  set  V2. 

2*NP+4  to  3*NP+3 
(XL0SS(V3,J) ,J=1,NP) 

Basic  transmission  loss  values  for  all  data 
points  for  variability  set  V3. 

3*NP+4  to  4*NP+3 
(XREL(V1,J) ,J=1,NP) 

Communication  reliability  values  for  all  data 
points  for  variability  set  VI. 

4*NP+4  to  5*NP+3 
(XREL(V2,J) ,J“1,NP) 

CoRinunication  reliability  values  for  all  data 
points  for  variability  set  V2. 

5*NP+4  to  6*NP+3 
(XREL(v3, J) ,J=1,NP) 

Communication  reliability  values  for  all  data 
points  for  variability  set  V3. 
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Table  B-4.  Variables  in  Common  Block  /OPTIONS/ 


i 

i 


Variable  Name 

TITLE 

IDATE 

ITIME 

MAP 

BLATl,  BLAT2, 
BLONl,  BLON2 

PRO 

PTAB 

PLOT 

TABL 

LB 

PD 

RSL 

REL 

XLAT 

XLON 

RLAT 

RLON 

XLOC 


Description 

Title  for  the  CSPM  run  (60  characters  maximum) . 
Date  of  the  CSPM  run. 

Time  of  the  CSPM  run. 

Flag  set  to  0 for  single  path  and  to  1 for  map 
option . 

The  latitude  and  longitude  boundaries  in  decimal 
degrees  of  the  map  area. 

Number  of  path  profiles  to  be  plotted. 

Number  of  path  profiles  to  be  ted^ulated. 

Flag  set  to  1 if  contour  plots  are  requested  and 
to  0 if  no  plots  are  to  be  produced. 

Number  of  paths  for  which  LB,  PD,  RSL,  and  REL 
will  be  tabulated. 

Flag  set  to  t if  basic  transmission  losses  are  to 
be  tabulated  or  plotted  and  to  0 otherwise. 

Flag  set  to  1 if  power  density  values  are  to  be 
tabulated  or  plotted  and  to  0 otherwise. 

Flag  set  to  1 if  received  signal  level  values  are 
to  be  tabulated  or  plotted  and  to  0 otherwise. 

Flag  set  to  1 if  communication  reliability  values 
are  to  be  tedsulated  or  plotted  and  to  0 
otherwise. 

Transmitter  latitude  in  decimal  degrees  north. 
Tr2msmitter  longitude  in  decimal  degrees  east. 
Receiver  latitude  in  decimal  degrees  north. 
Receiver  longitude  in  decimal  degrees  east. 

Name  of  tr^msmitter  location. 
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Variable  Name 

RL(X 

POL 

CLIM 

GND 

FREQ 


j PWR 

I 

I 


k 

[ 

i 

( 

I 


HGl 

HG2 

G1 

G2 

BPS 


SGM 

ENS 


RSLTH 

PUTIW 

HGIF 

HG2F 

VI,  V2,  V3 

NV 

RQ 


Table  B-4.  (Continued) 


Description 

Name  of  receiver  location. 
Polarization  of  transmitting  antenna. 
Climate  descriptor. 

Ground  descriptor. 

Carrier  frequency  in  MHz. 

Transmitter  power  in  dBm. 

Transmitter  antenna  height  in  meters. 
Receiver  antenna  height  in  meters. 
Transmitter  antenna  gain  in  dBi. 


Receiver  antenna  gain  in  dBi. 

Ground  dielectric  constant. 

Ground  conductivity. 

If  negative,  ABS(ENS  = NS  (surface  refractivity) , 
and  if  positive,  ENS  = HO  (surface  refractivity 
referenced  to  mean  sea  le'^el) . 


Received  signal  level  threshold  in  dBm. 

Transmitter  power  in  watts. 

Transmitter  antenna  height  in  feet. 

Receiver  antenna  height  in  feet. 

Three  sets  of  variability  statistics  (Particular 
QT,  QL,  Q values  from  one  set  of  variability 
statistics. ) 

Number  of  sets  of  variability  statistics  to  be 
used  in  the  CSPM  run. 

Contour  levels  specified  for  REL  plots. 
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Table  13-4.  (Continued)  ) 


I 


Vari^d5le  Name 
NRQ 

RESTART 

NERR 

XCDHS 

XCDD 

RCDMS 

RCDD 

MCmS 


Description 


Number  of  contour  levels  to  be  used  for  REL  plots. 

Flag  set  to  1 if  the  restart  option  has  been 
selected  and  to  0 if  not. 

Fumber  of  fatal  errors  found. 

Transmitter  coordinates  in  degrees,  minutes, 
seconds  display  code  format 
(XXX,XX,XX,D/XXX,>:X,XX,D)  . 

Transmitter  coordinates  in  decimal  degrees  display 
code  format  (XXX.XX,D/XXX.XX,D) . 

Receiver  coordinates  in  degrees,  minutes, 
seconds  display  code  format 
(XXX,XX,XX,D/XXX,XX,XX,D) . 


Receiver  coordinates  in  decimal  degrees  display 
code  format  (XXX.XX,D/XXX.XX,D) . 

Map  bounds  coordinates  in  degrees  minutes, 
seconds  display  code  format 

( XXX , XX , XX , D/XXX , XX , XX , D/XXX , XX , XX , D/XXX , XX , XX , D ) . 


MCDD 


Map  bounds  coordinates  in  decimal  degrees  display 
code  format  (XXX. XX, D/XXX. XX, D/XXX. XX, D/XXX. XX. D) . 


Teible  B-5.  Variables  in  Common  Block  /LOSSMOD/ 


Variable  Neune 

Description 

XMAX 

Distance  in  meters  from  the  transmitter  to 
east  longitude  boundary  of  the  map  area. 

the 

XMIN 

Distance  in  meters  from  the  transmitter  to 
west  longitude  boundary  of  the  map  area. 

the 

YMAX 

Dist2uice  in  meters  from  the  transmitter  to 
north  latitude  boundary  of  the  map  area. 

the 

YMIN 

Distance  in  meters  from  the  transmitter  to 
north  latitude  boundary  of  the  map  area. 

the 

Bl 

Bearing  in  decimal  degrees  to  the  NE  corner 
the  map  area. 

of 

B2 

Bearing  in  decimal  degrees  to  the  SE  corner 
the  map  area. 

of 

B3 

Bearing  in  decimal  degrees  to  the  SW  corner 
the  map  area. 

of 

B4 

Bearing  in  decimal  degrees  to  the  IW  corner 
the  map  area. 

of 

ZZD 

Standard  normal  deviate  for  time  variability  in 
attenuation  at  small  probabilities  (often 
associated  with  the  phenomonen  of  ducting) . 

in  Table  B-4  and  those  in  /LOSSMOD/  in  Table  B-5.  These  data  values  and 
option  flags  are  referenced  in  many  routines  in  each  overlay. 


The  execution  flow  of  the  (1,0)  overlay  is  shown  in  Figure  B-3.  INPUT, 
first,  initializes  necessary  flags  and  data  values  and  then  reads  the  title 
card  for  the  run.  It  then  checks  for  a RESTRT  file  by  attempting  to  read 
the  OPTIONS  record  from  it.  If  the  read  is  successful,  the  /OPTIONS/  common 
block  will  be  reset  to  the  values  used  when  the  RESTRT  file  was  created. 

Input  then  sets  the  restart  flag  and  identifies  the  options  that  may  not 
be  respecified  for  a restart  run  as  illegal.  If  an  EOF  is  read,  INPUT 
sets  its  flags  for  an  initial  run.  In  either  case,  INPUT  then  proceeds 
to  read  the  options  from  the  user  input  data  cards.  The  NXTVAR  routine 
is  called  to  read  the  data  cards  as  needed  and  return  the  next  syntactic 
elements  and  numeric  values.  It  is  described  in  Table  B-6.  As  each  new 
option  keyword  is  read,  INPUT  checks  it  against  a list  of  legal  keywords. 

It  also  checks  to  see  if  that  keyword  has  been  read  previously  during  this 
run,  if  it  is  illegal  with  restart  (if  this  is  a restart  run),  or  if  it 
conflicts  with  another  option  already  processed.  Once  this  preliminary 
checking  is  complete,  the  program  branches  to  the  section  of  INPUT  that 
reads  the  values  given  for  the  specified  option.  Here  NXTVAR  again  is 
called  to  return  the  fields  that  comprise  the  values  of  the  option.  The 
syntax  of  the  option  is  checked  in  addition  to  the  value  given.  Numeric 
values  are  checked  to  be  within  specified  limits  and  keywords  must  be  one 
of  a legal  list  of  options. 

Once  all  the  option  keywords  have  been  processed,  INPUT  cross-checks 
all  the  selections  to  be  sure  that  no  information  is  missing  or  that  con- 
flicting options  have  not  been  chosen.  INPUT  issues  two  levels  of  error 
diagnostics.  Warning  messages  are  given  when  execution  can  proceed,  but 
the  user  should  be  made  aware  of  some  possibly  erroneous  input  data,  such 
as  an  option  specified  twice.  Fatal  error  diagnostics  are  issued  when 
the  input  data  is  incomplete,  in  conflict,  out  of  range,  or  contains  syntax 
errors.  Fatal  errors  will  cause  the  CSPM  run  to  be  aborted.  The  PRINT 
routine  is  called  to  print  the  input  data  summary  which  lists  all  the  options 
and  values  selected  by  the  user  or  by  default  that  will  be  used  for  this 
CSPM  run.  Finally,  if  a restart  run  has  not  been  requested,  the  OPTIONS 
record  is  written  to  file  LOSS  and  control  is  returned  to  the  CSPM. 
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Figure  B-3.  Continued 


Table  B-6.  Variables  in  Common  Block  /KXTVAE/ 


Variad^le  Fame 

Description 

IVAR 

Alphanxameric  keyword  read  is  returned  here  in  AlO 
format. 

XNUM 

Decoded  real  number  is  returned  here.  A -1.0  is 
returned  in  the  case  of  a null  field  where  the 
following  delimiter  is  a , or  a ) . 

IDEL 

The  delimeter  character  read  is  returned  here.  It 
must  be  one  of  the  characters  =,/()  or  -1  to 
signify  end-of-card. 

lERR 

The  error  flag  is  set  to  1 if  an  error  is  found  by 
MXTVAR,  -1  if  an  error  is  found  by  the  calling 
routine,  99  for  EOF  on  INPUT,  and  0 if  no  errors. 
If  a syntax  error  is  detected,  the  rest  of  the 
data  card  is  skipped. 

Table  B-7.  Varicibles  in  Common  Block  /PFILE/ 


Varicible  Ncune 

Description 

DIST 

Path  length  in  meters. 

XI 

Distance  increment  in  meters  between  profile 
points. 

ENP 

Number  of  profile  points  (real  value)* 

NP 

Number  of  profile  points  (integer  value)* 

PFL 

Array  containing  the  elevation  in  meters  at  each 
profile  point. 

*The  number  of  profile  points  given  here  includes  the  point  at  the 
transmitter  but  does  not  include  the  point  at  the  receiver.  To 
reference  the  elevation  at  the  transmitter  site,  use  PFL(l)  and  for 
the  elevation  at  the  receiver  site,  use  PFL(NP+1) . The  values  of 
ENP  and  NP  are  equal,  ENP  is  saved  here  to  avoid  the  conversions  of 
the  integer  NP  value  to  a real  value  to  be  used  in  calculations. 


The  flow  diagram  for  the  (2,0)  overlay  (figure  B-4)  explains  the  process 
of  extracting  path  profiles  and  calculating  ACR  and  statistical  values  with 
the  Longley-Rice  propagation  prediction  routines  for  a single  path  or  points 
throughout  a geographic  area.  The  variables  in  common  blocks  /PFIliE/  and 
/GEODYS/,  described  in  Tables  B-7  and  B-8,  are  used  to  communicate  with 
the  path  profile  extraction  routines.  Common  blocks  /EVPARS/,  /PHPARS/, 
and  /SYPARS/  contain  variables  used  to  describe  characteristics  of  the  system 
and  are  used  by  the  routines  that  calculate  path  parameter  information. 

These  common  blocks  are  described  in  Tables  b-9,  B-10,  and  B-11.  Notice 
that  some  of  the  same  values  that  appear  in  these  common  blocks  are  repeated 


in  the  common  blocks  used  with  the  Longley-Rice  routines,  /UlINP/,  /LRDAS/, 
/LRVAR/,  and  /LRRTH/  which  are  described  in  Tables  B-12,  B-13,  B-14,  and 
B-15.  Two  separate  cases  are  considered  in  LOSSMOD:  one  for  a single  path, 
and  emother  for  a geographic  area. 

If  a single  path  has  been  specified,  LOSSMOD  first  calls  PFLTPO  to 
extract  path  ele/ations  from  the  topographic  data  base.  Next,  the  path 
parameters  for  this  path  are  calculated  from  the  input  data  and  the  profile 
values  by  QPARS.  Next,  LOSSMOD  sets  a flag  to  tell  the  OUTPUT  routine  whether 
the  path  profile  is  to  be  plotted,  tabulated,  or  plotted  and  tabulated. 

Next,  the  PROFILE  record  containing  the  output  type  flag,  the  path  param- 
eter values,  and  the  profile  elevations  is  written  to  file  LOSS.  (See  Table 
B-2  for  the  structure  of  this  record) . LOSSMOD  then  checks  to  see  if  basic 
tramsmission  loss  has  been  requested  for  the  path.  If  so,  it  tr£msfers 
the  necessary  data  values  from  the  path  parameter  common  blocks  to  those 
used  by  the  Longley-Rice  routines.  These  comnon  blocks  have  not  been  combined, 
because  the  Longley-Rice  routines  can  be  replaced  more  easily  with  new  versions 
if  they  are  not  modified  in  the  CSPM.  The  climate-dependent  variables  are 
initialized  by  calling  QVARLR  with  the  climate  code  specified  in  the  /OPTIONS/ 
comnon  block.  Next,  QLRPP  is  used  to  initialize  necessary  variables  for 
the  Longley-Rice  routines  when  used  in  the  point-to-point  mode.  It  also 
calls  QACRLR  to  initialize  the  constemts  and  the  coefficients  for  the  dif- 
fraction region.  Next,  LOSSMOD  calls  QQVARLR  with  the  total  path  length 
to  calculate  the  ACR  and  statistical  values  for  the  path.  These  values 
are  returned  in  common  block  /LRVAR/.  This  information,  needed  to  calculate 
LB  for  the  path  is  written  to  file  LOSS  in  an  ACR  DATA  record  (see  Table  B-2) , 
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Figure  B-4.  Process  diagram  for  overlay  (2,0)  - LOSSMOD. 
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Table  B-8 


Variables  in  Common  Block  /EVPABS/ 


Variable  Meime 

Description 

HSYS 

The  average  of  the  transmitter  and  receiver  site 
elevations  above  mean  sea  level  in  meters. 

ENO 

NO  is  the  surface  refractivity  referenced  to  mean 
sea  level. 

ENS 

NS  is  the  surface  refractivity. 

RKFE 

The  effective  earth  radius  coefficient. 

GMFE 

The  curvature  of  the  effective  earth. 

EPS 

Ground  dielectric  const2mt. 

SGM 

Ground  conductivity. 

DLTH 

Interdecile  remge  of  the  terrain  elevation  in 
meters.  Delta  H is  a measure  of  the  terrain 
irregularity  along  a path. 

HMSL(l) 

Transmitter  site  elevation  above  mean  sea  level  in 
meters . 

fIMSL(2) 

Receiver  site  elevation  above  mean  sea  level  in 
meters. 
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Table  B-9.  Variables  in  Common  Block  /PHPARS/ 


VarieQ^le  Neune 

DS 

HSE(l) 

HSE(2) 

KSE(3) 

HS(1) 

HS(2) 

HE(1) 

HE(2) 

ISO 

HL(1) 


Description 


Path  length  in  meters.  j 

Value  returned  from  FFGND  is  the  effective  surface  ] 

elevation  in  meters  for  the  path  from  the  trans-  j 

mitter  to  its  horizon.  It  is  determined  by  ] 

averaging  the  middle  80%  of  the  profile  elevations  | 

for  the  path.  j 

I 

1 

Value  returned  from  FFGND  is  the  effective  surface  j 

elevation  in  meters  for  the  path  from  the  receiver  j 

to  its  horizon.  i 

Used  by  FFGND  to  calculate  the  effective  surface  i 

elevation  for  the  path  from  the  transmitter  j 

horizon  to  the  receiver  horizon  for  double  j 

horizon  paths . | 

! 

The  total  height  of  the  transmitter  above  mean  sea  j 

level  in  meters  calculated  by  summing  the  j 

structural  antenna  height  of  the  transmitter  and 
the  site  elevation. 

The  total  height  of  the  receiver  above  mean  sea 
level  in  meters  calculated  by  summing  the  i 

structural  antenna  height  of  the  receiver  and  the  ' 

site  elevation.  ’ 

The  effective  antenna  height  in  meters  of  the 
transmitter  which  is  the  meucimum  of  the  structural 
antenna  height  of  the  transmitter  or  the 
difference  between  HS(1)  and  HSE(l). 

The  effective  emtenna  height  in  meters  of  the 
receiver  which  is  the  maximum  of  the  structural 
antenna  height  of  the  receiver  or  the  difference 
between  MS (2)  and  HSE(2). 

Flag  set  to  0 if  path  is  line-of-sight,  1 if 
double  horizon,  and  -1  if  single  knife  edge.  This 
information  presently  is  not  used  in  the  Longley- 
Rice  model. 


Transmitter  horizon  elevation  in  meters. 


J 
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Table  B-9 . (Continued) 


Variadsle  Name 

Description 

HL(2) 

Receiver  horizon  elevation  in  meters. 

DL(1) 

Distance  from  the  transmitter  to  its  horizon  in 
meters. 

DL(2) 

Distance  from  the  receiver  to  its  horizon  in 
meters. 

THE(l) 

Transmitter  horizon  elevation  angle  in  radians. 

THE (2) 

Receiver  horizon  elevation  angle  in  radians. 

THEA 

Sum  of  the  horizon  elevation  angles. 

Table 

B-10.  Variables  in  Common  Block  /SYPARS/ 

Variable  Name 

Description 

FMHZ 

Carrier  frequency  in  MHz. 

WN 

Wave  number. 

JPOL 

Polarization  of  the  transmitting  antenna,  set  to 

0 for  horizontal  or  1 for  vertical  polarization. 

HG(1) 

Transmitter  structural  antenna  height  in  meters. 

HG(2) 

Receiver  structural  euitenna  height  in  meters. 

r 

I 

- - - 


Variable  Neune 


Description 


DH 

Terrain  irregularity  in  meters. 

I 

FMHZ 

Carrier  frequency  in  MHz. 

HG(1) 

Transmitter  structural  antenna  height  in  meters. 

i 

HG(2) 

Receiver  structural  antenna  height  in  meters. 

HMSL(l) 

Transmitter  site  elevation  in  meters. 

^IMSL(2) 

Receiver  site  elevation  in  m.eters. 

EPS 

Ground  dielectric  constant. 

SGM 

Ground  conductivity. 

ENO 

NO  is  the  surface  refractivity  referenced  to  me2m 
sea  level. 

ENS 

NS  is  the  surface  refractivity. 

KST(l) 

Unused. 

1 

i 

KST(2) 

Unused. 

1 

j 

KLIM 

Code  number  for  climate  descriptor. 

IPOL 

Polarization  of  the  transmitting  antenna,  set  to  0 
for  horizontal  and  to  1 for  vertical  polarization. 

Table  B-12.  Variables  in  Coranon  Block  /LRDAS/ 


Variable  Name 

Description 

HE(1) 

The  effective  antenna  height  in  meters  of  the 
transmitter;  value  is  the  meucimum  of  the  structural 
antenna  height  of  the  transmitter  and  the  dif- 
ference between  the  total  height  of  the  transmitter 
above  mean  sea  level  and  the  effective  surface 
elevation  for  the  transmitter. 

HE  (2) 

The  effective  antenna  height  in  meters  of  the 
receiver;  value  is  the  meucimum  of  the  structural 
antenna  height  of  the  receiver  and  the  difference 
between  the  total  height  of  the  receiver  above 
mean  sea  level  and  the  effective  surface  elevation 
for  the  receiver. 

THE(l) 

Transmitter  horizon  elevation  angle  in  radians. 

THE (2) 

Receiver  horizon  elevation  angle  in  rac  ins. 

THA 

Sum  of  the  horizon  elevation  angles. 

HZ(1) 

Transmitter  effective  antenna  height  that  includes 
a complex  factor  to  compensate  for  ground  wave 
effects. 

HZ{2) 

Receiver  effective  antenna  height  that  includes  a 
complex  factor  to  compensate  for  ground  wave 
effects. 

DL{1) 

Distance  from  the  transmitter  to  its  horizon  in 
meters. 

DL(2) 

Distance  from  the  receiver  to  its  horizon  in 
meters. 

DLA 

Sum  of  the  distances  in  meters  to  the  transmitter 
and  receiver  horizons. 

DLS(l) 

Distance  from  the  transmitter  to  its  horizon  in 
meters  over  a smooth  earth. 

DLS(2} 

Distance  from  the  receiver  to  its  horizon  in 
meters  over  a smooth  earth. 

DLSA 

1 

Sum  of  the  distances  in  meters  to  the  transmitter 
and  receiver  horizons  over  a smooth  earth. 
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Table  B-13.  VariaOoles  in  Corimon  Block  /LRVAR/ 


Variable  Name 

Description 

XLBF 

Free  space  loss  in  dB. 

ACRV 

A predicted  reference  value  of  attenuation  below 
free  space  in  dB. 

VMD 

The  median  time  deviation  above  reference  loss. 

SGC 

Steuidard  deviation  of  prediction  error. 

SGL 

Standard  deviation  of  predicted  attenuation  due  to 
location  variability. 

SGTM 

Standard  deviation  of  predicted  attenuation  for 
time  variabilities  less  than  the  median  value. 

SGTP 

Standard  deviation  of  predicted  attenuation  for 
time  variabilities  greater  than  the  median  value. 

SGTD 

Standard  deviation  of  predicted  attenuation  due  to 
time  variabilities  for  the  ducting  range. 

ZD 

Standard  normal  deviate  of  predicted  attenuation 
for  small  time  variabilities,  such  as  are 
associated  with  the  ducting  phenomenon. 

YD 


The  ducting  range  breakpoint,  ZD*SGTP. 


Variable  Name 


Description 


The  effective  earth's  radius. 

The  effective  earth's  radius  coefficient. 

The  inverse  of  the  effective  earth's  radius. 

The  surface  transfer  impedance  (complex) . 

Wave  number. 

A filtering  constant  for  diffraction  weights.  For 
known  paths,  such  as  in  the  CSPM,  it  is  set  to  10. 

A normalized  earth's  radius. 


Tadsle  B-15.  Variables  in  Common  Block  /LBREL/ 


Variable  Name 


Description 


STDIST 


XLOSS 


Distance  in  kilometers  from  the  transmitter  to  the 
first  data  point  along  a map  radial. 

Bearing  of  the  path  in  decimal  degrees  clockwise 
from  north. 

Number  of  data  points  on  this  path. 

Array  containing  basic  transmission  loss  values 
at  each  data  point  on  the  path  for  all  three  sets 
of  variability  statistics. 

Array  containing  communication  reliability  values 
for  each  data  point  on  the  path  for  all  three  sets 
of  variability  statistics. 

Array  used  for  AGR  data  values. 


emd  control  is  returned  to  CSPM. 

The  procedure  followed  for  calculating  values  for  a geographic  area 
, repeats  the  process  used  for  a single  path  for  each  data  point  in  the  area. 

LOSSMOD  calls  QVARLR  at  the  beginning  to  initialize  climate-dependent  variables 
for  the  map  area.  The  method  used  to  select  the  data  points  over  the  geographic 
area  is  described  in  Appendix  F.  Basically,  it  involves  determining  how 
many  radials  will  be  required  to  adequately  cover  the  area  and  then  using 
; the  system  described  in  Table  F-1  to  locate  data  points  along  each  radial. 

* LOSSMOD  determines  the  number  of  radials  to  be  used  by  calculating  the  maximum 

( 

|l  distance  from  the  transmitter  to  any  point  on  the  map  boundary.  The  number 

i of  radials  is  then  equal  to  the  largest  power  of  two  that  is  less  than  this 

I maximum  distance.  Once  the  number  of  radials  has  been  determined,  the  angular 

I increment  between  radials  can  be  calculated.  It  is  also  necessary  to  consider 

[ the  case  in  which  the  transmitter  is  not  interior  to  the  area  boundary, 

I 

I but  lies  on  a side  or  in  a comer  of  the  boundary.  In  this  case,  radials 

' are  not  generated  for  a full  360°,  but  the  initial  ^tnd  final  bearings  are 

> restricted,  so  that  all  radials  generated  will  be  within  the  area.  LOSSMOD 

^ loops  through  all  the  radials  to  be  used  for  the  area.  A radial  may  be 

I skipped  if  no  contour  plots  have  been  requested,  and  the  data  on  the  radial 

f 

I are  not  to  be  tcdsulated.  If  contour  plots  have  been  selected,  data  for 

each  radial  must  be  calculated  and  written  to  file  LOSS,  so  that  the  entire 
geographic  area  will  be  covered  with  data  points. 

Each  radial  is  considered  to  be  like  a series  of  single  paths,  one 
from  the  transmitter  to  each  data  point  along  the  radial.  The  distance 
between  the  transmitter  cuid  the  first  data  point  to  be  used  on  that  radial 
is  determined  according  to  the  method  shown  in  Table  F-1.  To  insure  that 
the  contour  plots  will  be  accurate,  close  to  the  area  boundary,  one  extra 
data  point  outside  the  boundary  is  calculated  for  each  radial.  This  point 
is  determined  by  using  QRTHETA  to  find  the  point  where  the  radial  intersects 
the  boundary,  to  find  the  distemce  from  the  tr2msmitter  to  that  point,  and 
then  to  round  that  dlstemce  to  the  nearest  kilometer  and  find  the  point 
along  the  radial  at  that  rounded  distance  from  the  tremsmitter  using  QRDESIC. 
Data  points  then  will  be  calculated  at  one  kilometer  intervals  along  the 
radial  from  the  starting  point  to  the  final  point  outside  the  box.  PFLTPO 
is  called  to  extract  the  path  profile  for  the  path  from  the  transmitter 
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to  the  final  data  point.  The  point  spacing,  for  path  elevations,  is  fixed 
at  250  m,  since  the  path  length  is  always  an  even  number  of  kilometers. 

Thus,  the  profile  can  be  used  in  sections  (four  points  to  a kilometer)  to 
provide  the  profile  data  for  any  one  of  the  subpaths  from  the  transmitter 
to  a data  point  without  calling  the  profile  extraction  routines  each  time. 

QPARS  is  used  to  calculate  the  path  parameters  for  the  entire  radial  path. 

If  the  profile  for  this  radial  is  to  be  tcibulated  or  plotted,  a PROFILE 
record  is  written  to  file  LOSS.  LOSSMOD  then  loops  through  all  the  subpaths 
for  this  radial  calling  QPARX  to  calculate  the  path  parameters  for  the  profile 
values  between  the  transmitter  and  each  data  point.  Next,  QLRPP  is  called 
to  initialize  the  diffraction  range  coefficients  for  the  path  and  to  set 
values  for  the  propagation  loss  calculation.  Finally,  QQVARLR  is  called 
with  the  distance  from  the  transmitter  to  the  data  point  to  calculate  the 
ACR  and  statistical  values  at  that  point.  These  are  stored  in  the  DATA 
array.  When  all  data  points  for  that  radial  have  been  processed,  the  ACR 
DATA  record  for  that  radial  is  written  to  file  LOSS.  This  process  continues 
until  the  data  values  have  been  calculated  for  all  the  map  radials.  The 
LOSS  file  then  contains  all  the  necessary  data  to  plot  or  tabulate  basic 
transmission  loss,  power  density,  received  signal  level,  or  communication 
reliability  for  the  geographic  area  and  control  returns  to  CSPM. 

Figure  B-5  explains  the  execution  flow  of  the  (3,0)  OUTPUT  overlay. 

When  an  initial  run  has  been  made  (not  a restart  run) , OUTPUT  rewinds  the 
LOSS  file  and  skips  the  OPTIONS  record  positioning  to  read  the  PROFILE  and 
ACR  DATA  records  that  follow.  If  this  is  a restart  run,  the  file  RESTRT 
is  read,  and  INPUT  already  will  have  read  the  OPTIONS  record,  so  no  positioning 
is  needed.  Next,  OUTPUT  calls  QERFI  to  calculate  the  standard  normal  deviates 
for  the  variability  sets  to  be  used.  It  also  calculates  the  power  density 
and  received  signal  level  consteuits  that  are  added  to  basic  transmission 
loss  to  calculate  power  density  2md  received  signal  level  values.  OUTPUT 
then  reads  each  record  from  file  LOSS  or  RESTRT  until  an  EOF  is  found. 

If  the  record  read  is  a PROFILE  record,  PRTAB  is  call^  to  tabulate  or  plot 
the  profile  data.  If  this  is  a restart  run,  all  PROFILE  records  are  skipped. 

If  the  record  read  is  an  ACR  DATA  record  and  the  run  is  for  a single  path, 
the  basic  transmission  loss  for  each  variability  set  specified  is  calculated 
and  printed,  and  control  returns  to  CSPM.  For  a geogr^hic  area,  the 
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ACR  DATA  record  contains  the  values  necessary  to  calculate  basic  transmission 


loss  for  every  data  point.  This  value  is  stored  in  the  array  XLOSS. 

If  communication  reliability  is  requested,  it  also  is  calculated  and  stored 
in  array  XREL  for  each  data  point  and  for  each  variability  set.  PRLOSS 
is  called  to  print  the  requested  basic  transmission  loss,  power  density. 


received  signal  level,  or  communication  reliability  values  if  tabulated 


data  have  been  requested  for  this  radial.  If  contour  plots  have  been 


requested,  the  basic  transmission  loss  and  communication  reliability  values 


are  written  to  file  CNTR  for  each  radial.  When  all  the  records  on  LOSS 


or  RESTRT  have  been  read,  OUTPUT  checks  the  plot  flag  to  see  if  contour 


plots  are  to  be  produced.  If  so,  CONTOUR  is  called  to  read  the  CNTR  data 


file  and  create  the  CALCOMP  plots.  When  control  is  returned  to  CSPM, 


the  final  dayfile  messages  are  issued  and  the  CSPM  run  is  complete. 

For  more  detailed  description  of  the  routines  involved  in  each  overlay. 


refer  to  documentation  given  in  section  B.2. 


B.2  Introduction 


Documentation  for  each  routine  used  in  the  CSPM  (except  the  CALCOMP 


plotting  software)  is  included  in  this  section.  Each  routine  is  described 


in  the  same  order  as  it  appears  on  the  UPDATE  program  library.  The  informa- 


tion given  for  each  routine  includes: 


(1)  a short  statement  of  purpose, 

(2)  a description  of  the  input  and  output  values, 

(3)  a list  of  routines  called. 


(4)  a list  of  calling  routines,  and 

(5)  a general  discussion  of  algorithms  used  and  execution  flow. 


The  routine  descriptions  that  follow  have  been  grouped  by  overlay  to  facili- 


tate coordinating  the  information  given  here  with  that  in  section  B.l. 


B.2.1  (0,0)  Overlay  Routines 


B.2. 1.1  CSPM 


Program  CSPM  is  the  main  controlling  routine  for  the  (0,0)  overlay. 

It  controls  execution  by  calling  the  secondary  overlays  as  they  are  needed 
and  issues  the  CSPM  dayfile  messages. 
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FILES 


The  files  used  by  CSPM  are  listed  in  Table  B-1.  CSPM  does  not  attach 
or  catalog  any  files.  The  user  must  include  the  appropriate  control  cards 
in  his  deck  to  attach  a RESTRT  file,  catalog  a LOSS  file,  or  dispose  a PLOT 
file  to  the  CALCOMP  plotter.  The  record  structures  for  files  LOSS  (RESTRT) 
and  CNTR  are  given  in  Tables  B-2  and  B-3. 

COMMON  BLOCKS 

CSPM  defines  two  common  blocks  that  are  used  to  store  options  and 
data  values  used  in  more  than  one  overlay.  Tables  B-4  and  B-5  describe 
the  variables  included  in  the  /OPTIONS/  and  /LOSSMOD/  common  blocks. 

ROUTINES  CALLED 

CSPM  calls  the  following  routines; 

INPUT 

LOSSMOD 

OUTPUT 

PROGRAM  DESCRIPTION 

CSPM  displays  information  about  the  number  of  fatal  errors  found, 
the  number  of  plots  generated,  and  the  execution  time  for  the  run  in  dayfile 
messages.  It  initially  calls  the  system  routine  SECOND  to  return  the 
starting  CP  time.  At  the  end  of  the  run,  it  again  calls  SECOND  to  obtain 
the  final  CP  time  and  then  displays  the  execution  time  for  the  run  in 
seconds  in  the  dayfile.  If  fatal  errors  are  encountered  during  input 

processing  or  topographic  data  extraction,  the  number  of  errors  is  displayed 

« 

and  the  run  aborted.  If  plotted  path  profiles  or  contour  plots  have  been 
requested,  the  number  of  plots  generated  by  this  CSPM  run  on  file  PLOT 
also  will  be  displayed.  When  the  run  is  ccxnplete,  the  final  dayfile  message 
is  issued  to  inform  the  user  that  the  CSPM  terminated  normally. 

CSPM  controls  the  execution  flow  by  calling  the  secondary  overlays 
in  the  following  order.  The  INPUT  (1,0)  overlay  is  called  to  read  the 
user  input  data  cards  and  check  for  a restart  file.  If  any  fatal  errors 
are  detected  during  input  processing,  CSPM  displays  the  number  of  errors 
in  a dayfile  message  and  aborts  the  run.  If  no  errors  are  found,  CSPM 
checks  the  restart  flag  to  determine  if  the  LOSS  (2,0)  overlay  should 

be  called.  If  a RESTRT  file  has  not  been  provided,  IXJSSMOD  is  called 
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to  extract  path  profiles  and  generate  ACR  and  statistical  data  for  a single 
path  or  data  points  over  a geographic  area.  If  fatal  errors  are  encountered 
during  the  profile  extraction  due  to  missing  data  on  the  topographic  data 
base,  CSPM  aborts  the  run.  When  all  the  data  have  been  calculated  and 
written  to  file  LOSS,  CSPM  calls  the  OUTPUT  (3,0)  overlay  to  read  the  LOSS 
file  and  produce  the  requested  tabulated  and  plotted  output.  If  this  is 
a restart  run,  OUTPUT  reads  the  data  from  file  RESTRT,  since  no  new  LOSS 
file  has  been  generated.  When  CSPM  has  finished  executing,  it  is  the  user's 
responsibility  to  dispose  the  PLOT  file  to  the  CALCOMP  plotter  and  save 
the  LOSS  file  for  use  in  future  restart  runs. 

B.2.1.2  HEADER 

Subroutine  HEADER  prints  the  two  line  page  heading  on  every  page  of 
the  CSPM  output. 

ARGUMENTS 

NLINES  - Number  of  lines  to  be  printed  by  the  calling  routine.  If  this 
number  is  greater  than  the  number  of  lines  left  on  the  page,  a 
page  eject  is  issued  and  the  page  header  is  printed.  A new  page 
may  be  requested  by  the  calling  program  by  setting  NLINES  = -1. 

ROUTINES  CALLED 
None. 

CALLED  BY 

HEADER  is  called  by  all  CSPM  routines  that  print  formatted  output. 

LOCAL  VARIABLES  ! 

NPAGE  - Last  page  number  printed.  j 

MAXL  - Maximum  number  of  lines  printed  per  page. 

NUM  - Number  of  lines  already  printed  on  the  current  page.  | 

SUBROUTINE  DESCRIPTION 

HEADER  can  be  used  in  two  ways.  Normally,  it  is  called  whenever  a 
routine  is  edsout  to  print  one  or  more  lines  that  should  appear  on  the  same 
page,  to  Insure  that  space  will  be  availcible.  However,  if  a table  is  being 
printed  that  extends  over  more  than  one  page,  it  is  more  efficient  for 
the  calling  routine  to  count  locally  the  number  of  lines  on  the  page  and  j 


} 
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only  call  HEADER  when  the  current  page  is  filled.  This  scheme  enables 
the  calling  routine  to  print  titles  for  the  columns  on  each  page.  HEADER 
is  called  with  NLINES  = -1  whenever  it  is  desirable  to  begin  printing  on 
a new  page. 


B.2.2  (1,0)  Overlay  Routines 

B . 2 . 2 . 1 INPUT 

Program  INPUT  is  the  controlling  routine  for  the  (1,0)  overlay.  It 
calls  NXTVAR  to  read  options  from  the  CSPM  input  data  cards,  checlcs  them 
for  errors,  stores  flags  and  data  values  in  common  bloc)c  /OPTIONS/  (see 
Table  B-4) , and  calls  PRINT  to  print  the  CSPM  input  data  summary.  It  also 
determines  if  a restart  file  is  available  and  if  so,  sets  the  proper  vari- 
ables. If  restart  has  not  been  selected,  INPUT  writes  the  OPTIONS  record 
(see  table  B-2)  to  file  LOSS. 

COMMON  BLOCKS 

INPUT  defines  the  common  bloc)c  /NXTVAR/  which  is  used  to  communicate 
with  the  routine  NXTVAR.  It  is  described  in  Table  B-5. 


ROUTINES  CALLED 

INPUT  calls  the  following  routines: 
DMS. 

HEADER. 

NXTVAR. 

PRINT. 

CALLED  BY 

INPUT  is  called  by  CSPM. 


LOCAL  VARIABLES 
KEYWORD 
NKEY 
KFLAG 


JTITLE 

JDATE 


Array  containing  list  of  legal  lieywords. 

Number  of  legal  Iceywords. 

Set  to  0 if  Iceyword  has  not  yet  been  read,  to  -1  if 
Iceyword  is  illegal  for  a restart  run,  to  1 if  Iceyword 
has  been  read. 

Title  of  restart  run. 

Date  of  restart  run. 


i 
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JTIME  Time  of  restart  run. 

IVAL  Array  equivalenced  to  integer  variables  in  /OPTIONS/, 

NAME  Array  equivalenced  to  Hollerith  variables  in  /OPTIONS/. 

VAL  Array  equivalenced  to  real  variables  in  /OPTIONS/. 

X Array  used  for  latitude-longitude  coordinates. 

GROUND  Array  containing  legal  values  for  option  GND. 

GNDDAT  Array  containing  values  of  SGM/EPS  for  GND  type. 

DRQ  Default  values  for  RQ. 

VLIM  Range  limits  for  real  variables  in  /OPTIONS/. 

KUNITS  Units  for  real  variables  in  /OPTIONS/. 

PROGRAM  DESCRIPTION 

INPUT  controls  the  reading  and  processing  of  user  input  options. 

It  first  initializes  necessary  flags  and  data  values  and  then  reads  the 
title  card  for  the  run.  INPUT  then  checks  for  a RESTRT  file  by  atten^ting 
to  read  the  OPTIONS  record  from  it.  If  the  read  is  successful,  the  /OPTIONS/ 
common  block  will  be  reset  to  the  values  used  when  the  RESTRT  file  was 
created.  INPUT  then  sets  the  restart  flag  and  identifies  the  options  that 
may  not  be  respecified  for  a restart  run  as  illegal.  If  an  EOF  is  read, 

INPUT  sets  the  flags  for  an  initial  run.  In  either  case,  INPUT  then  proceeds 
to  read  the  options  from  the  user  input  data  cards.  It  calls  NXTVAR  to 
read  data  cards  as  needed  and  return  the  next  syntactic  element  it  finds 
(keyword,  numeric  value,  or  delimiter) . As  each  new  option  keyword  is  read, 
it  is  checked  against  a list  of  legal  options.  INPUT  also  checks  to  see 
if  that  option  has  been  specified  previously  during  the  run,  if  the  option 
is  illegal  with  restart  (if  a restart  run  has  been  specified) , or  if  the 
option  conflicts  with  another  option  previously  read.  Once  this  preliminary 
checking  is  completed,  the  program  branches  to  the  section  of  INPUT  that 
reads  the  values  given  for  the  selected  option.  Separate  processing  sections 
are  used  for  the  following  groups  of  options: 

(1)  PATH,  PRO,  PTAB,  PLOT,  TABL,  LB,  PD,  RSL,  REL. 

(2)  MAP,  XMTRCOOR,  RCVRCOOR, 

(3)  XLOC,  RLOC,  POL,  CLIM,  GND. 

(4)  FREQ,  PWR,  HGl,  HG2,  Gl,  G2,  EPS,  SGM,  NS,  NO,  RSLTH. 

(5)  VI,  V2,  V3,  RQ. 
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NXTVAR  is  used  to  return  the  fields  that  correspond  to  the  values  of  the 
option.  The  syntax  of  the  option  is  checked  in  addition  to  the  value  given. 
Numeric  values  are  checked  to  be  within  specified  limits,  and  keywords  must 
be  one  of  a legal  list  of  options.  Once  all  the  option  keywords  have  been 
processed  from  the  user  input  data  cards,  INPUT  cross-checks  all  the  selections 
to  be  sure  that  no  information  is  missing  and  that  no  conflicting  options 
have  been  chosen.  INPUT  issues  both  warning  and  fatal  error  diagnostics. 
Warning  messages  are  given  when  execution  can  continue,  but  the  user  should 
be  made  aware  of  some  possibly  erroneous  input  data.  Fatal  error  diagnostics 
are  issued  when  the  input  data  is  incomplete,  in  conflict,  out  of  range, 
or  contains  syntax  errors.  Fatal  errors  will  cause  the  CSPM  run  to  be  aborted. 
The  PRINT  routine  then  is  called  to  print  the  input  data  summary  which  lists 
all  the  options  and  values  selected  by  the  user  together  with  the  default 
values  that  will  be  used  for  this  CSPM  run.  Finally,  if  the  run  is  not 
a restart  run,  the  OPTIONS  record  is  written  to  file  LOSS,  and  control  is 
returned  to  CSPM. 


B . 2 . 2 . 2 DMS 

Subroutine  DMS  converts  a latitude -longitude  coordinate  in  decimal 
degrees  north  and  east  to  a Hollerith  encoded  format  ready  to  print  in  positive 
degrees/minutes/seconds  and  decimal  degrees  with  the  proper  direction  North, 
South,  East,  and  West  appended. 

ARGUMENTS 

N Number  of  arguments  to  be  encoded  (2  or  4) . 

ZLAT  Latitude  in  decimal  degrees  north. 

ZLON  Longitude  in  decimal  degrees  east. 

YLAT  Second  latitude  if  map  boundary  is  input. 

YLON  Second  longitude  if  map  boundary  is  input. 

IDMS  Array  returned  with  positive  degrees/mlnutes/seconds  and 

proper  direction  for  latitude  cuid  longitude  values  input 
in  Hollerith  format. 

IDEC  Array  returned  with  positive  decimal  degrees  and  proper 

direction  for  latitude  euid  longitude  values  input  in 
Hollerith  format. 
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ROUTINES  CALLED 


None. 

CALLED  BY 

DMS  is  called  by  INPUT. 

PROGRAM  DESCRIPTION 

DMS  is  used  to  convert  the  latitude  and  longitude  values  for  the  transmit- 
ter and  receiver  coordinates  from  decimal  degrees  north  2md  east  to  an  encoded 
Hollerith  format  ready  for  printing.  The  latitude  and  longitude  boundaries 
for  a map  also  are  encoded  by  DMS  into  a print-ready  format.  DMS  initially 
takes  the  absolute  value  of  the  latitude  or  longitude  in  decimal  degrees 
and  determines  the  proper  direction  North,  South,  East,  or  West  to  be  used 
with  the  positive  value.  The  decimal  degrees  format  then  is  e '•coded  into 
array  IDEC.  DMS  then  converts  the  decimal  degrees  to  degrees/minutes/ seconds 
format  and  encodes  these  values  into  array  IDMS. 

B . 2 . 2 . 3 NXTVAR 

NXTVAR  is  used  to  read  the  free-format  user  input  data  cards.  It  reads 
successive  data  cards  until  an  EOF  is  found.  Each  card  is  printed  when 
it  is  read,  so  that  error  messages  pertaining  to  options  on  that  card  will 
immediately  follow  the  card  image.  On  each  call  to  NXTVAR,  the  next  keyword 
or  numeric  value  along  with  the  following  delimiter  is  returned.  If  the 
variable  read  is  not  of  the  type  requested,  an  error  message  is  printed 
and  the  error  flag  is  set.  All  80  columns  of  the  data  card  are  processed. 

Blanks  are  ignored.  The  five  characters  =,/()  and  the  end  of  the  card  are 
recognized  as  delimiters.  The  first  nonnumeric  character  following  a numeric 
field  also  is  considered  to  be  a delimiter.  Integer  values  are  returned 
as  real  numbers.  If  a syntax  error  is  found,  the  rest  of  the  data  card 
is  skipped.  Multiple  input  values  may  be  given  on  a single  data  card,  but 
a single  input  definition  may  not  sp2m  card  boundaries. 

ARGUMENTS 

ITYPE  The  type  of  variable  expected.  Set  to  one  for  keyword,  ! 

two  for  numeric  value,  and  three  for  delimiter.  If  the  | 

variable  field  read  does  not  match  the  type  expected  by  I 
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the  calling  routine,  an  error  message  is  printed,  and 
the  error  flag  is  set. 


r 


COMMON  BLOCKS 

Common  block  /NXT'/AR/  is  used  to  return  the  alphanumeric  keywords, 
numeric  values,  and  delimiters  read  as  well  as  an  error  flag.  The  variables 
in  this  common  block  are  described  in  Tedile  B-5. 

ROUTINES  CALLED 

NXTVAR  calls  the  routine  HEADER. 

CALLED  BY 

NXTVAR  is  called  by  INPUT. 

LOCAL  VARIABLES 

ICARD  Array  containing  the  input  data  card  being  processed  in 

80  Rl  format. 

ICOL  Pointer  set  to  the  position  of  the  last  character  used 

from  ICARD. 

JVAR  Characters  read  from  ICARD  are  packed  into  JVAR  until  a 

delimiter  is  found.  The  variable  is  encoded  later  into 
IVAR  if  it  is  a keyword,  or  decoded  into  XNUM  if  it  is 
determined  to  be  a numeric  value. 

The  following  flags  are  used  to  decide  the  type  of  the  variable  read 
and  to  diagnose  errors: 

JALPHA  Nonzero  if  a nonnumeric  character  is  packed  into  JVAR. 

JDECPT  Number  of  decimal  points  in  JVAR. 

JPLUS  Number  of  plus  signs  in  JVAR. 

JMINUS  Number  of  minus  signs  in  JVAR. 

JFIRST  First  nonblank  character  read. 

KNT  Number  of  characters  read  before  a delimiter  is  found. 

PRCXiRAM  DESCRIPTION 

NXTVAR  is  called  to  return  one  syntactical  element  at  a time  except 
in  the  case  of  an  alphanumeric  or  numeric  field  in  which  case  the  following 
delimiter  also  is  returned  to  the  calling  routine.  This  enables  the  INPUT 
progreim  to  look  only  at  the  synt<uc  of  the  options  and  their  values  without 
being  concerned  about  skipping  blanks,  reading  the  next  card,  or  determining 
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if  characters  in  a field  are  a keyword,  a real  number,  an  integer,  or  a 
delimiter.  NXTVAR  reads  a data  card  into  the  array  ICARD  in  80  Rl  format 
when  all  characters  on  the  previous  card  have  been  processed.  It  looks 
for  the  next  syntactic  element  by  taking  successive  characters  from  the 
last  column  position  processed  on  the  last  card  toward  the  end  of  the  card 
until  a delimiter  is  found  or  more  than  10  characters  are  processed.  These 
characters  are  packed  into  variable  JVAR.  While  each  character  from  ICARD 
is  being  processed,  certain  information  is  kept  about  the  type  of  characters 
found.  Blanks  are  ignored  and  are  not  put  into  JVAR.  If  a minus,  plus, 
or  decimal  point  is  found,  the  value  of  the  field  must  be  numeric.  Only 
one  of  each  of  these  characters  may  appear  in  one  field,  or  an  error  will 
be  diagnosed.  If  the  first  nonblank  character  in  a field  is  alpha,  then 
the  field  must  be  a keyword.  Keywords  are  returned  in  AlO  format  in  the 
variable  IVAR.  Niuneric  fields  are  decoded  into  a real  numeric  value  and 
returned  in  the  varicd>le  XNUM.  Delimiter  characters  that  are  found  are 
returned  in  Rl  format  in  the  variable  IDEL.  The  error  flag  lERR  will  be 
zero  if  no  errors  are  diagnosed,  and  one  if  an  error  is  found. 

B.2.2.4  print 

Subroutine  PRINT  prints  the  CSPM  input  data  summary.  All  the  values 
that  it  prints  have  been  stored  previously  in  coninon  block  /OPTIONS/.  They 
are  described  in  Table  B-4. 

ARGUMENTS 

None. 

ROUTINES  CALLED 

PRINT  calls  the  routine  HEADER. 

CALLED  BY 

PRINT  is  called  by  the  INPUT  routine. 

SUBROUTINE  DESCRIPTION 

PRINT  is  called  to  print  the  CSPM  input  data  summary  after  all  the 
user  input  options  have  been  read  and  processed.  The  data  values  and  option 
flags  selected  by  the  user  together  with  the  default  values  set  by  the  CSPM 
are  printed  in  the  input  data  summary,  so  that  the  user  will  have  a complete 
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record  of  all  values  used  in  this  CSPM  run.  All  the  necessary  data  for 
PRINT  to  print  are  stored  in  the  /OPTIONS/  common  block.  The  values  that 
are  printed  depend  on  the  type  of  output  selections  made.  The  format  also 
varies  between  runs  for  a single  path  and  those  for  analyzing  a geographic 
area . 

B.2.3  (2,0)  Overlay  Routines 

B.2.3.1  LOSSMOD 

Program  LOSSMOD  is  the  controlling  routine  for  the  (2,0)  overlay. 

It  controls  the  process  of  extracting  path  profiles  from  the  topographic 
data  base  and  calculating  ACR  and  statistical  values  with  the  Longley-Rice 
propagation  prediction  routines  for  a single  path  or  points  throughout 
a geographic  area.  If  plotted  or  tabulated  path  profiles  have  been  requested, 
the  path  elevations  are  written  to  file  LOSS  to  be  used  by  the  OUTPUT  (3,0) 
overlay.  If  basic  transmission  loss  has  been  requested  for  a single  path, 
or  basic  transmission  loss,  power  density,  received  signal  level,  or  communi- 
cation reliability  data  have  been  requested  for  contour  plots  or  tabulated 
output,  the  ACR  DATA  records  for  the  necessary  paths  also  are  written  to 
file  LOSS. 

COMMON  BLOCKS 

LOSSMOD  uses  many  common  blocks  to  communicate  values  with  other  routines 
It  references  the  input  data  values  and  option  flags  stored  in  the  /OPTIONS/ 
comnon  block  and  generates  the  values  to  be  passed  to  the  OUTPUT  overlay 
in  the  /LOSSMOD/  common  block.  These  are  described  in  Tables  b-4  and  B-5 
The  variables  in  common  blocks  /PPILE/  and  /GEODYS/  described  in  Tables 
B-7  and  B-8  are  used  to  communicate  with  the  path  profile  extraction  routines. 
Common  blocks  /EVPARS/,  /PHPARS/,  and  /SYPARS/  contain  variables  used  to 
describe  characteristics  of  the  system  and  are  used  by  the  routines  that 
calculate  path  parameter  information.  These  common  blocks  are  described 
in  Tables  B-9,  B-10,  and  B-11.  Notice  that  some  of  the  same  values  that 
appear  in  these  common  blocks  are  repeated  in  the  common  blocks  used  with 
the  liOngley-Rice  propagation  loss  routines,  /LRINP/,  /LRDAS/,  LRVAR/,  and 
LRRTH/  which  are  descrii)ed  in  Tables  B-12,  B-13,  B-14,  and  B-15.  Because 
the  Longley-Rice  propagation  loss  prediction  routines  that  are  used  in 
the  CSPM  may  be  replaced  at  a future  time,  it  is  considered  necessary  to 


keep  their  common  blocks  intact.  Consequently,  the  path  pa^cuneter  and 
system  information  contained  in  common  blocks  /OPTIONS/,  /EVPARS/,  /PHPARS/, 
and  /SYPARS/  simply  is  transferred  to  the  proper  position  in  the  Longley- 
Rice  model  common  blocks  before  those  routines  are  called. 

ROUTINES  CALLED 

LOSSMOD  calls  the  following  routines: 

PFLTPO 

DESIC 

QLRPP  (entry  point  in  QLRAP) 

QPARS  (and  the  alternate  entry  point  QPARX) 

ARDESIC 

QRTHETA 

QVARLR  (and  the  alternate  entry  point  QQVARLR) . 

CALLED  BY 

LOSSMOD  is  called  by  CSPM. 

PROGRAM  DESCRIPTION 

LOSSMOD  is  used  to  extract  path  profiles,  calculate  path  parameter 
information,  and  calculate  the  ACR  and  statistical  data  for  either  a single 
path  or  for  points  throughout  a geographic  area.  This  information  is  written 
to  file  LOSS  and  later  is  used  by  the  (3,0)  overlay,  OUTPUT,  to  produce 
tabulated  and  plotted  path  profiles  and  tcibulate  output  or  contour  plots 
for  basic  transmission  loss,  power  density,  received  signal  level,  and 
coinnunication  reliability. 

If  a single  path  has  been  specified,  LOSSMOD  first  calls  PFLTPO  to 
extract  the  path  profile  elevations  from  the  topographic  data  base.  Next, 
QPARS  is  called  to  calculate  the  path  parcuneters  from  the  input  data  in 
/OPTIONS/  and  the  path  profile  elevations  in  /PFILE/.  Next,  LOSSMOD  sets 
a flag  to  tell  the  OUTPUT  routine  whether  the  path  profile  is  to  be  plotted, 
tabulated,  or  plotted  and  t^UDulated.  Then  the  PROFILE  record  containing 
the  output  type  flag,  the  path  parameter  values,  and  the  profile  elevations 
is  written  to  file  LOSS.  See  Table  B-2  for  the  structure  of  this  record. 
LOSSMOD  then  determines  if  basic  transmission  loss  has  been  requested  for 
the  path.  If  so,  it  transfers  the  necessary  data  values  from  the  path 
parameter  common  blocks  to  those  used  by  the  Longley-Rice  propagation  loss 


routines.  The  climate  dependent  variables  are  initialized  by  calling  QVARLR 
with  the  climate  code  specified  in  the  /OPTIONS/  common  block.  Next,  QLRPP 
is  called  to  initialize  necessary  variables  used  in  the  propagation  loss 
routines  for  the  point-to-point  mode  and  to  calculate  the  ACR  coefficients 
for  the  diffraction  region.  IXDSSMOD  then  calls  QQVARLR  with  the  total 
path  length  to  calculate  the  ACR  and  statistical  values  for  the  path. 

These  values  are  returned  in  the  common  block  /LRVAR/.  This  information 
is  written  to  file  LOSS  as  an  ACR  DATA  record  in  the  format  described  in 
Table  B-2.  Control  then  is  returned  to  CSPM. 

The  procedure  followed  for  calculating  values  for  a geographic  area 
repeats  the  process  used  for  a single  path  for  each  data  point  in  the  area. 
LOSSMOD  first  calls  QVARLR  to  initialize  the  climate  dependent  variables 
for  the  area  being  studied.  The  data  points  to  be  used  throughout  the 
geographic  area  are  selected  by  the  method  described  in  Appendix  F.  Basically 
the  method  involves  determining  how  many  radials  are  required  to  adequately 
cover  the  area  and  then  using  the  system  described  in  Table  F-1  to  locate 
data  points  along  each  radial.  To  determine  the  number  of  radials  to  be 
used,  LOSSMOD  first  calculates  the  maximum  distance  from  the  transmitter 
to  any  point  on  the  area  boundary.  The  number  of  radials  to  be  used  is 
then  equal  to  the  largest  power  of  two  that  is  less  than  this  maximum  distance 
Once  the  number  of  radials  has  been  determined,  the  angular  increment  between 
radials  is  calculated.  If  the  transmitter  is  located  on  the  area  boundary, 
the  initial  and  final  angular  values  are  restricted,  so  that  only  those 
radials  that  are  interior  to  the  area  are  used. 

LOSSMOD  loops  through  all  radials  used  for  the  area.  A radial  is 
skipped  if  its  data  values  are  not  to  be  tabulated  and  no  contour  plots 
have  been  requested.  If  any  contour  plots  have  been  selected,  data  are 
generated  for  all  radials  so  that  the  spacing  of  data  points  will  be  uniform 
across  the  entire  area.  ACR  and  statistical  values  arc  generated  for  each 
data  point  along  the  radial.  This  is  done  by  considering  each  radial  to 
be  like  a series  of  single  paths  from  the  transmitter  to  each  data  point 
on  the  radial.  The  distance  between  the  transmitter  and  the  first  data 
point  to  be  used  on  that  radial  is  determined  according  to  the  method  shown 
in  Table  F-1.  To  Insure  that  the  contour  plots  will  be  accurate  close 

to  the  area  boundary,  one  extra  data  point  outside  the  boundary  is  used 
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for  each  radial.  The  location  of  this  point  is  determined  by  using  QRTHETA 
to  find  the  point  where  the  radial  intersects  the  boundary,  calculating 
the  distance  from  the  transmitter  to  that  point,  rounding  that  distance 
up  to  the  nearest  kilometer,  and  calling  QRDESIC  to  calculate  the  coordinates 
of  the  point  that  lies  on  the  radial  at  that  distance  from  the  transmitter. 
Data  points  are  distributed  then  at  one  kilometer  intervals  along  the  radial 
from  the  starting  point  to  the  final  point  outside  the  boundary.  For  any 
radial,  if  the  distance  to  the  starting  data  point  is  greater  than  the 
distance  to  the  boundary,  the  radial  is  skipped. 

PFLTPO  is  called  to  extract  the  path  profile  from  the  topographic 
data  base  for  the  path  from  the  transmitter  to  the  final  data  point  outside 
the  boundary.  The  point  spacing  on  the  profile  is  fixed  at  250  m,  since 
the  path  length  is  always  an  even  number  of  kilometers.  Thus,  the  profile 
can  be  used  in  sections  (four  points  to  a kilometer)  to  provide  the  profile 
data  for  any  one  of  the  subpaths  from  the  transmitter  to  a data  point. 

QPARS  is  called  to  calculate  the  path  parameters  for  the  entire  radial 
path.  If  the  profile  for  this  radial  has  been  selected  to  be  tctbulated 
or  plotted,  a PROFILE  record  is  written  to  file  LOSS. 

LOSSMOD  then  loops  through  all  the  subpaths  for  this  radial  calling 
QPARX  to  calculate  the  path  parameters  for  the  profile  values  between  the 
transmitter  and  each  data  point.  QPARX  is  an  entry  point  in  QPARS  that 
approximates  AH  for  the  subpath  rather  than  calculating  the  actual  inter- 
decile range  of  the  profile  elevations.  Next,  QLRPP  is  called  to  initialize 
the  diffraction  range  coefficients  for  the  path  and  to  set  values  for  the 
point-to-point  mode  propagation  loss  routines.  Finally,  QQVARLR  is  called 
with  the  distance  from  the  transmitter  to  the  data  point  (the  length  of 
the  subpath)  to  calculate  the  ACR  and  statistical  values  at  the  data  point. 
These  values  are  stored  in  the  DATA  array.  When  all  the  data  points  for 
this  radial  have  been  processed,  the  ACR  DATA  record  for  the  radial  is 
written  to  file  LOSS.  This  process  continues  until  the  data  values  have 
been  calculated  for  all  the  radials  throughout  the  area.  The  LOSS  file 
then  contains  all  necessary  data  to  plot  or  tcd>ulate  basic  transmission 
loss,  power  density,  received  signal  level,  or  communication  reliability 
for  the  geographic  area,  and  control  is  returned  to  CSPM. 
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B . 2 . 3 . 2 ELVTPO 

Function  ELVTPO  returns  the  elevation  in  meters  for  a point  within 
range  of  the  topographic  data  base.  The  elevation  is  determined  by  using 
bivariate  linear  interpolation  between  the  four  surrounding  data  points 
on  the  topographic  data  base. 

ARGUMENTS 

ALAT  Latitude  of  the  point  in  decimal  degrees  north. 

ALON  Longitude  of  the  point  in  decimal  degrees  east. 

lERR  Error  flag: 

(0)  no  errors. 

i 

(1)  parity  error  in  index  block. 

(2)  parity  error  in  data  block. 

(8)  missing  data  was  filled  in  with  an  averaging 
technique . 

(-1)  point  is  out  of  range  of  the  topographic  data 
base . 

(-2)  data  block  is  missing. 

(-3)  too  many  missing  data  points. 

' The  principle  here  is  that  errors  with  negative  codes  are  serious  and  the 

point  is  definitely  unattainable,  while  positive  error  codes  indicate  th^t 
the  data  were  obtained  under  suspicious  circumstances.  In  the  case  of 
I a positive  error  code  an  attempt  is  made  to  return  a correct  value,  but 

I for  a negative  code,  there  is  no  such  attempt  and  the  returned  value  is 

^ meaningless. 

Furthermore,  the  routine  is  written  in  such  a way  that  in  a series 
of  calls  (as  for  example  along  a great  circle  path) , the  error  return  need 
j be  checked  only  at  the  end.  The  normal  return  leaves  lERR  unaffected  which 

is  why  it  should  be  set  initially  to  0.  A positive  error  code  is  logically 

f 

* or'ed  with  lERR,  and  a negative  code  turns  off  all  other  attempts  at  making 

evaluations . 

COMMON  BLOCKS 

/PKDTPO/  This  common  block  is  defined  in  subroutine  LDDTPO. 

ROUTINES  CALLED 


ELVTPO  calls  LDDTPO. 
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CALLED  BY 


ELVTPO  is  called  by  PFLTPO. 

FUNCTION  DESCRIPTION 

Upon  entry,  the  routine  first  assures  that  the  proper  1°  by  1°  block 
is  loaded  in  /PKDTPO/  using,  if  necessary,  a call  to  LDDTPO.  It  then  unpacks 
the  elevations  at  the  four  corners  of  the  30"  by  30"  square  within  which 
the  given  point  falls.  The  returned  value  is  obtained  using  the  standard 
bilinear  interpolation  between  these  four  points. 

If  from  one  to  three  of  the  corner  values  are  missing  from  the  tape, 
they  are  filled  with  the  simple  average  of  the  remaining  corner  values, 
and  bit  number  3 of  lERR  is  set.  If  all  four  values  are  missing,  then 
lERR  is  replaced  by  -3.  If  upon  first  entering  the  value  of  lERR  is  negative, 
the  entire  routine  is  skipped. 

B . 2 . 3 . 3 FFGND 

Subroutine  FFGND  calculates  an  effective  surface  elevation  for  the 
terrain  between  the  transmitting  antenna  and  its  horizon  and  for  the  terrain 
between  the  receiving  antenna  and  its  horizon  for  paths  which  are  not  line- 
of-sight.  For  a path  that  is  determined  to  be  line-of-sight,  this  subroutine 
calculates  effective  surface  elevation  for  the  terrain  between  the  trans- 
mitting and  receiving  antennas.  Effective  surface  elevation  is  defined 
as  the  average  elevation  for  the  middle  80%  of  the  profile  points  comprising 
the  path  or  path  segment  (as  in  the  case  of  non  line-of-sight  paths) . 

ARGUMENTS 
None . 

COMMON  BLOCKS 

FFGND  uses  common  blocks  /PFILE/,  defined  in  table  B-7  and  /PHPARS/, 
defined  in  Table  B-10. 

ROUTINES  CALLED 
None . 

CALLED  BY 

FFGND  is  called  by  QPARS . 
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SUBROUTINE  DESCRIPTION 

FFGND  calculates  an  effective  elevation  for  the  foreground  of  an  antenna 
under  certain  conditions  for  which  use  of  actual  foreground  elevation  would 
result  in  propagation  loss  predictions  which  would  be  quite  different  from 
measured  loss  over  that  path.  The  calculation  of  effective  foreground 
elevation  is  accomplished  by  first  defining  four  variables  which  express 
the  number  of  profile  points  from  the  transmitter  to  other  locations  along 
the  path.  The  first  variable,  NL(1)e  0,  is  the  number  of  points  to  the 
transmitter  location.  The  second  variable,  NL(2) , is  the  number  of  points 
to  the  horizon  for  the  transmitting  antenna.  The  third  variable,  NL(3), 
is  the  number  of  points  (from  the  transmitter)  to  the  horizon  for  the  receiv- 
ing antenna.  The  fourth  variable,  NL(4),  is  the  number  of  points  (from 
the  transmitter)  to  the  receiver  location.  This  value  is,  in  fact,  the 
number  of  points  in  the  profile  array. 

These  four  variables  are  used  to  define  segments  of  the  path.  The 
path  segments  of  interest  for  a transhorizon  path  are  the  segments  from 
the  transmitting  antenna  to  its  horizon  and  the  segment  from  the  receiving 
antenna  to  its  horizon.  For  a line-of-sight  path,  the  entire  path  is  consid- 
ered. For  each  segment  and  the  entire  path,  the  middle  80%  of  the  segment/ 
path  elevations  are  used  to  calculate  an  average  elevation.  A check  is 
made  on  the  value  for  variable  ISO  provided  in  common  block  /PHPARS/. 

When  ISOO,  the  path  is  line-of-sight.  A value  of  1 or  -1  denotes  a transho- 
rizon path.  The  effective  surface  elevation (s)  then  is (are)  returned  to 
QPARS  using  variable  HSE. 


B.2.3.4  ISEARCH 

Function  ISEARCH  returns  the  value  of  the  smallest  index  at  which 
the  corresponding  element  in  the  array  JA  matches  the  word  K.  If  no  match 
is  found,  a zero  is  returned. 

ARGUMENTS 

K Contains  the  work  for  which  a match  in  JA  is  being  searched. 

N Number  of  words  to  be  searched  in  array  JA. 

JA  Array  of  words  to  be  searched. 
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TTTOSS 


ROUTINES  CALLED 
None . 

CALLED  BY 

ISEARCH  is  called  by  LDDTPO, 

FUNCTION  DESCRIPTION 

ISEARCH  is  an  integer  function.  It  loops  through  the  first  N elements 
of  array  JA  looking  for  a match  with  K.  If  a match  is  found,  ISEARCH  returns 
the  index  value  in  array  JA  at  which  the  matching  value  is  stored.  If 
no  match  is  found,  a zero  is  returned.  JA  and  K may  have  either  integer 
or  Hollerith  values. 


B . 2 . 3 . 5 LDDTPO 

LDDTPO  is  the  basic  subroutine  for  reading  the  topographic  data  tapes. 
Given  the  latitude  and  longitude  of  the  southwest  corner  of  a data  block, 
the  routine  finds  the  proper  record  and  reads  it  into  the  common  block 
/PKDTPO/.  The  data  are  still  in  packed  form,  but  enough  information  is 
supplied  to  make  extraction  of  the  desired  values  reasonably  easy.  Most 
jobs  require  shuttling  back  and  forth  between  two  or  more  records.  To 
increase  efficiency  and  to  reduce  wear  and  tear  on  the  topographic  data 
tape,  a mechanism  is  built  in  for  saving  and  retrieving  such  records. 

A scratch  file  on  disc  is  used  in  which  up  to  ten  records  are  stored. 

Normally,  LDDTPO  is  not  called  directly  by  the  user,  but  only  indirectly 
through  routine  ELVTPO.  This  routine  not  only  makes  the  proper  calls  to 
LDDTPO,  but  it  also  unpacks  the  data  into  immediately  usable  forms. 

LDDTPO  uses  local  file  names  TOPO,  for  the  magnetic  tape,  and  TOPS, 
for  the  scratch  file.  File  TOPS  is  a word-addressable  file  which  is  handled 
automatically,  hence  it  will  not  be  discussed  further.  The  statement 

CALL  LDDTPO (LAT, LON) 

causes  one  particular  record  to  be  read.  The  variables  LAT,  LON  are  the 
latitude  and  longitude  of  the  southwest  corner  of  the  desired  data  block. 

They  are  integers  and  give  degrees  measured  north  of  the  equator  and  east 
of  the  Greenwich  meridian.  In  CONUS  the  latitude  is  positive  and  the  longitude 
negative.  The  output  is  contained  in  the  common  block  /PKDTPO/  which  is 
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defined  later  in  this  section. 


FILES 

The  files  used  by  LDDTPO  have  been  identified  as  TOPO  and  TOPS.  File 
TOPO,  containing  the  topographic  data  on  magnetic  tape,  is  described  in 
Appendix  E.  As  mentioned,  TOPS  is  a word-addressable,  scratch  file  to 
which  blocks  of  data  read  from  TOPO  are  written  automatically. 

ARGUMENTS 


LAT 

Latitude  in  decimal  degrees  north  of  the  southwest  comer 

of  the  requested  data  block. 

LON 

Longitude  in  decimal  degrees  east  of  the  southwest  comer 

of  the  requested  data  block. 

CX)MMON  BLOCKS 

/PREFACE/ 

Contains  information  that  was  read  from  the  tape  header 

record  by  READPRF  for  the  topographic  data  tape  being  used 

(file  TOPO) . 

NNFO 

Number  of  words  of  information  contained  in  the  topographic 

data  tape  header  record. 

FMT 

Eight-word  array  containing  the  FTN  FORMAT  statement  used 

by  PRNTPRF  to  print  the  information  contained  in  the  header 

record . 

TITLE 

Four-word  array  containing  the  title  of  the  topographic 

data  tape. 

LATP 

Two-word  array  containing  latitude  range  bounds  for  the 

topographic  data  tape. 

LONP 

Two-word  array  containing  longitude  range  bounds  for  the 

topographic  data  tape. 

KODE 

Integer  giving  the  key  to  format  of  the  tape,  type  of  data, 

etc.  Must  be  •>!  to  be  compatible  with  the  CSPM. 

JNFO 

Eleven-word,  unused  array. 

/PKDTPO/ 

Contains  data  read  from  the  topographic  data  tape  for  a 

1°  by  1°  block. 

PKDN 

A floating  point  number  defined  as  1.000001/n,  where  n is 

the  number  of  data  points  packed  into  a word. 
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A floating  point  number  defined  as  nb,  where  b is  the 
number  of  bits  used  for  a datum  (2,3,4,5,6,7,8,10,  or  12). 
KXR  A coded  word  which  indicates  parity  errors  realized  in 

reading  a record.  Values  are: 

KXR  = 0,  No  parity  errors. 

1,  A parity  error  while  reading  the  index  record 
of  the  latitude  data  partition. 

2,  A parity  error  while  reading  the  block  data 
record. 

3,  Parity  error  in  both  records. 

MSKD  A mask  of  b I's,  right  justified. 

LPKD  The  number  b of  bits  per  datum. 

NPKD  Status  indicator  for  the  block  data  record.  Values  are: 

NPKD  = -2  The  block  is  missing  from  the  tape. 

= -1  The  block  is  outside  the  range  of  the  tape. 

= 0 The  block  is  ocean. 

= n>0  Data  for  the  block  are  on  the  tape  and 

packed  n elevations  per  word. 

MPKD  The  block  data  record. 

MPKD(l)  gives  the  latitude  and  longitude  of  the  block. 

MPKD (2)  gives  the  packing  information  required  to  decode 
the  data. 

MPKD{3)  gives  the  bias  to  be  applied  to  the  data. 

MPKD(7)  is  the  first  word  of  the  array  of  elevations. 
ROUTINES  CALLED 

LDDTPO  calls  subroutine  READPRF,  ISEARCH,  and  various  Record  Manager 
routines . 

CALLED  BY 

LDDTPO  is  called  by  function  ELVTPO. 

LOCAL  VARIABLES 

TOPO  - Local  name  for  the  file  containing  the  topographic  data  on 
magnetic  tape. 

TOPS  - Local  name  for  the  word-addressable,  scratch  file  to 
which  block  of  data  are  written  from  TOPO. 
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LATX 

LONX 

KODE 


WPRF 

LDRS 


JDRS 

IX 


The  latitude  boundaries  for  the  topographic  data  tape. 

The  longitude  boundaries  for  the  topographic  data  tape . 

A check  for  proper  topographic  data  tape.  In  the 
topographic  data  tape  preface,  KODE  = 1 indicates  format 
of  file  TOPO  is  compatible  with  LDDTPO. 

A logical  flag  used  to  determine  if  fits  for  TOPO  and  TOPS 
have  been  created  and  the  necessary  fields  intialized. 

The  index  for  file  TOPS,  Index  consists  of  latitude  and 
longitude  (of  southwest  corner)  for  blocks  that  have  been 
written  to  file  TOPS. 

An  integer  array  containing  usage  information  for  the 
blocks  of  data  contained  in  file  TOPS. 

An  integer  array  used  to  contain  the  index  records  for  each 
latitude  partition  read  from  file  TOPO. 


SUBROUTINE  DESCRIPTION 

Upon  entry  the  routine  first  checks  in  sequence  whether  the  desired 
block  is  already  in  /PKDTPO/,  whether  the  tape  reading  has  been  Initialized, 
or  whether  the  desired  block  is  stored  in  TOPS. 

In  the  initialization  process,  the  preface  record  is  read  from  the 
tape,  the  value  of  KODE  is  checked,  and  the  extremes  of  latitude  and  longitude 
are  noted.  The  scratch  file  TOPS  is  defined  to  be  a Word  Addressable  (WA) 
file  with  a fixed  length  record  of  24530  characters — the  length  of  the 
common  block  /PKDTPO/.  In  addition,  a counter  is  set  to  indicate  to  the 

routine  what  latitude  data  partition  currently  is  being  read. 

If  the  record  is  in  /PKDTPO/,  the  routine  immediately  returns.  If 
it  is  in  TOPS,  it  is  transferred  to  /PKDTPO/  and  the  routine  returns. 
Otherwise,  the  desired  block  is  checked  to  see  if  it  lies  within  the  range 

of  the  tape,  and  only  then  does  the  routine  begin  to  read  the  tape.  The 

latitude  data  partition  counter  is  examined,  and  if  it  does  not  correspond 
to  the  desired  latitude,  the  proper  number  of  partitions  are  skipped,  the 
index  record  is  read,  and  a second  counter  is  set  to  indicate  what  record 
of  the  partition  can  be  read  next.  The  index  is  examined  at  the  desired 
longitude  and,  depending  on  its  value,  the  proper  action  is  taken.  If 
the  index  is  positive  the  routine  skips  the  proper  number  of  records  and 
finally  reads  the  desired  one  into  the  array  MPKD.  It  computes  the  unpacking 
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aids  using  MPKD(2) . 

Each  time  a new  data  block  is  estedDlished  in  /PKDTPO/,  no  matter  what 
its  status  NPKD,  the  routine  will  write  the  entire  common  block  into  TOPS. 

The  first  ten  such  records  are  written  in  sequential  order.  After  that, 
one  of  the  previously  written  records  is  deleted  from  the  stack  and  replaced 
by  the  present  one.  The  one  deleted  is  chosen  to  be  that  record  which 
has  not  been  used  (or  reused)  for  the  longest  period  of  time. 

Only  one  TOPO  tape  may  be  used  in  any  one  job.  Furthermore,  the  user 
should  not  issue  any  outside  commands  that  will  affect  the  positioning 
of  the  file  TOPO.  If  he  does,  the  synchronization  within  LDDTPO  may  be 
lost. 

At  present  the  parity  checks  are  not  implemented.  The  user  must  rely 
on  such  information  as  he  can  obtain  from  the  dayfile. 

B.2.3.6  PFLTPO 

Subroutine  PFLTPO  extracts  a path  profile  from  the  TOPO  data  base. 
Originally  written  by  George  Hufford,  this  version  has  been  modified  for 
use  with  the  CSPM.  The  maximum  number  of  points  per  profile  has  been  increased 
from  601  to  1251  to  achieve  higher  resolution  for  long  paths. 

ARGUMENTS 

ALAT  - Transmitter  latitude  in  decimal  degrees  north. 

ALON  - Transmitter  longitude  in  decimal  degrees  east. 

BLAT  - Receiver  latitude  in  decimal  degrees  north. 

BLON  - Receiver  longitude  in  decimal  degrees  east. 

XXI  - Distcmce  increment  between  points  of  the  profile  in  meters 
to  be  used  as  an  initial  estimate. 
lERR  - Error  flag. 

COMMON  BLOCKS 

PFLTPO  uses  the  common  blocks  /OPTIONS/  and  /PFILE/  described  in  Tables 
B-4  and  B-7. 

ROUTINES  CALLED 

PFLTPO  calls  the  following  routines: 
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DESIC 


ELVTPO 

QDESIC 
CALLED  BY 

PFLTPO  is  called  by  LOSSMOD. 

SUBROUTINE  DESCRIPTION 

PFLTPO  determines  the  number  of  profile  points  to  be  used  from  the 
path  length  and  the  initial  distance  increment  XXI . The  path  length  in 
meters  and  the  bearing  from  the  transmitter  to  receiver  is  calculated  by 
DESIC.  If  this  profile  is  for  a map  radial,  this  distance  is  rounded  to 
the  next  kilometer.  The  number  of  intermediate  points  on  the  path  at  which 
the  elevation  is  to  be  extracted  from  the  TOPO  base  is  then  calculated 
by  dividing  the  path  length  by  XXI  and  rounding  up  to  the  next  integer 
value.  A maximum  of  1251  profile  points  is  allowed.  The  actual  spacing 
between  points  is  determined  by  dividing  the  path  length  by  the  number 
of  intervals  between  profile  points.  PFLTPO  then  loops  through  all  the 
profile  points,  calling  QDESIC  to  return  the  latitude  and  longitude  of 
the  intermediate  point  and  ELVTPO  to  extract  the  elevation  from  the  TOPO 
base.  These  elevations  are  stored  sequentially  in  array  PFL,  with  the 
transmitter  site  elevation  stored  in  PFL(l)  and  the  receiver  site  elevation 
in  PFL{NP+1) . 


B.2.3.7  QDESIC 

Subroutine  QDESIC  calculates  the  latitqde  and  longitude  for  the  end- 
point of  a path  given  the  latitude  and  longitude  of  the  initial  point, 
the  bearing,  and  the  path  length.  The  alternate  entry  point,  DESIC,  uses 
the  latitude  and  longitude  of  both  endpoints  of  a path  to  calculate  the 
bearing  and  path  length. 

ARGUMENTS 

ALAT  - Latitude  of  point  1 in  decimal  degrees  north. 

ALON  - Longitude  of  point  1 in  decimal  degrees  east. 

BLAT  - Latitude  of  point  2 in  decimal  degrees  north. 

BLON  - Longitude  of  point  2 in  decimal  degrees  east. 

- Bearing  of  path  in  decimal  degrees  clockwise  from  north. 

121 


AZIM 


DST  - Path  length  in  meters. 


1 


ROUTINES  CALLED 
None. 

CALLED  BY 

QDESIC  and  DESIC  are  called  by  LDDTPO  and  PFLTPO. 

ENTRY  POINTS 

DESIC  is  called  to  calculate  path  length  and  bearing  given  the 
latitude  and  longitude  of  both  endpoints  of  the  path. 

SUBROUTINE  DESCRIPTION 

QDESIC  and  DESIC  use  the  relationships  of  spherical  trigonometry  to 
calculate  the  requested  values. 


I 

! B.2.3.8  QDLTH 


1 Subroutine  QDLTH  computes  DLTK,  the  terrain  irregularity,  for  a given 

' path  profile.  Terrain  irregularity  is  defined  as  the  interdecile  range 

j of  the  deviations  of  the  elevations  at  points  along  the  profile  from  the 

i corresponding  values  on  a least  squares  linear  approximation  to  the  profile. 

i 

I . ARGUMENTS 

I None . 

I 

COMMON  BLOCKS 

; QDLTH  uses  the  common  blocks  /PFILE/  and  /EVPARS/  described  in 

' Tables  B-7  and  3-9. 

ROUTINES  CALLED 

QDLTH  calls  QTILE. 

CALLED  BY 

QDLTH  is  called  by  QPARS. 

SUBROUTINE  DESCRIPTION 

QDLTH  initially  calculates  the  coefficient  for  a least-squares  linear 
approximation  to  the  profile  elevations.  It  then  calculates  the  deviation 
of  each  point  on  the  profile  from  the  least-squares  approximation.  QDLTO 
^ calls  QTILE  to  partially  order  this  array  of  deviation  values.  The 
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interdecile  range  then  is  defined  as  the  difference  between  the  10th  and 
90th  percentiles  of  the  array  of  deviations.  DLTH  is  an  asymptotic  value 
determined  as  a function  of  the  distance  and  the  interdecile  range. 

B.2.3.9  QPARS 

Subroutine  QPARS  calculates  the  path  parameters  used  by  the  Longley- 
Rice  routines  in  the  point-to-point  mode. 

ARGUMENTS 
None . 

COMMON  BLOCKS 

QPARS  uses  the  common  blocks  /PFILE/,  /EVPARS/,  /PHPARS/,  and  /SVPARS/ 
described  in  Tables  B-7,  B-9,  B-10,  and  B-11. 

ROUTINES  CALLED 

QPARS  calls  QDLTH  and  FFGND. 

CALLED  BY 

QPARS  is  called  by  LOSSMOD. 

ENTRY  POINTS 

QPARX  is  called  when  an  approximation  to  DLTH  is  to  be  calculated 
instead  of  calling  QDLTH. 

SUBROUTINE  DESCRIPTION 

QPARS  calls  QDLTH  to  calculate  the  terrain  irregularity,  DLTH.  HSYS 
is  the  average  of  the  transmitter  and  receiver  site  elevations.  ENS  may 
be  input  by  the  user  or  calculated  from  ENO  and  HSYS.  QPARS  also  calculates 
the  horizon  elevation  anagles,  distances  to  the  horizon,  elevations  at 
the  horizons,  and  effective  antenna  heights.  The  flag  ISO  is  set  to  0 
for  a line-of-sight  path  and  ft  0 otherwise. 


B . 2 . 3 . 10  QRTHETA 

^ Subroutine  QRTHETA  calculates  the  latitude  and  longitude  of  the  point 

I where  a radial  from  the  transmitter  at  bearing  AZIM  intersects  the  map 

i boundary . 

i 

I 
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ARGUMENTS 

AZIM  - Bearing  of  the  radial  in  decimal  degrees. 

ZLAT  - Latitude  of  the  intersection  point  in  decimal  degrees  north. 

ZLQN  - Longitude  of  the  intersection  point  in  decimal  degrees  east. 

COMMON  BLOCKS 

QRTHETA  uses  the  common  blocks  /OPTIONS/  and  /LOSSMOD/  described  in 
Tables  B-4  and  B-5  respectively. 

ROUTINES  CALLED 
None. 

CALLED  BY 

QRTHETA  is  called  by  LOSSMOD. 

SUBROUTINE  DESCRIPTION 

QRTHETA  uses  the  relationships  of  sperical  trigonometry  to  calculate 
the  point  where  a map  radial  at  bearing  AZIM  intersects  the  latitude/longitude 
boundaries  of  the  map. 

B.2.3.11  QTILE 

Subroutine  QTILE  partially  reorders  the  elements  of  an  array  so  that 
A(J),  J=l,  IR  are  all  greater  than  or  equal  to  all  A(J),  J=1R,  NN.  In 
particular  A(IR)  will  have  the  same  value  it  would  have  if  the  array  were 
completely  sorted  in  descending  order. 

ARGUMENTS 

NN  - Number  of  elements  of  the  array  to  be  ordered. 

A - Array  to  be  ordered. 

IR  - Array  elements  are  to  be  sorted  so  that  A(J),  J=l,  IR  are  all 
greater  than  or  equal  to  A(J),  J=IR,  NN. 

COMMON  BLOCKS 
None. 

ROUTINES  CALLED 
None. 

CALLED  BY 

QTILE  is  called  by  QDLTH. 
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SUBROUTINE  DESCRIPTION 


QTILE  is  called  to  partially  sort  the  elements  of  array  A.  Elements 
are  interchanged  until  the  first  IR  elements  are  the  largest.  This  routine 
is  used  to  order  values  for  QDLTH  so  that  it  can  determine  the  10th  and 
90th  quantiles. 


B.2.3.12  READPRF 

Subroutine  READPRF  is  called  to  create  the  FIT  (File  Information  Table) 
for  reading  the  TOPO  tape. 

ARGUMENTS 

FIT  - Array  containing  the  record  manager  FIT  and  I/O  buffer  for 
reading  the  TOPO  tape. 

LFN  - Logical  file  name  of  the  TOPO  tape. 

MRLW  - Maximum  record  length  in  words  for  the  TOPO  tape. 

COMMON  BLOCKS 

READPRF  uses  the  common  block  /PREFACE/  described  in  the  documentation 
for  LDDTPO. 

ROUTINES  CALLED 

READPRF  calls  record  manager  routines  to  set  fields  in  the  FIT. 

CALLED  BY 

READPRF  is  called  by  LDDTPO. 

ENTRY  POINTS 

PRNTPRF  is  called  to  print  the  information  contained  in  the 
PREFACE  record  on  the  TOPO  tape. 

SUBROUTINE  DESCRIPTION 

READPRF  calls  FILESQ  to  set  up  a FIT  for  the  sequential  tape  file 
TOPO.  It  opens  the  file  for  INPUT  and  stores  the  maximum  record  length 
in  the  FIT.  It  uses  a 1065  word  buffer  beginning  at  FIT(36).  READPRF 
reads  the  PREFACE  record  from  the  TOPO  tape. 


B.2.3.13  Longley-Rice  Propagation  Loss  Model  Routines 

The  routines  ACRLR,  ADHLR,  AKNFLR,  ALOSLR,  DEELR,  EIGLER,  FHZLR,  HOFLR, 
SCATLR,  QLRAP,  and  QVARLR  are  also  a part  of  the  (2,0)  overlay.  These 
routines  with  XLBLR  in  the  (3,0)  overlay  comprise  the  Longley-Rice  propaga- 
tion loss  model  used  in  the  point-to-point  mode  by  the  CSPM.  The  Longley- 
Rice  model  is  discussed  briefly  in  Appendix  C of  this  report  and  documented 
in  the  reports  by  Longley  and  Rice  (1968)  and  Rice,  et  al.  (1967)  to  which 
the  reader  is  referred  for  additional  detail. 

B.2.4  (3,0)  Overlay  Routines 

B . 2 . 4 . 1 OUTPUT 

Program  OUTPUT  is  the  controlling  routine  for  the  (3,0)  overlay. 

It  reads  both  PROFILE  and  ACR  DATA  records  from  the  LOSS  or  RESTRT  files. 
OUTPUT  calls  PRTAB  to  produce  the  requested  path  profile  output.  If  this 
is  a single  path  and  basic  transmission  loss  (LB)  has  been  requested,  OUTPUT 
calculates  and  prints  the  basic  transmission  loss  for  the  path.  For  a 
geographic  area,  OUTPUT  uses  the  ACR  and  statistical  data  read  for  each 
radial  to  calculate  basic  transmission  loss  and  communication  reliability 
values  at  each  data  point  along  the  radial.  These  values  are  used  by  PRLOSS 
to  print  the  requested  tabulated  output  and  are  written  to  file  CNTR  to 
be  read  by  CONTOUR  to  produce  the  selected  contour  plots. 

COMMON  BLOCK 

OUTPUT  uses  the  /OPTIONS/,  /LOSSMOD/,  /LRVAR/,  and  /LBREL/  common 
bloc)cs.  These  are  described  in  Tables  B-4,  B-5,  B-14,  and  B-16,  respectively. 

ROUTINES  CALLED 

OUTPUT  calls  the  following  routines: 

CONTOUR 

HEADER 

PLOTF  (CALCOMP  initialization  routine) 

PQFF  (CALCOMP  termination  routine) 

PRLOSS 

PRTAB 
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QERF 

QERFI 

XLBLR 

SMODE  (CALCOMP  routine) . 

CALLED  BY 

OUTPUT  is  called  by  CSPM. 

PROGRAM  DESCRIPTION 

If  this  is  an  initial  run,  the  PROFILE  and  ACR  DATA  records  for  OUTPUT 
to  read  are  on  file  LOSS.  To  properly  position  the  LOSS  file,  OUTPUT  rewinds 
the  file  and  skips  the  OPTIONS  record.  If  this  is  a restart  run,  the  PROFILE 
and  ACR  DATA  records  are  read  from  RESTRT.  Since  INPUT  will  already  have 
read  the  OPTIONS  record,  the  RESTRT  file  does  not  need  to  be  repositioned. 
OUTPUT  initially  calls  QERFI  to  calculate  the  standard  normal  deviates 
for  the  variability  sets  used.  It  also  calculates  the  constants  that  are 
added  to  basic  transmission  loss  values  to  calculate  power  density  and 
received  signal  level.  If  a geographic  area  is  being  analyzed  instead 
of  a single  path,  OUTPUT  will  print  the  CSPM  Output  Data  Summary  for  a 
geographic  area.  OUTPUT  then  reads  each  record  from  the  LOSS  or  RESTRT 
file  until  it  finds  an  EOF.  If  the  record  read  is  a PROFILE  record,  PRTAB 
is  called  to  tabulate  or  plot  the  profile  data.  If  the  record  read  is 
an  ACR  DATA  record  and  the  run  is  for  a single  path,  OUTPUT  calculates 
and  prints  the  basic  transmission  loss  for  the  path.  For  a geographic 
area,  OUTPUT  uses  the  ACR  and  statistical  data  for  the  radial  to  calculate 
basic  transmission  loss  and  communication  reliability  values  for  each  point 
along  the  radial.  If  tabulated  output  has  been  selected  for  this  radial, 
PRLOSS  is  called  to  print  LB,  PD,  RSL,  and  REL.  If  contour  plots  have 
been  requested,  these  data  are  written  to  file  CNTR.  When  all  the  records 
on  file  CNTR  have  been  read,  OUTPUT  checks  the  plot  flag  to  see  if  contour 
plots  have  been  requested.  If  they  have,  CONTOUR  is  called  to  read  the 
CNTR  file  and  produce  the  CALCOMP  plots. 


I 

I 
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B.2.4.2  CONTOUR 


Subroutine  CONTOUR  produces  all  of  the  contour  plots. 

ARGUMENTS 

XLMAX  - Array  of  maximum  loss  value  for  each  variability  set. 

XLMIN  - Array  of  minimum  loss  value  for  each  variability  set. 

PDC  - The  constant  for  computing  power  density  from  loss  data. 

RSLC  - The  constant  for  computing  received  signal  level  from  loss  data. 

ROUTINES  CALLED 
PTITLE . 

■ ) 

Plotting  Routines. 

CALLED  BY 

S 

OUTPUT. 


LOCAL  VARIABLES 

HEIGHT  - The  height  of  the  page  to  contain  a plot  in  plotter  units. 
WIDTH  - The  width  of  the  page  to  contain  a plot  in  plotter  units. 

SIZE  - The  vertical  size  of  a plotted  letter  in  plotter  units. 

WL  - The  width  of  a letter  in  plotter  units. 

TYPE  - Array  of  specified  plot  types 

1 = Basic  Transmission  Loss 

2 = Power  Density 

3 = Received  Signal  Level 

4 = Reliability 

NTYPE  - The  number  of  types  of  plots  in  TYPE  array. 

C - Array  of  levels  to  be  contoured. 

NC  - Number  of  levels  to  be  contoured. 

XLOSSl  - Array  of  still  active  loss  data. 

XRELl  - Array  of  still  active  reliability  data. 

XL0SS2  - Array  of  loss  data  from  radial  which  was  just  read  in. 

XREL2  - Array  of  reliability  data  from  radial  which  was  just  read  in. 
COOR  - Linked  array  of  coordinates  for  contour  lines. 

P - Array  of  pointers  into  COOR  list. 

NCH  - Number  of  groups  of  points  which  form  continuous  lines  in  COOR. 


128 


SUBROUTINE  DESCRIPTION 


Subroutine  CONTOUR  checks  the  flags  for  contour  plots  and  produces  all 
specified  plots  in  a single  call.  PTITLE  is  called  for  each  plot  to  calculate 
the  levels  to  be  contoured,  and  to  write  a title  and  contour  level  key  on 
the  plot. 

Loss  and  Reliability  data  for  the  plots  is  read  from  logical  unit  10 
one  radial  at  a time.  It  calculates  the  endpoints  of  all  possible  contour 
lines  between  the  last  two  radials  before  reading  in  the  next  one.  The  end- 
points are  calculated  by  linear  interpolation  between  data  points.  The  end- 
points are  stored  as  linked  lists  where  each  list  forms  a continuous  contour 
line. 

Once  all  the  radials  have  been  read  in,  the  contours  are  drawn.  A mark 
is  put  at  the  transmitter  location  with  a circle  of  radius  1 km  drawn 
around  it.  CONTOUR  then  goes  on  to  the  next  plot. 

B . 2 . 4 . 3 PRLOSS 

Subroutine  PRLOSS  prints  tabulated  CSPM  evaluation  data  for  selected 
radial  paths.  Any  combination  of  basic  transmission  loss,  power  density, 
received  signal  level,  and  communication  reliability  values  may  be  printed 
for  all  variability  sets. 

ARGUMENTS 

2 

PDC  Power  density  constant  used  to  calculate  PD  in  dBm/m 

from  basic  transmission  loss,  PD  = PDC  - LB. 

RSLC  Received  signal  level  constant  used  to  calculate  RSL  in 

dBm  from  basic  transmission  loss,  RSL  = RSLC  - LB. 

XDLTH  Terrain  irregularity  in  meters  for  this  radial . 

XENS  Surface  ref r activity  for  this  radial. 

DIST  Path  length  in  meters. 


COMMON  BLOCKS 

PRLOSS  uses  conmon  block  /OPTIONS/  described  in  Table  B-4. 


CALLED  BY 


PRLOSS  is  called  by  OUTPUT. 

SUBROUTINE  DESCRIPTION 

PRLOSS  prints  a short  summary  heading  called  CSPM  Evaluation  Data 
which  gives  the  radial  bearing,  length  in  meters,  and  path  parameter  informa- 
tion. Next,  the  column  headings  for  the  requested  types  of  output  (LB, 

PD,  RSL,  REL)  are  printed.  Column  headings  for  each  specified  variability 
set  are  also  printed.  For  each  data  point  along  the  radial,  the  requested 
power  density  and  received  signal  level  values  are  calculated  from  LB. 

Then  the  values  of  LB,  PD,  RSL,  and  REL  for  that  subpath  are  tabulated 
together  with  the  subpath  length  in  kilometers.  If  the  data  must  be  con- 
tinued on  subsequent  pages,  PRLOSS  calls  HEADER  to  title  the  new  pages 
and  prints  the  appropriate  column  headings  on  each. 

B.2.4.4  PROFILE 

Subroutine  PROFILE  draws  the  path  profiles. 

ARGUMENTS 

PFL  An  array  of  profile  elevations  in  meters. 

NP  The  number  of  points  in  PFL  array. 

D The  length  of  the  path  in  meters. 

THETA  The  bearing  of  the  path  (if  a map  radial) . 

EERC  The  effective  earth  radius  for  the  path. 

CALLING  ROUTINES 

PRTAB  calls  it  for  each  profile. 

ROUTINES  CALLED 

XAXIS  to  draw  the  horizontal  axes. 

Plotting  routines. 

LOCAL  VARIABLES 

ZSF  Vertical  scale  factor. 

ZINIT  Minimum  elevation  of  vertical  eucis  in  meters. 

DINC  Distance  increment  along  horizontal  axis  in  plotter  units. 

XINC  Distance  increment  along  horizontal  axis  in  kilometers. 

ZAXIS  Height  of  vertical  axis  in  plotter  units. 

Length  of  horizontal  axis  in  plotter  units. 

130 


XLNGTH 


SUBROUTINE  DESCRIPTION 


1 


E 


PROFILE  writes  the  run  title  in  the  upper  right-hand  corner  of  the 
plot  and  the  data  and  time  in  the  upper  left-hand  corner.  Across  the 
top  it  writes  a title  which  depends  on  whether  it's  a single  path  or  a 
map  radial.  It  writes  the  transmitter  and  receiver  (single  path  only) 
coordinates  in  degrees,  bearing  (map  radial  only),  the  path  length  in 
Itilometers  and  miles,  and  the  effective  earth  radius  in  kilometers.  If 
the  size  of  the  earth's  bulge  for  the  path  is  larger  than  two-thirds  of 
the  maximum  difference  in  elevation  along  the  path,  then  a flat  earth  is 
drawn.  For  a curved  earth  a parabolic  approximation  is  used.  Vertical 
axes  are  in  both  meters  and  feet.  Horizontal  axes  are  in  both  kilometers 
and  miles. 


j 

i 


! 


i 

\ 
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B . 2 . 4 . 5 PRTAB 

Subroutine  PRTAlJ  prints  tabulated  path  profile  data  for  single  paths 
or  selected  radial  paths  within  a geographic  area.  It  also  prints  the 
heading  information  for  plotted  profiles  and  calls  PROFILE  to  produce  the 
CALCOMP  plots. 


ARGUMENTS 

IPLT  Plot  or  tabulate  flag.  If  IPLT  = -1,  the  heading  is 

printed  and  PROFILE  is  called  to  plot  the  path  profile. 
If  IPLT  = 0,  the  heading  is  printed  and  the  tabulated 
path  profile  data  are  printed.  If  IPLT  = 1,  the  heading 
is  printed,  the  ted^ulated  path  profile  data  are  printed, 
and  PROFILE  is  called  to  plot  the  profile. 

AZIM  Bearing  decimal  degrees  clockwise  from  north. 

XDLTH  Delta-H  (terrain  irregularity)  in  meters  for  this  path. 

ENS  Surface  refractivity  for  the  path. 

HSE  HSE(l)  transmitter  effective  surface  elevation  in 

meters  above  MSL. 

HSE (2)  receiver  effective  surface  elevation  in  meters 
above  MSL. 

DL  DL(1)  >>  distance  in  meters  from  the  transmitter  to  its 


horizon. 
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DL(2)  = disteince  in  meters  from  the  receiver  to  its 
horizon. 

HL  HL(1)  = transmitter  horizon  elevation  in  meters. 

HL(2)  = receiver  horizon  elevation  in  meters. 

D Path  length  in  meters. 

XI  Distance  increment  in  meters  between  profile  points. 

NP  Number  of  profile  points  - not  including  the  receiver 

location. 

PFL  Array  containing  NP+1  profile  elevations  in  meters. 

COMMON  BI/XKS 

PRTAB  uses  the  /OPTIONS/  common  block  described  in  Table  B-4. 

ROUTINES  CALLED 

PRTAB  calls  the  following  routines: 

HEADER. 

PROFILE. 

CALLED  BY 

PRTAB  is  called  by  OUTPUT. 

SUBROUTINE  DESCRIPTION 

Subroutine  PRTAB  calls  HEADER  to  title  a new  page  and  then  checks 
to  see  if  the  profile  data  to  be  tabulated  are  for  a single  path  or  a radial 
path.  If  it  is  a single  path  case,  PRTAB  prints  the  CSPM  Output  Data  Summary 
for  a single  path.  For  a radial  path,  PRTAB  prints  an  abbreviated  heading 
called  Topographic  Profile  Data.  Both  contain  information  about  the  path 
length  and  path  parameters.  The  profile  elevations  and  the  distance  of 
the  point  from  the  transmitter  are  printed  in  both  English  and  International 
units.  If  more  than  one  page  is  required  to  tabulate  all  the  profile  eleva- 
tions, PRTAB  calls  HEADER  to  title  the  page  and  also  to  print  the  column  ' 

headings  on  each  continuation  page.  When  all  data  have  been  printed,  the  PRO 
flag  is  checked  to  see  if  a plotted  profile  has  been  requested.  If  so, 

PROFILE  is  called  to  create  the  CALCOMP  plot.  i 

I 
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B.2.4.6  PTITLE 


Subroutine  PTITLE  writes  the  title  and  contour  level  key  on  the 
contour  plots.  It  also  computes  the  contour  levels. 

ARGUMENTS 
XLMAX 
XLMIN 
PDC 
RSLC 

TYPE 


T 

V 


HEIGHT 
WIDTH 
SIZE 
WL 

LETTER 
C 
NC 

yc2 

COMMON  BLOCKS 

PTITLE  uses  the  common  block  /OPTIONS/  described  in  Table  B-4. 

ROUTINES  CALLED 

Plotting  routines. 
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The  maximum  loss  values  for  each  variability  set. 

The  minimum  loss  values  for  each  variad>ility  set. 

The  constant  for  computing  power  density  from  loss  data. 

The  constant  for  computing  received  signal  level  from  loss 
data . j 

Integer  array  containing  the  contour  plots  which  were 
requested: 

1 = Basic  Transmission  Loss 

2 = Power  Density 

3 “ Received  Signal  Level 

4 = Communication  Reliability. 

Index  to  the  TYPE  array  (integer) . 

Which  variability  set: 

1 = VI 

2 = V2 

3 = V3. 

Height  of  page  containing  the  plot. 

Width  of  page  containing  the  plot. 

Vertical  size  of  plotted  character. 

Width  of  a plotted  character. 

Array  of  letters  for  contour  level  key. 

Array  of  contour  levels  for  CONTOUR. 

Number  of  contour  levels. 

Y-coordinate  of  upper  edge  of  contour  map. 


CALLING  ROUTINES 


CONTOUR. 

LOCAL  VARIABLES 

CONTR  Array  of  contour  levels  for  the  key. 

SUBROUTINE  DESCRIPTION 

PTITLE  writes  the  run  title  in  the  upper  right-hand  corner  of  the 
plot  and  the  date  and  time  in  the  upper  left-hand  corner.  Across  the  top 
it  writes  the  type  of  contours  being  plotted  and  the  variability  set  which 
was  used  to  compute  the  data  for  the  plot.  After  determining  the  levels 
to  be  contoured,  PTITLE  writes  a key  at  the  bottom  of  the  plot.  For  the 
Basic  Transmission  Loss,  Power  Density,  and  Received  Signal  Level  contours 
the  levels  are  computed  so  that  they  will  be  evenly  spaced  multiples  of  5 
with  a maximum  of  ten  levels.  For  these  three  types  of  plots  CONTOUR  does 
the  plotting  from  basic  transmission  loss  data  so  PTITLE  sends  it  the 
corresponding  contour  levels.  The  reliability  contour  levels  are  specified 
by  the  user.  PTITLE  removes  any  duplicates  and  sorts  them  into  descending 
order . 

B . 2 . 4 . 7 QERF 

Function  QERF  converts  a standard  normal  deviate  to  a probability 
value. 

ARGUMENTS 

X Standard  normal  deviate  for  time  availability  (ZT) , location 
variability  (ZL) , or  confidence  (ZC) . 

ROUTINES  CALLED 
None . 

CALLED  BY 

QERF  is  called  by  OUTPUT. 

FUNCTION  DESCRIPTION 

QERF  uses  a Chebyshev  approximation  to  the  standard  normal  complementary 
probability  due  to  C.  Hastings,  Jr.,  in  "Approximations  for  Digital  Com- 
puters," Princeton  University  Press,  1955.  It  has  a maximum  error  of  7.5 
k E-8.  To  avoid  underflow  QERF  is  truncated  at  |x|  = 10.  This  is  the 

I 13^ 


"q  error  function"  returning  the  probability  that  the  value  X is  exceeded 
under  a standard  normal  distribution. 

B.2.4.8  QERFI 

Function  QERFI  converts  a probability  value  to  a standard  normal  deviate. 
It  is  the  inverse  of  QERF. 

ARGUMENTS 

Q - Probability  value  for  time  availcibility  (QT) , location 
variability  (QL) , or  confidence  (QC) . 

ROUTINES  CALLED 
None . 

CALLED  BY 

QERFI  is  called  by  OUTPUT. 

FUNCTION  DESCRIPTION 

QERFI  returns  the  standard  normal  deviate  as  a function  of  the  complemen- 
tary probability  truncated  at  0.00001  and  0.99999.  It  uses  the  approximation 
due  to  C.  Hastings,  Jr.,  in  "Approximations  for  Digital  Computers,"  Princeton 
University  Press,  1955,  and  has  a maximum  error  of  4.5  E-4. 

B.2.4.9  XAXIS 

Subroutine  XAXIS  draws  the  two  horizontal  axes  on  the  path  profile. 

The  axes  are  curved  according  to  the  effective  earth  radius. 

ARGUMENTS 

2 The  Z-coordinate  of  the  starting  point  of  the  axis.  The 

X-coordinate  is  an  implied  zero. 

LABEL  A character  string  title  for  the  axis. 

NL  The  number  of  characters  in  the  title  if: 

NL>0  - title  will  be  written  above  the  axis, 

NL<0  - title  will  be  written  below  the  aucis. 


AXLEN  The  length  of  the  axis. 

AXINC  The  distance  between  tick  marks. 


DINC  The  increment  for  the  data  values  to  be  written  by  the 

tick  marks.  The  first  value  is  always  zero. 

F A constant  for  computing  the  amount  of  curvature. 

COMMON  BLOCKS 
None. 

ROUTINES  CAI.LED 

Plotting  routines. 

CALLING  ROUTINES 
PROFILE. 

SUBROUTINE  DESCRIPTION 

XAXIS  draws  a curved  cixis  according  to  the  factor  "F" . "F"  is 

computed  for  a parabolic  approximation  of  the  curvature  of  the  earth  using 
its  effective  radius.  The  axis  and  tick  marks  are  labeled  as  specified 
by  the  parameters. 

B . 2 . 4 . 10  XLBLR 

Function  XLBLR  calculates  basic  transmission  loss  (LB)  from  the  ACR 
values  input  in  the  common  block  /LRVAR/.  The  value  returned  is  the  loss 
that  is  not  exceeded  with  confidence  QC  in  at  least  QL  of  the  locations 
for  at  least  QT  of  the  time. 

ARGUMENTS 

ZT  - The  standard  normal  deviate  for  time  availability  (QT) . 

ZL  - The  standard  normal  deviate  for  location  variability  (QL) . 

ZC  - The  standard  normal  deviate  for  confidence  (QC) . 

YT  - Variability  in  attenuation  due  to  long-term  fading.  This  value 
is  returned. 

YL  - Variability  in  attenuation  due  to  random  path  selection.  This 
value  is  returned. 

YC  - Variability  in  attenuation  due  to  prediction  confidence.  This 
value  is  returned. 

COMMON  BLCCKS 

XLBLR  uses  comroon  block  /LRVAR/  defined  in  Table  B-14. 
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ROUTINES  CALLED 


None. 

CALLED  BY 

XLBLR  is  called  by  OUTPUT. 

FUNCTION  DESCRIPTION 

XLBLR  calculates  the  basic  transmission  loss  for  the  path  as  the  sum 
of  the  free-space  loss  and  the  calculated  reference  modified  by  the  variabil- 
ity due  to  long-term  fading,  random  path  selection,  and  prediction  confidence. 

B.3  CSPM  Installation  Instructions  for  a SCOPE  3.4  System 

The  CSPM  source  is  being  delivered  as  an  UPDATE  PL  on  a 7-track,  800 
bpi  density,  unlabeled,  SI  format  tape.  The  following  control  cards  should 
be  used  to  create  and  execute  the  CSPM  overlays. 

CSPM,T77,MT1. 

REQUEST (NEWPL, *PF) 

REQUEST (CSPM, *PF) 

REQUEST , TAPE , VSN=CSPM , HY , MT . 

UPDATE ( P=TAPE , F , N, *=/) 

CATALOG (NEWPL, CSPMPL, ID=XXX) 

FTN(I,R=3) 

ATTACH (CLIB,ID=XXX)  Library  of  CALCOMP,  routines  plotter 

MAP(P) 

RFL, 100000. 

LDSET (PRESET=NGINDEF) 

LDSET ( LIB=CLIB/FORTRAN/SYS 10) 

LDSET,USE=$CLSV.WA$/$GET.FIT$/$REW.SQ$. 

LDSET,USE»$OPEN.WA$/$PUT.WA$/$GET.WA$/$CLSF.WA$. 

LDSET, USE=  $GET . W$/$GET . R$/$GET . D$ . 

LDSET, USEP-$INCOM=$/SYSTEM/SySTEMC . 

LOAD(LGO) 

NOGO. 

CATALOG (CSPM, ID=XXX) 

RETURN (TAPE) 
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LABEL (TOPO,MT,HY,VSN=ZZZZ)  Topographic  data  tape 
CSPM. 

7/8/9 


CSPM  TEST  CASE 

XMTRCOOR= ( 31 , 29 , 00 , N/110 , 20 , 00, W) , XMTRLOC=HUACHUCA , 
MAP=(31,27,00,N/31,34,00,N/110,26,30,W/110,17,30,W) , TABL=16, PLOT=0 
FREQ=394.675,POL=V,HG1=20,HG2=6FT,LB,PD,RSL,REL,PTAB=16 
Vl=(.9,.9, .9) ,V2=(.9, .9, .5) ,V3=(.9, .5, .5) 

RSLTH=-85.5,RQ=( .99, .90, .50, .10) ,G2=10 

GND=PCX)R,PWR=100W,G1=15 

6/7/8/9 


APPENDIX  C.  LONGLEY-RICE  PROPAGATION  LOSS  MODEL 
C.l  General  Summary 


The  Longley-Rice  propagation  loss  model  is  the  best  general  purpose 
model  available  at  the  Institute  for  Telecommunication  Sciences  for  pre- 
dicting long-term  (hourly)  median  radio  transmission  loss  at  VHF  and  higher 
frequencies  over  irregular  terrain.  The  method  is  based  on  well-established 
propagation  theory  which  has  been  documented  in  considerable  detail  (see  Rice 
et  al.,  1967;  Longley  and  Rice,  1968;  and  Longley,  1976).  The  model  has  been 
tested  against  a large  number  of  propagation  measurements  and  used  exten- 
sively to  provide  information  cibout  curea-wide  situations  employing  line-of- 
sight,  diffractior',  and  troposcatter  types  of  communication  systems.  Ihe 
model  frequently  is  used  to  analyze  broadcast  and  mobile  system  applications 
where  only  moderate  detail  is  )cnown  about  the  propagation  path.  For  such 
analyses,  statistical  descriptors  of  terrain  characteristics  (i.e.,  terrain 
irregularity,  surface  refractivity , etc.)  are  used.  The  CSPM,  however,  uses 
the  Longley-Rice  propagation  loss  model  in  the  point-to-point  mode.  In  this 
mode,  a particular  propagation  path  may  be  considered  for  which  a path 
terrain  profile  is  generated  and  used  to  provide  the  details  of  terrain 
char ;c1 eristics  along  the  propagation  path.  That  is,  the  path  is  determined 
to  be  line-of-sight,  diffraction,  or  troposcatter;  and  appropriate  influences 
upon  the  radio  signals  are  supplied  mathematically. 

The  prediction  method  has  been  developed  for  operation  on  a high  speed, 
digital  computer,  whereby  predictions  have  been  tested  against  data  for  wide 
ranges  of  frequency,  antenna  height  and  distance,  and  for  all  types  of 
terrain  from  very  smooth  plains  to  extremely  rugged  mountains.  The  data 
base  includes  more  than  500  long-term  recordings  throughout  the  world  in 
the  frequency  range  40  to  10,000  MHz,  and  several  thousand  mobile  record- 
ings in  the  United  States  at  frequencies  from  20  to  1,000  MHz.  Some 
experience  has  shown  that  the  model  will  overestimate  the  loss  for  good 
line-of-sight  paths. 

The  "heart"  of  the  propagation  loss  model  is  the  calculation  of  median 
values  of  reference  attenuation  relative  to  free-space  loss  as  a function 
of  distance.  This  modification  to  the  free-space  loss  ta)ces  into  considera- 
tion whether  the  path  is  line-of-sight,  diffraction,  or  tropospheric  scatter. 
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For  radio  line-of-sight  paths,  the  calculated  reference  attenuation 
factor  is  based  on  two-ray  theory  and  an  extropolated  value  of  diffraction 
attenuation.  For  transhorizon  paths,  the  reference  value  is  either  diffrac- 
tion attenuation  or  forward  scatter  attenuation,  whichever  is  smaller. 

Further  considerations  due  to  the  influences  of  season  and  climate  upon 
propagation  loss  also  are  included. 

Without  defending  derivation  of  the  expression,  which  is  contained 
in  the  references  cited  above,  an  expression  Cein  be  written  for  basic  trans- 
mission loss,  with  quantities  expressed  in  decibels,  as 

L = + A - V (C-1) 

D fs  cr 

where : 

is  the  calculated  free-space  loss, 

A is  the  calculated  reference  attenuation  factor  relative  to  free- 
cr 

space  loss  and  depends  on  whether  the  path  is  line-of-sight, 
diffraction,  or  tropospheric  scatter,  and 
V is  an  empirically  determined  adjustment  to  the  calculated 

reference  attenuation  factor  to  provide  a median  attenuation 
factor,  denoted  A(0.5),  for  specific  seasons  and  climates. 

Median  basic  transmission  loss  can  be  expressed,  then,  as 

U (0.5)  = L-  + A - V(0.5)  (C-2) 

D fs  cr 

= + A(0.5) 

f s 

where  the  (0.5)  denotes  median  conditions.  Discussion  on  the  time,  location, 
and  prediction  confidence  variabilities,  to  allow  calculations  of  basic 
transmission  loss  for  other  than  median  conditions,  is  contained  in  the 
next  section  of  this  appendix.  Other  terms  in  (C-1)  will  be  discussed 
briefly  in  the  remainder  of  this  section. 

The  general  expression  for  free-space  propagation  loss  is 

L,  = 10  log,-  ((4Trfd)/c)^,  (C-3) 

f S 10 

where  d is  the  distance  (path  length) , f is  the  carrier  frequency  of  the 
radiated  signal,  and  c is  the  propagation  velocity  in  free  space.  When 
frequency  is  expressed  in  megahertz  and  distance  is  expressed  in  )cilometers, 
the  free-space  loss  expression  becomes 
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= 32.45  + 20  log,„(fd) 
fs  10 


(C-4) 


As  indicated  earlier,  the  reference  attenuation  factor  depends  upon 
whether  the  path  is  line-of-sight,  diffraction,  or  tropospheric  scatter. 
The  graph  in  Figure  C-1  shows  how  varies  with  distance. 
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Figure  C-1.  Calculated  attenuation  factor  A^,^  versus  path 
length  for  line-of-sight,  diffraction,  and 
tropospheric  scatter  regions. 


In  the  free-space  region  A^^  is  0 since  there  is  no  attenuation  above 
free-space  loss.  For  a line-of-sight  path,  where  a weighted  average  of 
two-ray  optics  theory  and  diffraction  theory  is  applied,  the  reference 
attenuation  is  calculated  as 

‘los  ■ *0  ♦ * "2 

where 


Aq  is  a weighted  average  attenuation  above  free- space  loss  computed 
at  a distance  d^  using  techniques  defined  by  Longley  and  Rice, 
1968; 

dg  is  a distance  selected  to  approximate  the  greatest  distance  at 
which  the  attenuation  above  free-space  loss  is  0 dB,  and 
and  are  constants  determined  by  techniques  defined  by  Longley 
and  Rice,  1968. 

The  reference  attenuation  factor  for  a diffraction  path  is 


A = A + A - m (d  -d) , where  m = (A  -A  )/(d  -d  ) (C-6) 

dfo4d4  d 4343 

and  A^^  is  an  estimate  of  reference  attenuation  relative  to  free-space 

attenuation  due  to  surface  clutter, 

A^  and  A^  are  diffraction  reference  attenuation  factors  computed  at 

distances  d^  and  d^  respectively,  using  the  relationship 

A^  = (l-w)Aj^  + w A^, 

A^  is  an  estimate  of  reference  attenuation  relative  to  free-space 
attenuation  in  the  far  diffraction  region  for  highly  irregular 
terrain, 

A^  is  an  estimate  for  reference  attenuation  relative  to  free-space 
attenuation  in  the  far  diffraction  region  for  smooth  terrain, 
d^  and  d^  are  distances  selected  well  beyond  the  horizon,  using 
techniques  defined  by  Longley  and  Rice,  1968,  at  which 
diffraction  reference  attenuation  factors  are  calculated,  and 
w is  a weighting  factor  empirically  determined  by  techniques  defined 
by  Longley  and  Rice,  1968. 

Reference  attenuation  relative  to  free-space  attenuation  for  a forward 
scatter  path  is 


and 


A = A_  - m d_  + m d,  where  m = (A.-A,. ) / (d,-d^) 
SDSSs  s 6565 


(C-7) 


A-  and  A-  are  factors  of  scatter  reference  attenuation  above  free-space 
5 6 

loss  computed  at  distances  d_  and  d-  respectively,  using 

5 D 

estimating  techniques  defined  by  Longley-rice,  1968;  and 

d and  d,  are  distances  selected  to  be  200  km  and  400  km  greater 
b o 

respectively,  than  the  sum  of  the  distances  from  each  antenna  to 
its  respective  horizon. 
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For  transhorizon  paths,  the  reference  attenuation  relative  to  free- 
space  attenuation  is  defined  to  be  the  diffraction  reference  attenuation 
factor  for  d i d < d and  the  tropospheric  scatter  reference  attenuation 

Ls  X 

factor  for  d < d where  d is  the  sum  of  the  smooth  earth  horizon  distances 
X Ls 

and  d is  the  distance  at  which  diffraction  and  scatter  reference  attenuation 

X 

factors  are  equal. 

The  specification  of  atmospheric  refractivity  is  allowed  as  either  the 
surface  refractivity,  N^,  or  the  surface  refractivity  referenced  to  mean 
sea  level,  N^.  Values  throughout  the  world  of  surface  refractivity 
referenced  to  mean  sea  level  (N^)  may  be  found  in  the  CCIR  (1974)  report 
or  the  monograph  by  Bean  et  al.  (1960).  However,  the  propagation  loss 
model  requires  the  surface  refractivity  (N^) . When  surface  refractivity 
referenced  to  mean  sea  level  is  specified,  it  is  converted  to  surface  refrac- 
tivity as  shown  by  (D-3)  in  Appendix  D.  A commonly-used  value  for  surface 
refractivity  is  = 301  N-units,  which  corresponds  to  an  effective  4/3 
earth  radius . 

C.2  Discussion  of  Statistical  Variabilities 

The  median  basic  transmission  loss,  a function  of  distance,  described 

in  the  preceding  section  is  combined  (by  addition  processes  when  values 

are  expressed  in  decibels)  with  still  other  parameter  values  which  account 

for  the  variabilities  in  transmission  loss  due  to  long-term  fading  (time 

availability),  path-to-path  variations,  and  estimating  confidence.  These 

operations  provide  predicted  basic  transmission  loss  which  will  not  be 

exceeded  for  at  least  the  specified  time  (q^)  over  at  least  the  required 

fraction  of  paths  (q  ) with  stated  confidence  (Q) . The  loss  variability 

L 

paraimeters  commonly  are  referred  to  as  Y-parameters  with  the  time  availa- 
bility influence  denoted  by  Y^,  the  location  variability  influence  for 

specified  time  availability  conditions  denoted  by  Y , and  the  estimating 

L» 

confidence  influence  for  specified  time  availedaility  and  location  variability 
conditions  denoted  by  Y^ . 

The  Y-parcuneters  express  attenuation  variability  in  decibels  and  are 

considered  to  be  approximately  normally  distributed  with  zero  means  and 
2 2 2 

variances  of  a , a , and  a respectively.  They  each  are  written  in  terms 
I*  L ^ 
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of  the  standard  ^^{q)  where  0 < q < 1 and  their  standard  deviations  o^, 

a , and  o . The  standard  normal  deviate  is  expressed  in  terms  of  the  error 

^ ^ -1 
function  exf (x)  and  its  inverse  erf  (x)  as 

Z^(q)  = /2  erf"^(2q-l)  (C-8a) 

and 

q = 0.5  + 0.5  erfCZ^/vT)  . (C-8b) 

The  standard  deviation  for  each  variability  has  been  determined  empirically 
from  measured  data  cited  in  the  report  references  mentioned  earlier  in  this 
appendix . 

The  long-term  time  variability  ediout  the  median  attenuation  is  approxi- 
mately normal  and  calculated  as  follows: 

If  Zp(q^)  <0,  = 0^  ZQ(q^)  (C-9) 

If  ^T  ° 

If  Z^(q^)  > 0 and  Zp(qy)  > ZD  , 

Z„(q  0 + YD 
T TD  0 T 

If  ZQ(qj,)  > 0 and  Z^Cq^)  s ZD  , 

^T  “ °TP  ^O^'^T^ 

where  q^  is  the  fraction  of  time  when  attenuation  on  the  average  will  not 
exceed  L^d.O)  - Y^; 

a__  is  the  standard  deviation  for  the  minus  range; 

TM 

is  the  standard  deviation  for  the  ducting  range; 
o^p  is  the  standard  deviation  for  the  plus  range; 

ZD  is  the  standard  normal  deviate  of  predicted  attenuation  for  small 

time  variabilities,  such  as  are  associated  with  ducting  phenomenon; 
and 

YD  is  the  ducting  range  breakpoint,  ZD  * SGTP. 

Note  that  the  empirically  determined  values  of  . o^,  and  o_„  are  not 
r 1 TM  TD  TP 

constant  for  all  climates  and  effective  path  lengths.  Neither  are  the 
vari2d}ilities  for  q^  < 0.5  symmetrical  with  those  for  q^  > 0.5. 

Path-to-path  variability  in  attenuation  for  given  time  availability 
is  expressed  by 
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(C-10) 


where  q is  the  fraction  of  paths  for  which  attenuation  will  not  exceed 
L 

the  value  calculated  by  L,  (0.5)  - y„  - Y, . The  standard  deviation  o,  has 

b T L L 

been  determined  empirically  to  be  a function  of  terrain  irregularity  and 
frequency  and  is  calculated  by 


10(l-.8  (If) 

13  * (1-.8  (i;) 


(C-11) 


where 


X is  the  wavelength  in  meters  of  the  attenuated  radio  frequency 
signal; 

D is  the  path  length  in  meters;  and 

Ah  is  the  interdecile  range  of  terrain  elevations  in  meters  along 

the  path. 

Variability  due  to  prediction  confidence  for  specified  conditions  of 
time  availability  and  location  variability  is  determined  from  the 
relationship 


2 T 

Y = O + r + 

C C 


7.8  + (Zq(Q))  24  + (Zq(Q))' 


{C-12) 


where  0 < Q < 1 is  the  value  for  prediction  confidence  and  is  the  standard 

deviation  of  the  normal  distribution  of  the  measured  data  for  the  given  q.^, 

and  q^^.  The  value  of  is  calculated  by 

— d 

e 


, 100 

= 3e  +5 


(C-13) 


where  d is  an  effective  distance  in  )tilometers  calculated  as  a function  of 
e 

the  actual  great  circle  path  length,  the  smooth-earth  horizon 
distances,  and  the  difference  in  smooth-earth  horizon  distances 
at  which  diffraction  and  scatter  attenuation  factors  are  equal. 

The  prediction  confidence,  Q,  is  then  the  fraction  of  measured  data  points 
for  which  attenuation  will  not  exceed  the  value  given  by  Lj^(0.5)  “ ~ 
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APPENDIX  D.  DESCRIPTION  OF  THE  PATH  PROFILE  CAPABILITY 


! 


The  capability  is  required  within  the  CSPM  to  generate  propagation 
path  profiles  and  extract  from  them  the  path  parauneter  data  needed  as  input 
to  the  Longley-Rice  propagation  loss  model  operating  in  the  point-to-point 
mode.  If  desired,  the  path  profile  can  be  obtained  as  an  output  from  the 
CSPM  for  either  a single  geodesic  or  selected  radial  paths  extending  (at 
uniform  bearing  angle  increments)  from  a common  origin  into  the  geographical 
area  of  interest. 

The  propagation  path  is  defined  in  one  of  two  ways.  A single  geodesic 
path  is  defined  by  specifying  the  geographical  coordinates  of  the  path 
endpoints;  many  radial  paths  are  defined  by  specifying  the  coordinates 
of  the  origin  and  the  latitudes  and  longitudes  to  form  a geographic  "box" 
containing  the  origin.  (See  Appendix  F for  discussion  on  selection  of 

radial  paths.)  In  this  second  method,  the  CSPM  uses  the  required  number  | 

of  radial  paths  to  compute  the  bearing  angle  increment  between  radial  paths.  | 

I 

For  each  radial  path,  the  model  also  determines  the  intersection  of  the  j 

radial  path  with  the  geographical  area  boundaries.  Endpoints  for  each 
radial  path  then  are  known,  analogous  to  the  input  data  required  for  the 
single  path. 

Output  path  profile  data  for  selected  radial  paths  may  be  obtained 
by  specifying  the  number  (which  must  oe  a power  of  two)  of  radial  paths 
for  which  data  are  desired.  The  radial  paths  are  selected  with  uniform 

bearing  angle  separation  depending  upon  the  number  of  paths,  i.e.,  four 

o o 

paths  provide  90  separation,  eight  paths  provide  45  separation,  etc. 

Knowing  the  path  endpoints,  the  CSPM  accesses  the  digitized  terrain 
data  file  to  obtain  terrain  elevation  data  at  uniform  distance  intervals 
along  the  path.  In  general,  the  points  along  the  path  at  which  elevations 
are  required  will  not  coincide  (geographically)  with  the  points  at  which 
elevation  data  are  stored  in  the  terrain  data  file.  Note  that  these  data 
are  recorded  at  geographic  locations  spaced  by  30"  in  latitude  and  30" 
in  longitude.  This  means  the  data  points  are  separated  by  distances  ranging 
from  etbout  926  m down  to  alsout  583  m,  depending  upon  latitude.  (See  Appendix 
E for  further  details  of  the  topographic  data  file.)  Data  at  four  coordinate 
locations  which  define  the  30"  by  30"  "box"  containing  the  coordinate  point 
of  interest  on  the  propagation  path  are  retrieved  from  the  terrain  file. 

147 


A standard  bilinear  interpolation  is  performed  to  estimate  terrain  elevation 
at  the  coordinate  point  on  the  path  profile. 

Path  profile  point  intervals  along  the  radial  paths  are  250  m,  since 
a variable  distance  interval  would  provide  considerable  difficulty  when 
the  path  profile  data  were  used  to  calculate  basic  transmission  loss  data 
at  1 km  intervals  for  the  plotted  output  options  of  basic  transmission 
loss,  power  density,  received  signal  level,  and  communication  reliability. 

The  fixed  interval  does  mean  that,  in  general,  a data  point  will  not  occur 
at  the  intersection  of  the  radial  path  and  the  area  boundary.  This  situation 
is  overcome  by  generating  the  radial  path  profile  to  a distance  correspond- 
ing with  the  next  integer  number  of  kilometers  beyond  the  intersection. 

The  output  plots  of  contoured  data  simply  are  restricted  to  the  area  within 
the  defined  boundaries.  A further  question  might  arise  about  profile  data 
generated  at  250  m intervals  when  the  topographic  data  are  separated  by 
intervals  ranging  from  about  583  m to  about  926  m.  It  is  evident  that 
the  topographic  data  are  used  redundantly  when  generating  the  path  elevation 
data.  Since,  in  general,  every  elevation  on  the  path  profile  is  an  interpola- 
ted value  however,  it  is  considered  that  some  smoothing  of  the  path  profile 
achieved  through  redundant  use  of  the  basic  terrain  data  is  desirable. 
Intervals  along  a geodesic  are  approximately  250  m for  paths  less  than 
312  km  in  length.  For  paths  greater  than  312  km  in  length,  interval  distance 
along  the  geodesic  is  determined  by  dividing  the  total  path  length  into 
an  integer  number  of  uniform  intervals  with  the  restriction  that  the  total 
niomber  of  points  along  the  path  must  not  exceed  1251. 

The  parameter  values  required  by  the  Longley-Rice  model  in  the  point- 
to-point  mode  include: 

Elevation  of  and  distance  to  the  horizon  "seen"  by  each  antenna  from 
its  path  endpoint. 

Terrain  irregularity  for  the  path. 

Effective  surface  heights  of  terrain  in  the  foreground  of  each  antenna. 

Atmospheric  refractivity  near  the  surface  of  the  earth  (frequently 
termed  surface  refractivity) , 

The  effective  earth's  radius  considering  the  surface  refractivity 
for  the  radial  being  considered. 


These  parameters  are  computed  for  each  path  (each  point  in  the  area  linked 
with  the  origin  constitutes  a path,  therefore  many  paths  can  occur  for 
a single  radial)  as  part  of  the  path  profile  and  parameter  generating  capabil 
ity  of  the  CSPM.  The  path  parcuneter  data  are  used  directly  by  the  Longley- 
Rice  model  with  no  provision  for  the  data  being  output  from  the  CSPM. 

The  horizon  "seen"  by  each  antenna  is  determined  first  by  calculating 
the  elevation  angle  between  the  transmitting  antenna  and  receiving  antenna. 
Defining  the  elevation  angle  to  be  0,  we  have 

= -e^  , (D-i) 

where  the  subscripts  1 and  2 denote  the  transmitter  and  receiver  respectively 
To  determine  the  horizon  for  the  transmitting  antenna,  the  next  step  is 
to  calculate  the  elevation  angle  for  each  successive  point,  beginning  at 
the  first  point  beyond  the  transmitter,  along  the  path  profile.  If  each 
calculation  of  elevation  angle  produces  a value  less  than  0^^,  the  path 
is  line-of-sight.  When  elevation  angles  larger  than  0^^  are  encountered, 
the  point  corresponding  to  the  largest  value  for  the  elevation  angle  is 
defined  to  be  the  horizon.  The  elevation  associated  with  that  point,  of 
course,  is  the  horizon  elevation,  and  the  great  circle  distance  between 
the  point  and  the  transmitting  antenna  location  is  the  distance  to  the 
horizon.  Working  from  the  receiving  antenna  location  and  following  the 
same  steps,  the  horizon,  horizon  elevation,  and  distance  to  the  horizon 
for  the  receiving  antenna  are  determined. 

The  terrain  irregularity  Ah  (expressed  in  meters)  is  computed  once 
for  each  entire  radial  path  as  the  interdecile  range  of  terrain  elevations 
at  each  point  along  the  path.  Terrain  irregularity  for  the  subpaths  of 
a radial  path  is  computed  using  the  relationship 

Ah(d)  = Ah  [1  - 0.8  exp(-0.02d)]  , (D-2) 

where  Ah(d)  is  the  terrain  irregularity  in  meters  for  the  subpath  and  d 
is  the  subpath  length  in  kilometers. 

The  effective  surface  height  for  terrain  in  the  foreground  of  each 
antenna  must  be  provided  to  the  propagation  loss  model.  To  do  so,  the 
path  is  divided  into  segments  when  transhorizon  propagation  will  occur. 
Transmitting  antenna  location  to  the  horizon  for  that  antenna  is  a path 
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segment,  and  receiving  antenna  location  to  the  horizon  for  the  receiving 
antenna  is  another  segment.  The  entire  path  is  considered  for  line-of- 
sight  conditions.  The  effective  terrain  elevation  for  each  segment  or 
the  entire  path,  if  appropriate,  is  calculated  as  the  average  elevation 
for  the  middle  80%  of  the  path  segment  or  entire  path. 

A CSPM  user  may  have  measurements  of  the  minimum  monthly  mean  surface 
refractivity,  denoted  by  N^,  for  the  propagation  paths  he  wishes  to  study. 

On  the  other  hand,  he  may  need  to  rely  on  use  of  maps  showing  the  parameter 
minimum  monthly  surface  refractivity  referenced  to  mean  sea  level,  denoted 
by  Nq,  and  sometimes  termed  "reduced  surface  refractivity."  A specification 
of  surface  refractivity  referenced  to  mean  sea  level  is  converted  to  surface 
refractivity  by  the  relationship 

N = exp  (-0.1057  h ),  (D-3) 

s 0 s 

where  h is  path  elevation  in  kilometers  cibove  mean  sea  level  and  is  computed 
s 

in  the  CSPM  as  the  average  value  for  terrain  elevations  at  the  transmitting 
and  receiving  antenna  locations  (the  path  endpoints) . The  surface  refrac- 
tivity is  used  to  compute  an  effective  earth's  radius  which  allows  for 
regional  differences  in  average  atmospheric  conditions.  The  effective 
earth's  radius  in  kilometers  is  defined  as 

a = 6370  [1-0.04665  exp  (0.005577  N^)]  (D-4) 

where  the  actual  earth's  radius  is  taken  to  be  6370  km. 
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APPENDIX  E.  DISCUSSION  OF  THE  TOPOGRAPHIC  DATA  BASE 


The  topographic  data  base,  maintained  on  magnetic  tape,  contains  digitized 

terrain  elevations  for  the  geographical  area  bounded  by  23°N  to  51°N  latitude 
o o 

and  60  W to  130  W longitude — the  continental  United  States  (CONUS)  essentially. 
(Data  for  a portion  of  western  Europe  are  expected  to  be  added  to  the  data 
base  in  the  near  future.)  The  data  base  provides  terrain  elevation  eibove 
mean  sea  level  at  every  30"  of  latitude  and  30"  of  longitude.  The  data 
originally  were  read  from  1/250,000  scale  maps  developed  by  the  Defense 
Mapping  Agency  (DMA) , then  adapted  by  the  Electromagnetic  Compatibility 
Analysis  Center  (ECAC)  for  automatic  processing.  The  ECAC  data  have  been 
reorganized  slightly  to  fit  better  the  constraints  of  the  CDC  6000  series 
computer. 

The  basic  unit  of  data  is  a "bloclt"  covering  an  area  1°  by  1°  with 
corners  at  integer  values  of  latitude  and  longitude.  It  is  referenced 
by  its  southwest  corner.  Within  it,  the  data  can  best  be  pictured  as  arranged 
in  a two-dimensional  array  with  each  index  ranging  from  0 to  120.  If  Z 
represents  such  an  array,  then  Z(IY,IX)  is  the  ground  elevation  at  latitude 
LAT  + IY/120,  longitude  LON  + IX/120,  where  LAT,  LON,  are  the  coordinates 
of  the  southwest  corner.  Coordinates  are  expressed  in  degrees  and  decimal 
parts  of  a degree,  and  latitudes  are  expressed  as  negative  in  the  southern 
hemisphere.  Longitudes  west  of  the  reference  meridian  also  are  negative. 

Note  that  within  CONUS,  only  positive  latitudes  and  negative  longitudes 
are  encountered. 

The  array  size  for  one  block  is  121  by  121  with  a total  of  14,641 
elements.  The  data  for  each  element  are  packe'd  into  60-bit  words  according 
to  a format  that  produces  variable  length  data  records.  The  elevations 
are  expressed  as  coded  integer  multiples  of  20  ft  (6.1  m).  The  minimum 
elevation  value  within  each  data  block  is  decreased  by  20  ft  (6.1  m),  and 
that  value  becomes  a biased  reference  elevation  for  that  block.  The  biased 
reference  elevation  is  subtracted  from  each  elevation  datum  within  the 
block,  ^tnd  that  elevation  difference  then  is  expressed  as  a positive  integer 
multiple  of  20  ft  (6.1  m).  (The  zero  value  is  used  to  indicate  missing 
data.  Next,  it  is  determined  whether  2,  3,  4,  5,  6,  7,  8,  10,  or  12  bits 
are  required  to  represent  the  largest  of  these  elevation  differences. 

Finally,  the  reduced  data  are  packed  consecutively  from  left  to  right  in 
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words  30,  20,  15,  12,  10,  8,  7,  or  6 elements  per  word,  respectively. 

(In  the  case  of  8-  or  7-bit  elements,  the  four  low  order  bits  of  each  word 
are  unused  and  set  to  zero  for  that  reason.)  The  topographic  data  tape 
organization  and  file  structure  are  shown  in  Figure  E-1 . Figure  E-2  shows 
how  the  data  tape  organization  relates  to  the  geographical  area. 

A 2-bit  element  allows  a total  elevation  difference  range  of  only 
40  ft  (12.2  m) . However,  if  each  coded  elevation  difference  integer  within 
a data  block  could  be  expressed  in  this  way,  only  489  words  would  be  required 
for  the  coded  and  packed  elevation  difference  data.  At  the  other  extreme, 
if  the  elevation  difference  were  such  that  the  coded  integers  required 
12-bits  per  datum,  such  a code  would  allow  an  elevation  difference  of  81,880 
ft  (24,960  m)  and  2929  words  would  be  required  to  express  the  coded  and 
packed  elevation  difference  data  within  the  block.  A 10-bit  element  will 
allow  an  elevation  difference  of  20,440  ft  (6230  m)  which  should  be  sufficient 
for  any  data  block  within  the  CONUS.  When  10-bits  per  datum  are  used, 

2441  words  are  required  for  the  data  of  that  block. 

Consider  an  example  of  decoded  terrain  elevation  where,  for  a particu- 
lar block,  the  minimum  elevation  is  2360  ft  (719.3  m) , the  maximum  elevation 
is  8740  ft  (2664  m) , and  the  elevation  difference  code  is  247.  (One  might 
note  that  10  bits  per  datum  would  be  required  to  code  the  elevation  differ- 
ence data.)  The  biased  reference  elevation  would  be  2340  ft  (713.2  m) 
for  the  block.  Elevation  at  the  point  would  be  247  x 20  + 2340  = 7280  ft 
(2219  m) . 

The  alphanumeric  descriptor  contained  in  word  six  of  the  coded  elevation 
difference  record  (illustrated  in  Figure  E-1)  nearly  always  will  show  a 
letter  "A"  which  indicates  the  data  source  to  be  the  DMA  or  its  predecessor 
agency,  the  Army  Map  Service. 


Partition 

1 

Data 

for 

23°N 

latitude 

Partition 

2 

Data 

for 

24°N 

latitude 

Partition 

28 

Data 

for 

50°N 

latitude 

Record  1 
Record  2 
Record  3 


Record  71 


Latitude  and  index  to  data  blocks. 

Coded  elevation  differences  for  block  70  (see  Figure  F-2) 
Coded  elevation  differences  for  block  69  (see  Figure  F-2) 


Coded  elevation  differences  for  block  1 (see  Figure  F-2) . 


(Fewer  records  will  be  used  when  data  for  a block  are  missing  or  the  block 
is  ocean.) 


Word  1 
Word  2 
Word  3 


Word  71 


Integer  value  of  latitude  for  that  partition. 
♦Index  for  block  70. 

♦Index  for  block  69. 


♦Index  for  block  1. 


♦Index  value  is  zero  if  data  are  missing  for  the  block. 
Index  value  is  ~1  if  block  is  ocean.  A positive  integer 
declares  the  numt>er  of  the  record  containing  data  for 
that  block. 


Figure  E-1.  Organization  of  topographic  data  tape. 

153 


Coded  Elevation  Differences  Record  Structure 


Word  1 Latitude  and  Longitude  of  SW  corner  of  the  block. 
Bits  30-59:  Latitude  (North  = +) 

Bits  0-29:  Longitude  (East  = +) 

Word  2 Bits  30-59:  Number  of  missing  data  points 
Bits  12-29:  Mask 

Bits  6-11:  Number  of  bits  per  datum 
Bits  0-5:  Number  of  data  points  per  word 
Word  3 Integer  value  of  biased  reference  elevation  in  feet. 

Word  4 Integer  value  of  minimum  elevation  in  feet. 

Word  5 Integer  value  of  maximum  elevation  in  feet. 

Word  6 Alphanumeric  descriptor  of  data  source. 

Word  7 Packed  and  coded  elevation  differences,  6 to  30 

points  per  word,  ordered  south  to  north  and  west  to 
east  starting  in  the  SW  corner  of  the  block. 


(Sufficient  words  for  14,641  points) 


Figure  E-1.  Continued 
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Figure  E-2.  Illustrated  relationships  between  geographical  area  and 
data  tape  org2mization  for  CONUS. 


APPENDIX  F.  DISCUSSION  OF  OUTPUT  DATA  AND 
CONTOURING  METHODOLOGY 

The  variety  of  calculated  output  data  (basic  transmission  loss,  power 
density,  received  signal  level,  and  communication  reliability)  provided 
by  the  CSPM  allow  the  CSPM  to  be  applied  to  a communication  system  problem 
at  nearly  any  stage  of  the  development  cycle  for  the  system.  For  example, 
basic  transmission  loss  can  be  calculated  when  only  operating  frequency, 
antenna  heights,  and  polarization  are  known.  When  the  transmitter  output 
power  and  transmitting  antenna  gain  are  known,  the  model  can  provide  power 
density.  Additional  detail  about  the  system  and  its  application,  such 
as  the  gain  and  height  of  each  antenna,  polarization,  and  the  transmitter 
output  power,  will  allow  calculations  of  received  signal  level.  Finally, 
knowledge  of  or  cui  assumption  about  received  signal  level  threshold  for 
satisfactory  performance  allows  an  output  termed  "communication  reliability" 
which  provides  em  overall  estimate  of  satisfactory  performance  as  a function 
of  given  time  availability  and  path  difference  conditions. 

From  Appendix  C,  we  know  the  basic  transmission  loss  that  is  not  exceeded 
a fraction  of  the  time  at  a fraction  q^  of  the  locations  with  a prediction 
confidence  Q is 

L^(q„*q. rQ)  = L-  + A - V - - Y^  - Y^,  (F-1) 

D ^T  fs  cr  T L C 


where  L,  is  the  calculated  free  space  loss, 
rs 

A is  the  calculated  attenuation  factor  relative  to  free-space 
cr 

attenuation, 

V is  eui  en^irically  determined  adjustment  factor  to  the  reference 
attenuation  factor  for  specific  seasons  and  climates,  and 

Y , Y , and  Y are  the  attenuation  variabilities  imposed  upon  median 

X Xi  Ck 

basic  transmission  loss  to  account  for  long-term  fading  effects, 
random  path  differences,  emd  prediction  confidence,  respectively. 
It  will  be  useful  later  in  this  appendix  to  recognize  that 

expresses  the  median  basic  transmission  loss. 
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The  power  per  unit  area  incident  upon  a surface  at  a distance  R from 


the  transmitting  antenna,  termed  power  density,  is 


= P^  g/(4TTR  ) 


(F-3) 


where  P is  the  power  density  in  W/m  , 
d 

P^  is  the  transmitter  output  power  in  watts, 

g is  the  gain  of  the  transmitting  antenna  (g=l  for  isotropic 
radiator) , and 

R is  the  distance  from  the  transmitting  antenna  in  meters. 

In  the  CSPM,  it  is  convenient  to  consider  relative  values  of  power  density, 

2 

expressed  as  decibels  relative  to  a milliwatt  per  square  meter  (dBm/m  ) , 
and  to  use  the  statistical  power  density  values.  The  frequency  dependence 
of  basic  transmission  loss  must  be  removed  from  the  value  when  following 
such  a scheme.  The  expression  for  power  density  in  decibels  relative  to 
a milliwatt  per  square  meter  then  can  be  written  as 


P^  = 20  log^pf  ^T  ^1  ■ 


{F-4) 


where  f is  the  carrier  frequency  in  megahertz, 

P^  is  the  transmitter  output  power  in  decibels  relative  to  a milliwatt, 

and 

is  the  gain  of  the  transmitting  antenna  in  decibels  relative  to 
an  isotropic  antenna  (dBi) . 

The  received  signal  level,  RSL,  expresses  the  power  delivered  by  a receiving 
antenna  to  the  input  of  its  receiver.  This  quantity  takes  into  consideration 
the  statistical  variations  in  signal  attenuation  and  is  calculated  as 

RSL(q^,q^Q)  = S ^2  " ^pg  " 

where  P^  and  are  as  defined  for  (F-4) , 

G^  is  the  gain  of  the  receiving  antenna  in  decibels  relative  to  an 
isotropic  antenna  (dBi) , and 

L is  the  antenna  aperture- to-medium  coupling  loss  in  decibels. 

pg 

= 0 dB,  for  ^2,  ^ ^2  ^ dBi. 

L = 1 dB,  for  25  < G,  + i.  50  dBi. 

pg  12 

= 0.07  exp  0.055  (G^+G^)  dB,  for  50  < G^^  + G^  < 100  dBi. 
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since  (F-1)  expresses  the  basic  transmission  loss  that  is  not  exceeded, 

as  a function  of  q , q , and  Q,  (F-5)  is  an  expression  of  the  received 
7 L 

signal  level  that  is  exceeded  q^  fraction  of  the  time  at  q^  fraction  of 
the  locations  with  prediction  confidence  Q. 

One  can  define  satisfactory  system  performance  to  be  provided  when 

RSL(q^,q^,Q)  > RSL^^  or  RSl,(q^,q^,Q)  - RSL^  > 0,  (F-6) 

where  is  a received  signal  level  threshold  for  satisfactory  performance 

of  the  system  being  considered.  Subtracting  RSL  from  each  side  of  (F-5) 

TH 

and  substituting  (F-1)  and  (F-3) , 

■ '■t  * =1  + =2  - 

As  a mathematical  convenience,  let  S.  = P + G,  + - L + V„  + Y,  - RSL„. . 

0T12pgTL  TH 

Then  (F-7)  becomes 


- MI-™  ■ =0  - ^0  * ''C-  "’-®’ 

Apply  the  conditions  of  (F-6)  by  equating  the  right-hand  side  of  (F-8) 
to  zero,  and  utilize  the  expression  for  given  by  (C-11)  in  Appendix  C 
to  obtain 

Vo'=’  ■ =0  - h,0- 

Now  use  (C-7b)  in  Appendix  C to  solve  for  Q as 


s . 0.5  . 0.5  erf  [(S„-L^„)/f7  oj. 


(F-10) 


Equation  (F-10)  expresses  the  probability  that  received  signal  level,  a 
function  of  q^  and  q^,  will  exceed  the  threshold  for  received  signal  level. 
That  is,  (F-10)  expresses  the  communication  reliability  as  a probability 
that  received  signal  level  will  exceed  the  received  signal  level  threshold 
a specified  fraction  q^  of  the  time  at  a specified  fraction  q^  of  all  possible 
locations . 

Any  area  within  the  limits  of  the  topographic  base,  up  to  a maximum 
size  of  2°  latitude  by  2°  longitude,  may  be  specified  for  plotting  equal- 
value  contours  of  basic  transmission  loss,  power  density,  received  signal 
level,  and  communication  reliability  about  a specified  point  within  the 
area.  Appendix  C discusses  the  Longley-Rice  propagation  loss  model  which 
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is  used  in  the  point-to-point  mode  to  statistically  estimate  signal  attenuation 
to  many  points  within  the  area.  These  estimates  of  signal  attenuation 
are  used  to  calculate  the  data  which  can  be  contoured. 

It  is  of  interest  to  discuss  the  algorithm  used  to  select  the  points, 
each  used  with  the  specified  point  to  form  the  many  point-to-point  paths 
for  signal  attenuation  calibrations.  Only  two  basic  schemes  seem  likely 
as  candidates.  One  scheme  we  shall  call  the  "grid  method";  the  other  scheme 
is  the  "radial  method."  Simply  stated,  the  grid  method  requires  the  geographic 
area  of  interest  to  be  covered  by  grid  lines  with  appropriate  separations, 
and  each  grid  line  intersection  defines  a point  which  with  the  transmitter 
location  forms  a propagation  path.  The  principal  advantages  of  the  grid 
method  are  the  ease  of  defining  the  points  and  the  uniform  distribution 
of  the  points.  The  principal  disadvcintage  of  such  a scheme  is  that,  in 
general,  a unique  path  profile  must  be  generated  for  the  path  to  each  point. 
With  1 km  separation  between  points,  more  than  45,000  path  profiles  could 
be  required.  In  simple  terms,  the  radial  method  requires  the  selection 
of  a number  of  radials,  with  uniform  azimuth  angle  separation,  extending 
from  the  transmitter  location  to  intersection  with  the  area  boundary. 

Points  then  are  selected  at  fixed  intervals  along  each  radial  to  form  propaga- 
tion paths  with  the  origin.  The  disadvantage  of  this  method  is  that  uniform 
distribution  of  the  points  is  impossible,  although  an  algorithm  will  be 
discussed  which  achieves  approximately  uniform  distribution.  A very  attrac- 
tive feature  of  the  radial  method  is  the  necessity  to  generate  path  profiles 
only  for  each  radial,  thereby  utilizing  a single  (radial)  path  profile 
for  m£uiy  propagation  paths  along  that  radial.  An  innovative  scheme  which 
provides  approximately  uniform  distribution  of  points  involves  disregarding 
points  near  the  origin  along  certain  radials.  Within  a 45°  sector,  each 
time  the  azimuth  angle  increment  between  radials  is  divided  in  half,  the 
distcuice  from  the  origin  is  doubled  at  which  points  along  the  radial  are 
used  first  for  establishing  propagation  paths  with  the  origin.  For  radials 
at  0°,  45°,  90°,  etc.,  points  separated  by  1 km  starting  at  a point  1 km 
from  the  origin  are  used.  For  radials  at  22.5°,  67.5°,  112.5°,  etc.,  points 

separated  by  1 km  starting  at  2 km  from  the  origin  are  used.  For  radials 
o o 

at  11.25  , 33.75  , 56.25  , etc.,  poihts  separated  by  1 km  starting  at  4 km 
from  the  origin  are  used.  Table  F-1  shows  the  relationships  between  azimuth 
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angle  increment,  distance  from  origin  at  which  new  radials  begin,  and  total 
number  of  required  radials  in  360°.  For  example,  if  a particular  analysis 
involves  maximum  radial  distances  of  100  km.  Table  F-1  shows  that  only 
512  radials  will  be  required,  angular  separation  between  radials  will  be 
0.70°,  and  propagation  paths  will  begin  at  distances  of  1,  2,  4,  8,  16, 

32,  and  64  km  from  the  origin.  Figure  F-1  shows  the  starting  point  methodology. 

Table  F-1.  Tabular  Information  on  the  Radial  Method 

for  Developing  Data  to  be  Plotted  as  j 

Equal-Value  Contours 


Radial  Separation 
Increment,  Degrees 

Minimum  Path 
Length,  (2*^) , km 

n 

Total  Number  of 
Radials  in  360 

90 

1 (defined) 

- 

4 

45 

1 (defined) 

0 

8 

22.5 

2 

1 

16 

11.25 

4 

2 

32 

5.63 

8 

3 

64 

2.81 

16 

4 

128 

1.41 

32 

5 

256 

0.70 

64 

6 

512 

0.35 

128 

7 

1024 

0.18 

256 

8 

2048 

Figure  F-1  is  a pictorial  presentation  for  a 45°  sector  of  essentially 
the  same  information  given  in  Table  F-1.  Note  from  the  figure  that  the 
separation  between  data  points  along  each  radial  is  always  1 km,  and  from 
radial  to  radial  the  separation  varies  from  0.79  to  1.57  km.  The  restriction 
of  2°  latitude  by  2°  longitude  for  the  geographical  area  limits  maximum 
radial  length  that  must  be  considered.  If  the  origin  were  located  at  one 
comer  of  the  area,  the  radial  extending  to  the  opposite  comer  would  be 
only  slightly  more  than  300  km  maximum  in  length. 

In  Appendix  D,  it  is  stated  that  the  interval  for  path  profile  points 
is  fixed  at  250  m for  radial  paths.  That  interval  is  a convenient  modulus 
for  the  1 km  Interval  at  which  basic  transmission  loss  data  are  calculated, 

as  discussed  in  this  appendix.  Yet  one  quickly  realizes  that,  in  general, 
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Angular  Separation  Between  Radials,  Degrees 


Figure  F-1.  Pictorial  presentation  of  minimum  path  lengths,  angular  radial 
separations,  and  point  separations  for  radial  method. 


the  radial  path  length  resulting  from  the  intersection  of  the  radial  with 
the  area  boundary  will  not  be  integer.  This  problem  is  solved  by  rounding 
up  to  the  next  integer  value  (in  kilometers) , computing  the  geographic 
coordinates  for  the  radial  path  endpoint  beyond  the  area  boundary,  and  using 
that  path  in  the  generation  of  the  radial  path  profile  and  calculations 
of  data  to  be  contoured.  The  contouring  subroutine  simply  is  restricted 
to  plotting  within  the  area  boundaries. 

Output  plots  of  basic  transmission  loss,  power  density,  and  received 

signal  level  may  contain  up  to  ten  contoured  values.  Program  logic 

determines  the  number  of  and  the  values  for  the  contour  levels  in  the 

following  way.  The  range  in  values  throughout  the  geographical  area  for 

a particular  parameter,  i.e.,  basic  transmission  loss  for  q = .95,  q = .95, 

T L 

and  Q = .50,  is  determined.  That  range  is  divided  by  5,  10,  15,...  to 
obtain  the  number  of  subintervals,  limited  to  ten  maximum  with  smallest 
subinterval  value.  Subinterval  values  are  rounded  up  to  values  divisible 
by  five  and  these  become  the  values  at  which  contours  are  plotted.  The 
CSPM  user  specifies  the  values  up  to  four  possible,  at  which  contours  of 
communication  relicibility  are  to  be  plotted. 
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